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Chapter 1 



Introduction 



Quantum Chromodynamics, as a gauge theory of strong interaction, is one 
of the major achievements of particle physics in the last 50 years. However, 
at sixties of the last century, particle physics was not compelling enough. 
Quantum Electrodynamics produced increasingly predictions confirmed ex- 
perimentally. Few years later, with the electroweak gauge theory of the 
Standard Model, also the weak interactions were included in the well de- 
fined framework of Quantum Field Theories. At that time, among all the 
fundamental interactions, the strong interactions still remained outside that 
description. With the invention of Bubble and Spark chambers, experimen- 
tal physics discovered a large number or particles called hadrons. At sixties, 
a glance at the hadron spectroscopy show that the increasing number of 
hadrons could be classified successfully by the quark model of Gell-Mann 
and Zweig. However, the dynamics behind the quark model was still a mys- 
tery. The interpretation of strong interactions in the framework of Quantum 
Field Theories was useful but only at the level of a toy model and Gell-Mann 
himself, who suggest this approach, declared the quarks to be purely math- 
ematical entities without any physical meaning. Hard work on the exper- 
imental side was needed to finally conclude with the acceptance of quarks 
been elementary particles which, encoded in a SU{3)c gauge symmetry, end 
up in the so-called QCD Lagrangian yielding the proper description of the 
strong interactions. 

The QCD Lagrangian is described then in terms of quarks and glu- 
ons which, due to confinement, can not be detected in isolation, but only 
colorless combinations of them, the hadrons. At high energies, however, 
asymptotic freedom renders the strong coupling to be small enough that 
perturbation theory becomes the appropriate tool to work with. At en- 
ergies below 1 GeV, confinement enters in the game binding quarks and 



gluons together and invalidating the perturbative technics. Models like the 
Constituent Quark Model, numerical simulations such as lattice QCD or 
analytical approximations to QCD are among the non-perturbative technics 
used nowadays to explore the non-perturbative regime. 

Also Effective Field theories have successfully contribute to a better un- 
derstanding of the properties of QCD, specially in its low-energy regime. 
In particular Chiral Perturbation Theory as a systematical organization of 
the physics in powers of momenta and quark masses resulted in a good tool 
for describing the low-energy phenomenology, but requires the knowledge of 
several low-energy constants to improve on the results. 

On the other hand, the large-A'^c expansion standed out as a promis- 
ing analytical approach capable of dealing with the complexities of non- 
perturbative QCD regime while offering a simple description of the physics. 
However, a limitation to doing phenomenology was found due to the fact 
that there is still no solution for large- A'^c QCD. On the other hand, the phe- 
nomenological approach of resonance saturation works quite well for QCD 
Green's functions. It was then realized that this successfully result could 
be encompassed at once as an approximation to large- A^c QCD consisting 
in keeping only a finite (instead of the original infinite in large- A^c) set of 
resonances in Green's functions, the Minimal Hadronic Approximation. 

The present work is devoted in a first part to study how can we under- 
stand the Minimal Hadronic Approximation an the limitations that appear 
when dealing with high-energy matching conditions with a finite number of 
resonances and how reliable is to extract information on individual mesons 
from it. We suggest then that resonance saturation in QCD can be under- 
stood in the large- A^c limit from the mathematical theory of Pade Approxi- 
mants to meromorphic functions. Due to the success of Pade Theory in the 
framework of large- A^c QCD thanks to the convergence theorems described 
in the literature, we also apply this technic in a second part to other open 
questions, such as what is the insight behind the unitarization processes 
at low energies, how to work out the possible ambiguities appearing in the 
approximations to the vacumm polarization of a heavy quark and, finally, 
how to go on improving the Vector Meson Dominance when dealing with 
experimental data. 

The outline of the present work will be as follows: in chapter [2] we give 
a brief reminder of the basics of Quantum Chromodynamics together with 
its low energy chiral realization (Chiral Perturbation Theory) and its high 
energy expansion through the Operator Product Expansion. We then intro- 
duce the 1/Nc expansion technic to deal with extensions of the chiral La- 
grangian to intermediate energies with the explicit introduction of resonance 



fields (what is called Resonance Chiral Lagrangian). Since this lagrangian 
needs an infinite number of resonances and in practice this request can not 
we easily attended, we also show the approximations that has been used in 
the past years, mainly the Minimal Hadronic Approximation, which is based 
in the phenomenological success of resonance saturation. 

Chapter [3] introduces the notion of Fade Approximant in a brief but 
precise mathematical way, with several examples to illustrate its usefulness 
and, more important, all the theorems of convergence that will we needed 
to tackle all the raised question. 

Chapter H] is focused on the application of Fade Theory to meromorphic 
functions since in the framework of large- A'^c QCD it turns out that reso- 
nance saturation are rational functions which encompass any saturation with 
a finite number of resonances. In particular the consequences of the Fom- 
merenke's theorem of convergence applied to rational functions are studied 
and it is shown that while this rational approximants may reliably describe 
a Green's functions in the Euclidean, the same is not true for the Minkowski 
regime. As an example, we estimate the value of low-energy constants for 
the {VV — AA) two-point correlator. 

Chapter [5] is devoted to study the reliability of the Fade Theory applied 
to Stieltjes functions when unitarizing low energy amplitudes or when ex- 
tracting threshold parameters from the vacuum polarization function of a 
heavy quark. 

In chapter [6] we illustrate how Fade Approximants can go beyond the 
vector meson dominance approach with the case of the pion vector form 
factor in the spacelike region and estimate the constants that appear in its 
low-energy realization. It is shown how Fade Approximants are useful to 
incorporate information from high energies in a straightforward way. 

Finally, we end with the conclusions. 



Chapter 2 



Quantum Chromodynamics 



The Standard Model of Particle Physics (SM) is an extremely successful 
theory of physics at the subnuclear scalqj. Starting from three generations 
of quarks and leptons subject to a local SU{3)c x SU{2)l x f7(l)y symmetry, 
it describes all observed phenomena known today (without including gravity 
effects) within theoretical and experimental errors. The first factor SU{2>)c 
accounts for the strong interactions, whose discussion will be the main sub- 
ject of this chapter. The remaining SU(2)l x U{1)y factor is responsible 
for the electroweak interactions. 

Our understanding of the electroweak and strong interactions is good 
and is based on a well defined mathematical theory. Actually, in the elec- 
troweak sector (for example Quantum ElectroDynamics, QED) successful 
mathematical tools are used, i.e.. Perturbation Theory (PT) as an expansion 
of the coupling constant. However, the strong sector (Quantum ChromoDy- 
namics) has a regime where PT can be applied and another regime where 
it can not. The reason is because the PT breaks down in certain regime 
called non-perturbative Quantum ChromoDynamics. Therefore, some ef- 
forts have been made towards the development of new tools that allow reli- 
able computations in this non-perturbative regime, for example the lattice 
gauge theories or analytical approximations. In the present chapter, we will 
present the main properties of the strong interactions, while focuss in both 
low-energy and high-energy regimes. We will also present several attempts 
done in the literature to try to link both regimes in a systematical way. In 
that sense, we will review beyond Quantum ChromoDynamics, Chiral Per- 
turbation Theory, the Operator Product Expansion technics, the Large-A'^c 



""^We will not present a detailed discussion about Standard Model, but just a brief 
overview of the main points, focused on the strong interaction sector. For comprehensive 
and general accounts, we refer the reader to [T] [2] . 



limit in Quantum ChromoDynamics and the Resonance Chiral Theory. 

2.1 Quantum Chromodynamics: Introduction 

Quantum Chromodynamics (QCDJj is nowadays the estabhshed theory 
of the strong interactions, a fundamental force describing the interactions of 
the quarks and the gluons. Mathematically, is a SU{3) Yang-Mills theory of 
color-charged entities. Color, in this framework, is a new quantum number 
corresponding to a SU{3)c symmetry, c standing for color. We would like to 
present some evidence of its existence before writing down the corresponding 
Lagrangian. 

The Particle Data Tables [5] reveals the richness and variety of the 
hadronic spectrum. In the sixties, the rapidly increasing number of hadrons 
discovered could be classified successfully by the quark model (with the three 
quarks u, d and s) [6], showing the existence of a deeper level of elementary 
constituents of matter, the quarks. Assuming that mesons are qq states 
(quark-antiquark) , while baryons are qqq states, the quark model allow us 
to classify the entire hadronic spectrum as a Periodic Table of Hadrons. 
However, the quark picture faces a basic problem suggesting the presence 
of a new degree of freedom: a fundamental state of composite system is ex- 
pected to have angular momentum zero, a baryon (A"'""'") with total angular 
momenta J = 3/2 corresponds to uuu with the three quark-spins aligned 
into the same direction and with the relative angular momenta equal to zero. 
In that case, the wave function is symmetric and our state obeys a wrong 
Fermi-Dirac statistics. A new degree of freedom, a new quantum number 
for quarks, the color, can solves this problem because our state can be to- 
tally antisymmetric in these color indices (at least 3 colors are needed for 
making an antisymmetric state) respecting the generalized Pauli principle 
with spatially symmetric wave function. 

However, assuming that the number of colors Nc = 3 for each quark and 
antiquark, we obtain 9 x 9 = 81 combination of qq only nine of which had 
been found. In order to avoid the existence of non-observed extra states with 
non-zero color, one needs to further postulate that all asymptotic state are 
colorless, singlets under rotations in color space. This postulate, called the 
confinement hypothesis, implies the non-observability of free quarks since 
they carry color. 

The decisive clue to manifest the existence of the color of the quarks 
came from experiment. Started by the MIT-SLAC collaboration at the end 



^For a pedagogical reviews see for example [3], [1] 



of the sixties, deep inelastic scattering of leptons on nucleons and nuclei 
produced unexpected results. Whereas at low energies the cross sections 
were characterized by baryon resonance production, the behavior at large 
energies and momentum transfer suggested that the nucleons seemed to 
consist of noninteracting partons. Obvious candidates for the partons were 
the quarks but this idea led to a seeming paradox. How could the quarks 
be quasi-free at high energies and yet be permanently bound in hadrons, a 
low-energy manifestation? 

That the strength of an interaction could be energy dependent was not 
really new to theorists. In Quantum ElectroDynamics, the vacuum acts 
like a polarisable medium leading to the phenomenon of charge screening, 
but the effective charge increases with energy. In the case of QCD, the 
deep inelastic experiments seemed to suggest the clue to understand that 
phenomena. 

The deep inelastic scattering measurements actually requires the exis- 
tence of electrically neutral as well as charged constituents of the proton. 
Charged partons can be identified with the colored quarks and the neutral 
partons with gluons. These two identifications are needed to have direct 
evidence that quantum chromodynamics is the correct physical theory of 
strong interactions. Its essential properties are: 

• Quarks carry color as well as electric charge; there are three colors: 
Red, Green and Blue. 

• Color is exchanged by eight bicolored gluons, which are massless and 
have spin 1. 



• 



Color interactions are assumed to be like the electromagnetic once (in 
a quark-gluon interaction use the rules of QED with the substitution 

\/a -^ \/aI)- 

• Gluons themselves carry color charge, and so they can interact with 
other gluons. 

• At short distances (large— Q^), as is sufficiently small so that we can 
compute color interactions using the perturbative techniques. 

The strong interactions, then, must have the property of asymptotic 
freedom (Groos and Wilczek [7], Politzer 1973 [8]): the interactions in the 
strong sector should become weaker at short distances, so that quarks behave 
as free particles for Q^ — )• oo. This also agrees with the empirical observation 
that the free-quark description of the ratio 



_ a{e~^e — )• hadrons) 
cr(e+e — > ^+/i j 

works better at higher energies. Thus, the interaction between a qq pair 
looks hke some kind of rubber band. If we try to separate the quark form 
the antiquark the force joining them increases. At some point, the energy 
on the elastic band is larger than 2mq', so that it becomes energetically 
favorable to create an additional q'q' pair; then the band breaks down into 
two mesonic systems, qq' and q'q, each one with its corresponding half- 
band joining the quark pair. Increasing more and more the energy, we can 
only produce more and more mesons, but quarks remain always confined 
within color-singlet bound states. Conversely, if one tries to approximate 
two quark constituents into a very short-distance region, the elastic band 
loses the energy and becomes very soft; quarks behave then as free particles. 
The ratio in Eq. ()2.ip is then given by the sum of the quark electric charges 
squared: 



Nf 



/=1 



iNc = 2 , (iV/ = 3 : n, d, s) 

fiVc = f , {Nf=A : u,d,s,c) (2.2) 

^-Nc = ^, iNf=5 : u,d,s,c,b) 



Notice that strong interactions have not been taken into account, only 
the confinement hypothesis has been used. 

The hadronic decay of the r lepton provides additional evidence for 
Nc = 3. The branching rations for the different channels are expected 
to be approximately: 

Br^iepton = Br{T' -> ^l' i>i) ~ 2 + N ^ 5 ^^■^'' 

for Nc = 'i which should be compared with the experimental average [S], for 
example. B^^^ = (18.01 ± 0.18)%. 

Finally, we would like to comment also on the anomalies as another 
compelling reason to adopt N^- An anomaly is a global symmetry which is 
broken by quantum effects. In particular, the theoretical prediction for the 
decay vr'' — )• 77, which occurs through a triangular quark loop related by the 
anomaly, needs A'^c = 3 to predict the experimental value. Also, anomalous 
triangle diagrams involving electroweak gauge bosons could spoil renormal- 
ization of the Standard Model. This potentially dangerous contributions 
are canceled due to a subtle combined action of quarks and leptons in each 
of the three generations, provided that quarks appear with three different 
colors, Ac = 3. 



2.1.1 QCD Lagrangian 

Strong interactions, then, bind quarks inside the atomic nuclei due to the 
mediation of gluons, the gauge bosons of the color group. The interaction 
is described by a Yang-Mills quantum field theory under the non-abelian 
gauge group SU{3)c of color. Taking the proper normalization for the gluon 
kinetic term, we finally have the invariant QCD Lagrangian: 

Cqcd ^-\ CrC^. + Yl ^f (^^''^M - mf) qf . (2.4) 

/ 

We required the Lagrangian Eq. ()2.4p to be also invariant under local 
SU{3)c transformations, 6a = Oa{x). To satisfy this requirement, we needed 
to change the quark derivatives by covariant objects. Since we have now 8 
independent gauge parameters, 8 different gauge bosons Ga{x), also the 
gluons are needed: 



D% 



d'^-igs^G^'^i^) 



qf ^ [d^^-igsG^^{x)]qf. (2.5) 



Notice that we have introduced the compact matrix notation 

where A*^ (a = 1,2, ...,8) denote the generators of the fundamental 
representation of the SU{3)c algebra (the Gell-Mann matrices A"). The 
matrices A'^ are traceless and satisfy the commutation relations 



A^A* 



2ir^\\ (2.7) 



with f the SU{3)c structure constants, which are real and totally 
antisymmetric. 

The gauge transformation of the gluon fields is more complicated than 
the one obtained in QED for the photon. The first part of Eq. (j2.4p and the 
non-commutativity of the SU{3)c matrices gives rise to an additional term 
involving the gluon fields themselves. The gluon fields belong to the adjoint 
representation of the color group. Note also that there is a unique SU{3)c 
coupling gs, which is called the strong coupling constant. The existence of 
self-interactions among the gauge fields is a new feature that was not present 
in QED; it seems then reasonable to expect that these gauge self- interactions 
could explain properties like asymptotic freedom and confinement, which do 
not appear in QED either. 



The asymptotic freedom, due to the fact that the strong couphng is 
small at high energies, let quarks be unbounded and behave as free parti- 
cles. In that regime, the parton model is explained and perturbation theory 
techniques are allowed. 

Confinement, however, is found at low energies where quarks become 
more and more tightly bounded due to the increasing of the value of the 
strong coupling constant, entering in a non-perturbative regime. This is 
also known as infrared slavery. Confinement makes the quark-gluon picture 
transform to the hadronic picture we observe in particle accelerators. 

The picture is then quite intricate: with quarks and gluons we are able 
to write down the QCD Lagrangian Eq. ()2.4p but we do not know how 
to evolve it to its dual description in terms of the asymptotic hadronic 
states. Several technics were developed to tackle this problem under effective 
field theories and symmetry arguments. Among them, Chiral Perturbation 
Theory provides the most general framework to deal with asymptotic states 
at low energies, and will be the remaining of the next section. 

2.2 Chiral Perturbation Theory 

Chiral Perturbation Theory (ChPTJfj provides a systematic framework for 
investigating strong-interaction processes at low energies, as opposed to a 
perturbative treatment of QCD at high momentum transfer in terms of the 
"running coupling constant". The basis of ChPT is the global SU{3)l (^ 
SU(3)Ft ^ U{l)v symmetry of the QCD Lagrangian in the limit of massless 
u, d, and s quarks. This symmetry is assumed to be spontaneously broken 
down to SU{3)v (E> U{l)v giving rise to eight massless Goldstone bosons. 
Actually, below the resonance region (E < Mp), the hadronic spectrum 
only contains an octet of very light pseudoscalar particles (vr, -ff, 77), whose 
interactions can be easily understood using this symmetry consideration. 

The ChPT formalism is then based on two key ingredients: the chiral 
symmetry properties of QCD and the concept of Effective Field Theory 

(EFT) [ma US]. 

Effective field theories are the appropriate theoretical tools to take ex- 
plicitly into account the relevant degrees of freedom for the process under 
investigation, i.e., those states with m « A, while the heavier excitations 
with M >> A are integrated out from the action. The information on the 
heavier degrees of freedom is then contained in the couplings of the result- 
ing low-energy Lagrangian, the Low-Energy Constants (LECs). Although 



^For a pedagogical review see for example [9], [TO], [TT], [12) . [13) 
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effective field theories contain an infinite number of terms at a given order in 
the energy expansion, the low-energy theory is specified by a finite number 
of couplings; this allows for an order-by-order renormalization [17j . 

On the chiral symmetry side, the other ingredient of Chiral Perturbation 
Theory, it turns out that in the limit of massless n, d and s quarks the 
QCD Lagrangian Eq. (|2.4p has a global symmetry spontaneously broken. 
Actually, 

^QCD = -^G^.Gr + iqil^D^qL + iqRl^D^qR (2.8) 

and is invariant under the global G = SU {N f) l<^ SU {N f) r®U {l)v (^U {I) a 
transformations of the left and right handed quarks in flavour space: 

qi -^ gLqi: qn -> gnqn, gL,R e SU{Nf)L,R . (2.9) 

The two quark chiralities live in separate flavour spaces implying that all 
previous flavor symmetries get duplicated in two chiral sectors {SU{Nf)L 
and SU{Nf)R spaces). On the other hand, under U{l)v all quarks have the 
same change in phase and corresponds to the baryon number while under 
U{1)a the right and left-handed quarks have the opposite change in phase. 
The symmetry is called chiral because it acts differently on the left and 
right-handed quarks. The U{1)a is only a symmetry of the classical action, 
not of the full quantum theory of QCD, it is an anomaly and, therefore, the 
divergence of the associated current does not vanish. It is nonzero by a total 
divergence but instantons allow for this to have a physical effect. In sum, 
the final symmetry of QCD in the limit where all quarks are massless, and 
taken in to account only the meson sector, is thus 

G^ = SU{3)l SU{3)r (2.10) 

The Noether currents associated with the chiral group Gy. are 

JT = Qx1^Y'^x, {X = L,R; a=l,...,8). (2.11) 

where Aq are the Gell-Mann's matrices with Tr{XaXb) = 26ab- The cor- 
responding Noehter charges Q^ = J d^xJ^{x) satisfy the familiar commu- 
tation relations 

m,Q'Y]=idxYfabcQx (2.12) 

which were the starting point of the Current Algebra methods of the six- 
ties. This chiral symmetry, which should be approximately good in the 
light quark sector (u,d,s), is however not seen in the hadronic spectrum. 
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Although hadrons can be classified in SU{3)v representations, degenerate 
multiplets with opposite parity do not exist. Moreover, the octet of pseud- 
scalar mesons happens to be much lighter than all the other hadronic states. 
To be consistent with this experimental fact, the ground state of the the- 
ory (the vacuum) should not be symmetric under the chiral group. The 
SU{3)l iX" SU{3)r symmetry spontaneously breaks down to SU{3)l+r and, 
according to Goldstone's theorem, an octet of pseudoscalar massless bosons 
appears in the theory (vr, K, rj). 

2.2.1 Effective chiral Lagrangian at lowest order 

The Goldstone nature of the pseudoscalar mesons implies strong con- 
straints on their interactions, which can be most easily analyzed on the 
basis of an effective Lagrangian. Since there is a mass gap separating the 
pseudoscalar octet from the rest of the hadronic spectrum, we can build an 
effective theory containing only the Goldstone modes. Our basic assumption 
is the pattern of Spontaneous Chiral Symmetry Breaking (SxSB) where the 
quark condensate 

{0\uu\0) = {0\dd\0) = (0|ss|0) / (2.13) 

is its natural order parameter. This pattern reads: 

G^ = SU{3)l SU{3)r ^H = SU{3)v (2.14) 

Let us denote (p'^{a = 1,...,8) the coordinates describing the Goldstone 
fields in the coset space G^/H, and choose a coset representative ^((/>) = 
{CLCi4')jCR{4')) £ Cx- The change of the Goldstone coordinates under a 
chiral transformation g = {gLidR) G G^ is given by 

id^) ^ gLiL{(p)h\(t^,g), iR{(p) ^ gRiR{(t^)h\(p,g) (2.15) 

where h{(p,g) £ H is a, compensating transformation which is needed 
to return to the given choice of coset representative ^. Since the same 
transformation h{(p, g) occurs in the left and right sectors, we can get rid of 
it by combining the two chiral relations above into the simpler form 

U{<P) ^ gRU{<p)gl, [/(</.) ^ ^Rimii^^)- (2.16) 

Without lost of generality, we can take a canonical choice of coset rep- 
resentative such that ^R = Cl ^ ^(<A)- The 3x3 unitary matrix 
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U{(l)) = u{<pf = expl —V2^ \ (2.17) 

gives a very convenient parametrization of the Goldstone fields 



75-° + 75^« -^ ^' 



Z^T^v^^-l -" -75^+75^8 ^ I (2-18) 



To get a low-energy effective Lagrangian realization of QCD for the light- 
quark sector (u, d, s), we should write the most general Lagrangian involving 
the matrix U{(p), which is consistent with chiral symmetry. The Lagrangian 
can be organized in terms of increasing powers of momentum (or number of 
derivatives) : 

£e// = 5^^2n (2.19) 

n 

In the low-energy domain we are interested in, the terms with minimum 
number of derivatives will dominate. To lowest order, the effective chiral 
Lagrangian is uniquely given by the term (a proof can be found in [18j ) 



£2 = ^{^,,U^^^'U) (2.20) 

(where {A) = trA) or expanding f7((/>) in power series in the gold field (p 

C2 = lid^cfd^^) + y^((</.9;.</<)(</<9^</.)) + 0(///o'). (2.21) 

With this expansion one obtains the Goldstone kinetic terms plus a tower 
of interactions involving and increasing number of pseudoscalars. 

This effective field theory technique becomes much more powerful if one 
introduces couplings to external classical fields. Let us consider an extended 
QCD Lagrangian, with quark couplings to external Hermitian matrix-valued 
fields Vfj,, ttfj,, s, p: 

^QCD = Cqcd + qi^{vpi + l^a^)q - q{s - i^^p)q (2.22) 

The external fields will allow us to compute the effective realization of 
general Green functions of quark currents in a very straightforward way. 
Moreover, they can be used to incorporate the electromagnetic and semilep- 
tonic weak interactions, and the explicit breaking of chiral symmetry through 
the quark masses: 
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^M 


= Vf,+ 


a-ti 


^M 


= v^- 


Clfi 


S - 


= M + 





eQA^ + ... 

^^^/^ + 7Ji^(^^^+ + ^--) + --- (2.23) 



Here, Q and M denote the quark-charge and quark-mass matrices, re- 
spectively, 



Q = -diag(2,-l,-l), X = diag(m„,md,m^) , (2.24) 

and T+ is a 3 x 3 matrix containing the relevant Cabibbo-Kobayashi- 
Maskawa factors 



'MS 



/ Vud Vu. 
T+ = I . (2.25) 

Y 

The Lagrangian ()2.22p is invariant under the local SU{3)l ® SU{3)ii 
symmetry and in presence of external sources it turns out that the gauge 
fields v^ and a^ can only appear through the covariant derivatives 

D^U = d^U-ir^U + iUl^, D^U^ = d^U^ + iU^r^-il^U\ (2.26) 
and through the field strength tensors 

F^"" = d^r -d^l^-il^Xi F^" = d^'r^ -^''r^'-^[r^',r'']. (2.27) 

At lowest order in momenta, the more general effective Lagrangian con- 
sistent with Lorentz invariance and (local) chiral symmetry is of the form 

m 

C2=^-f{D^U^D^U + U^X + X^U), (2.28) 

where 

X = 2Bo{s + ip) (2.29) 

and Bq is a constant which is not fixed by symmetry requirements alone. 

At leading order in the chiral expansion we are left with only two low- 
energy constants, /o and Sq, which can be related to QCD parameters 
through a matching procedure, which yields 
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,(5£2 



•^° " "';^^°'^'''^^^^^ " ^" ^ ^^^'^' ^^-^"^ 



i<»iti»>-i'-"'- '"'----"- "-"" "' 



Bo = T2(0|^|0) = -^2(01^910), {0\qq\0){2GeV) ~ -[(280±30)MeV]= 



/o can therefore be identified with the pion decay constant, and the 
parameter Bq is proportional to the quark condensate, which takes into 
account the effect of non- vanishing quark masses. This very first term was 
aheady written down by Weinberg [19], and to lowest order it reproduces 
the current algebra results. 

Taking s = M with M = diag{mu,m(i,ms), and p = the x term in 
Eq. (|2.28p gives rise to a quadratic pseudscalar-mass term plus additional 
interaction proportional to the quark masses. Expanding in powers of (j), 
one has (dropping irrelevant constants): 

^2Bo{MiU + t/t)) = i?o| - {M<P') + ^{M^') + o(|J) } (2.31) 

The explicit evaluation of the trace in the quadratic mass term provides 
the relation between the physical meson masses and the quark masses: 



(2.32) 



where 



M^2,= 


= 2mSo, M^o 


= 2rhBo - 


e + 


Oie% 


M'k. 


= (m„ + ms)Bo, 


M^o = 


{rrid 


+ ms)Bo 


K = 


-- Km + 2ms)Bo + e + 0{e^), 


! 






= ^{mu + md), 


^ Bo{ 


ruu - 


-rudf 




'- 4" 


(m. 


-m) ■ 



(2.33) 

Chiral symmetry relates the magnitude of the meson and quark masses 
to the size of the quark condensate. The Eq. (|2.32p imply the old Current 
Algebra mass ratios (Gell-Mann, Oakes and Renner 1958, Weinberg 1977), 

2m {rriu + mg) {md + mg) (2m + 4m,s) 

and (up to 0{m,u — m,d) corrections) the Gell-Mann Okubo (1962) mass 
relation, 

3M2^ = 4M|- - Ml (2.35) 

The lowest-order chiral Lagrangian Eq. (j2.28p encodes in a very compact 
way all the Current Algebra results obtained in the sixties. The nice feature 
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of the chiral approach is its elegant simplicity. Moreover, as we will see 
in the next section, the effective field theory method allows us to estimate 
higher-order corrections in a systematic way. 

2.2.2 ChPT at 0(/) 

Eq. ()2.28p comes from the generating functional Z defined as: 

^iZ[v,a,s,p\ ^ (olTe^J^'^'^^^^'^lO). (2.36) 

This vacuum-to- vacuum transition amplitude generates the Green functions 
of the vector, axial-vector, scalar and pseudoscalar quark currents built out 
of the three flavours n, d and s. Z[v, a, s,p] is associated with the lagrangian 
of Eq. (|2.22|) and admits an expansion in powers of the external momenta 
and of quark masses, Z = Z2 + Z4 + Zq... = J2n ^n- As a consequence of 
chiral symmetry and its spontaneous breakdown, Z2 coincides in the meson 
sector at leading order in ChPT with the classical action 

Z2= f d^xC2iU,v,a,s,p). (2.37) 

where £2 is the same Eq. (j2.28p . 

At O(p^), the generating functional consists of a contribution to ac- 
count for the chiral anomaly, the one-loop functional originating from the 
lagrangian (|2.28p and an explicit local action of order p^. The most gen- 
eral Lagrangian, invariant under parity, charge conjugation, the local chiral 
transformations, and able to generate a local action of order p"^ is given by 

m 

6^\U,DU) = Li{D^U^D''Uf + L2{D^U^D^U){D''U^D''U) + 

+ L■i{D^U^D^'UD^U^D^U) + U{D ^U^ D^'U) {U^ x + X^U) + 

+ HiD^U^D^'UiU^x + X^U)) + L^{U^x + X^Uf + 

+ Lj{Uh - xW + Lsix^Ux^U + U^xUh) - 

+ Hi{Fti^,F^'' + FLf,,F^'')+H2{x^x) (2-38) 

where F^^^ = div" ± a") - d{v^' ± a^') - i[v^' ±a^',v'' ± a""]. 

The terms proportional to Hi and H2 do not contain the pseudoscalar 
fields and are therefore not directly measurable. Thus, at O(p^) we need ten 
additional coupling constants Li to determine the low-energy behavior of 
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the Green functions, the low-energy constants (LECs). These constants pa- 
rameterize our ignorance about the details of the underlying QCD dynamics 
and are needed in order to make reliable phenomenological predictions. As 
with any other effective field theory, these LECs play the role of coupling 
constants and contain the information coming from the integration of the 
heavy degrees of freedom not explicitly present in the Chiral Lagrangian. 
We will exemplify this process in subsection 12.5.11 with meson resonances. 

However, the increase in the number of operators as one goes to higher 
orders in the chiral expansion, together with the poorly-known associated 
low-energy couplings, makes computations difficult already at 0{p'^). In 
principle, all the chiral couplings are calculable functions of Aqqd and the 
heavy-quark masses. At the present time, however, our main source of 
information about these couplings is low-energy phenomenology. At 0{p'^) 
there are 2 LECs, at 0{p'^) 10 [l7j, and at 0{p^) the number of constants 
becomes more than a hundred [20^ [21] [22] . In the electroweak sector this 
proliferation of constants appears already at 0{p^) [23l[Mj- Although in 
principle these low-energy constants may be computed on the lattice, in 
practice this has only recently [2S] been accomplished in a few cases for the 
strong Chiral Lagrangian at 0{p'^). 

Since ChPT is a quantum field theory, the next step in our description is 
take into account quantum loops with Goldstone-boson propagators in the 
internal lines. The chiral loops generate non-polynomial contributions with 
logarithms and threshold factors, as required by unitarity. 

ChPt is an Effective field theory, then have to be provided with a power 
counting rule. Actually, in our case, at 0{p'^), the diagrams that contribute 
are dictated by the relation ([H]): 

d = 2 + Y,Nn{n-2) + 2NL, n = 2,4,6,... (2.39) 

n 

where Nn is the number of vertices coming from 0{p^) operators, and 
Nl is the number of loops. 

At one loop (in C2), the ChPT divergences are O(p^) and are therefore 
renormalized by the low-energy couplings in equation (j2.38p : 



Li = L\{iJL) + Ti\, Hi = Hl{iJL) + Ti\, (2.40) 

where 

^=£{d^-5N(^')+m + i]}^ (2^41) 
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and the evolution under the renormaUzation group is then given by 



L[M = L[(Mi) + ^log(^) (2.42) 



lOTT^ /X2 

where Fj take the values due to the explicit calculation of the one-loop 
generating functional Z4 [T7] gives: 



^1-32' 


r2 = ^, r3 = o, 


^-^ 




P 3 


r« = i44<^-«' 


r7 = o, 


^8-48' ^^-4' 


rio = - 


1 

"4' 


f-i 


f-l- 


Decay constants at 0{p'^) 











For illustrative purpose we also show the O(p^) calculation of the meson- 
decay constants in the isospin limit {rriu = ra^ = rh) [ITJ: 



fK= / {1 -!/..- i/... - !/.,s + '-p'Llij.) + '-^^^Lli^^)} , 
U= / |l - 3A.X + '-^Llif,) + M4+l^LK;u)| , (2.44) 



M'p ,^JMI 



where 

masses could be found in [5] and Li in Table 12. 3[ 

Electromagnetic Form Factors 

As a second example, we present the main results for the electromagnetic 
form factors (that will we used in chapter 6) . That is a particular case where 
higher-order local terms in the chiral expansion are important. Ignoring 
those will lead to a wrong result compared to the experimental value. 

At Oip^) the electromagnetic coupling of the Goldstone bosons is just 
the minimal one, obtained through the covariant derivative. The next-order 
corrections generate a momentum-dependent form factor: 



^(r^)t%^ + -; i^f(?^) = ^ <'"-"" 



Ffiq') = 1 + '-{rX<i' + •••; 4i<f) = -JrX^' + ■ (2-45) 



The meson electromagnetic radius {r^)'y gets local contributions from 
the Lq term, plus logarithmic loop corrections (|26j): 
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(^^)r^i^-.^f2>oJ^V.o.r#)), (2.40) 



'Mk 



(-^)f =-T7^1og 



2.3 QCD Sum Rules, the Operator Product Ex- 
pansion 

As a consequence of asymptotic freedom the theoretical results obtained 
from QCD can be compared easily with the experimental situation for the 
hard processes, which is, at short distances the effective coupling constant 
as becomes small and the interaction can be treated perturbatively. On 
the other hand, as we have said, any comprehensive theory of the strong 
interaction must include long distance dynamics as well. Actually, quark 
interaction within hadrons is strong, since it binds quarks into insepara- 
ble groups. However, nowadays we do not have any accurate quantitative 
framework within QCD for dealing with the strong interaction regime and 
the evaluation of the hadron spectrum. 

An interesting approach was started in 1979 by Shifman, Vainshtein and 
Zakharov, assuming that confinement exists (|27j). The effects of confine- 
ment can be described through the use of a few parameters, the so called 
condensates, and allows to obtain many hadronic properties through an ap- 
propriate use of sum rules. One of the main ingredients of this approach is 
the Operator Product Expansion (OPE), proposed by Wilson in 1969 (|^H]). 
who hypothesized that the singular part, as x — )• y, of the product A{x)B{y) 
of two operators is given by a sum over other local operators 

Aix)Biy) ^ Yl ^c^'i^ - y)C{y) (2.47) 

c 

where Fq {x — y) are singular c-number functions. This OPE exists 
for the free scalar and spinor field theories and for renormalized interacting 
fields to all orders in perturbation theory. In every case, they are valid 
for any elementary or composite local fields: A and B can be elementary 
scalar or spinor fields or local currents or the stress-energy tensor or any 
local Wick product in a free field theory. Dimensional analysis suggests that 
F^^{x—y) behaves for x — >• y like power dc—dA—ds oix—y, where do is the 
dimensionality of the operator O in powers of mass or momentum. Since do 
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increases as we add more fields or derivatives for an operator O, the strength 
of the singularity of Fq {x — y) decreases for operators C of increasing 
complexity. The remarkable thing about the operator product expansion 
is that it is an operator relation. It is the decrease of the singularity in 
Eq. (j2.47p with operators C{y) of increasing complexity that makes this 
expansion useful in drawing conclusions about the behavior of the product 
A{x)B[y) for x ^ y. Renormalization effects modify the power counting: for 
asymptotic free theories F^^ [x — y) behaves like the power dc — dA — dB — j 
oi X — y suggested by dimensional analysis where 7 is a constant. 

The corresponding statement in momentum space is that for A; — )• 00, 

/ d'^xe-''''-'A{x)B{0) -^ Yl yc^{k)C{0) (2.48) 

•^ c 

and correspondingly 

/ d^xe-''^'-'T{A{x)B{0)} -^ ^ U^^{k)C{0) (2.49) 

"^ c 

where V(j (k) and U^ {k) are functions of k^ that for large k decrease 
increasingly rapidly for more and more complicated terms in the series. 

Let us show an example, following [29], to illustrate this process which 
will be useful to perform our own calculation in section 14.2.11 Ref. |29) 
suggests to calculate the p meson mass using the QCD sum rules. The first 
point begins with the isovector part of the electromagnetic current: 

JA* = l(M(x)7^n(x) - d{x)-i^'d{x)) (2.50) 

These composite operators are the ones which appear in the OPE Eq. ()2.47p . 
Let us now introduce the two point function: 

H^'^iq) = ijd'xe'''-^n I T {jl^{x)Jl^{Q)) \ Q) , (2.51) 

where | 0) is the physical vacuum of the theory. Since Eq. ()2.50p is a 
conserved current we can write 



nr(5) = (^'^9'^-'zV'^)np(g') (2.52) 



with 



np(5') = - p-i)g2 /^''^^^^"(^ I T {ji;,)ix)Jip),{0)) I O) , (2.53) 
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Dimension 


Operators 





1 


4 


mAqaix)qaix) 

Fi:'{x)F-,{x) 


6 


qa{x)Tqo,{x)qp{x)Tqp{x) 
qc,{x)T{\'')^pqp{x)q^{x)T{X')^sqs{x) 

fabcFr{x)Ft,{x)F',P{x) 



Table 2.1: Lowest dimension scalar operators. 



where D is the number of space-time dimensions. Notice that the di- 
mensions of the operators appearing here are d{Jp) = M and dilip) = 
M . Now, we are in the stage to use the Operator Product Expansion, 
then we can write: 



hm Iip{q^) = - .^ \. ^ V(0 I O„(0) I n) fd^'xe^'^-^Cnix) (2.54) 

where On(0) are the local operators that appears in Eq. (|2.47p . For 
large values of q^ the behavior of the integral appearing in the r.h.s. of this 
equation is: 



(q 



2^[d{On)+2-D]/2 



(2.55) 



where it will be necessary to consider in (|2.54p only operators Onix) such 



that 



d(On) <{D-A) + 2N 



(2.56) 



where A^ is used to identify the terms in Ilp{q'^) that decreases faster 
than (g^)~^ in the limit q'^ — )• oo. 

We will be interested in the lowest dimension scalar operators and these 
are shown in table 12.11 

The operators On are conveniently classified according to their Lorentz 
spin and dimension d. We will consider only spin-zero operators since only 
these contribute to the vacuum expectation value. Within the standard per- 
turbation theory only the unit operator would survive in Eq. ()2.54p . but the 
non-perturbative effects induce non- vanishing vacuum expectation values for 
other operators as well and they are the so called condensates. Therefore, 
the non-perturbative effects of QCD introduce power corrections of type 
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l/(g^) , A^ > 1, to the perturbative calculation. Since QCD is asymptot- 
ically free, the calculation of the coefficients C„ is reliable. The expansion 
coefficients in Eq. ()2.54p are calculated as a series in Ug expansion. The 
final results concern a perturbative and a non-perturbative contributions. 
The lowest order perturbative contribution is: 

n(p)(?^) = -^(i + v)i"(-^') (2.57) 

The non-perturbative contributions contains the quark condensate and 
the gluon condensate: 

11 1 

^(p)^"!^) = 2rvF ^"'" <uu>+md<dd>) + ^^.2\2 (^^a"'{x)F;i^{x)) 

(2.58) 
always up to terms of order a^. From the second order perturbative 
result due to the contribution of the condensate of four quark fields, we 
have (using isospin symmetry, < uu >=< dd >): 

112 1 

n(p) (Q^) = ^™^ < un >2 — (2.59) 

Taken into account all these contributions, i.e, Eqs. (j2.57p . (|2.58|) and 
(|2.59p . we get from QCD the following information (using —q^ = Q^): 



Hp)l^^ = -i(l + v)^^(^') + ^^' 



n(p)(Q^) = -TT^ 1 + — ln(Q2) + -T^imu <uu> +md < dd >X2.60) 



■^(^F™C(-))-^-.<.n>^i, 



where we have neglected constants, higher order terms in the OPE, 

2 

radiative corrections to the condensates and ^ terms. Except for the first 
term in Eq. (|2.60p . this is an expansion in terms of Q~^ when Q^ — )• oo. 

This way to proceed here illustrated will be useful in chapter H] and in 
further calculations of the OPE for two-point Green's functions. 

2.4 The Large- iVc limit 

QCD, as a non-Abelian gauge theory based on the gauge group SU{3) can 
be understood 2Q1 EI] by studying a gauge theory based on the gauge group 
SU{Nc) in the hmit iVc — )• oo Q One might think that letting A'^c — ^ oo 
would make the analysis more complicated; or SU{Nc) gauge theory could 



For a pedagogical review see for example [32l [33] 
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not be related to QCD because the difference between Nc ^ oo and Nc = 
3. We will soon see that SU{Nc) gauge theory simplifies in the A'^c — t- oo 
limit, that the true expansion parameter is 1/Nc, not Nc- In fact, the 1/Nc 
expansion is equivalent to a semiclassical expansion for an effective theory 
of color singlet mesons and also for baryons|34]. Results for QCD will we 
obtained from Nc ^- oo limit by using 1/Nc = 1/3, with good agreement 
with experiment in both meson and baryon sectors. A way to understand 
the validity of this approximation is done by 't Hooft and later Witten 
|30t l3T| and are decomposed in two ways. The first one related to why 
the perturbation theory is successful in QED. In this theory, the typical 

2 

expansion parameter is |;^ where e = 0.3 is the electric charge. If the typical 
expansion parameter had turned to be 47re^, perturbation theory would not 
have been very successful for e as large as 0.3. The other way asks us how 
small X must be for a series ^ a„x" to be dominated by the first few terms. 
The answer depends entirely on how large are the coefficients a„. If the 
coefficients are very small, x = 1/3 can be considered a small number. 

Therefore, the large Nc expansion [301 ET] stands out as a very promising 
analytic approach capable of dealing with the complexities of nonperturba- 
tive QCD while, at the same time, offering a relatively simple and manage- 
able description of the physics. This description includes some aspects that 
are not understood in the QCD context such as: 

1 . The suppression in hadronic physics of the qq sea; the fact that mesons 
are approximately pure qq states; the absence, or at least suppression, 
of qqqq exotic states. 

2. The Zweig's rule; the fact that mesons come in nonets of flavor SU(3); 
the decoupling of glue states. 

3. The fact that multiparticle decays of unstable mesons are dominated 
by resonant two body final states, when these are available. 

4. The Regge phenomenology and the success of a phenomenology that 
describes the strong interactions in terms of tree diagrams with ex- 
change of physical hadrons. 

For instance, mesons are qq states with no width, the OZI rule is exact 
and there is even a proof of spontaneous chiral symmetry breaking [35] (more 
details will be discussed in sec l2.4.1]l . Furthermore, interest in studying 
SU(N) gauge theories in the large Nc limit has increased due to the discovery 
of a duality of some highly supersymmetric gauge theories to gravity [36] . 
although the real relevance of this connection for QCD still remains to be 
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seen. However, in spite of all this, the fact that no solution to large- A'^c QCD 
has been found still poses a serious limitation to doing phenomenology. For 
instance, in order to reproduce the parton model logarithm which is present 
in QCD Green's functions in perturbation theory, an infinity of resonances 
is necessary whose masses and decay constants are in principle unknown. 

2.4.1 QCD in the Large-iVc limit 

When describing QCD, we will see that the coupling constant has been 
chosen to be g/\/Nc, rather than g, because this will lead to a theory with 
a sensible (and non-trivial) large N^ limit (further details will be presented 
in the next subsection). The field strength is then: 

-t^ fiij ^ t)^Ay — Oi/A^ + z— ^=- [A^, Ai^\ , 
and the Lagrangian is 

C = --TiF^^F'"' + J2^k {il^D^ - nif) qu- (2.61) 

/=i 

The large Nc limit will be taken with the number of flavors Np fixed. It 
is also possible to consider other limits, such as Nc ^ oo with Np/Nc held 
fixed [37]. 

The Lagrangian of the theory stands as in Nc = 3 QCD (after all, it is still 
a Yang-Mills Lagrangian) except for the fact that we have changed the gauge 
group. That means that each quark field, since they sit in the fundamental 
representation, appears as an A^c-plet. Gluons, on the contrary, live in the 
adjoint representation and enlarge their number to iV^ — 1. In practice, 
it is common to approximate this to Nc', in other words, we are skipping 
the tracelessness constraint (we are taking U{Nc) instead of SU{Nc)). The 
missing gluon is numerically unimportant at sufficiently large Nc- Besides, 
it can be shown that the abelian factor is indeed suppressed at large Nc- 
Notice that the theory we are looking at differs from Nc{= 3) QCD in that 
there exist far more gluons than quarks (the former scale with A*"^ while the 
latter only with Nc). 

Our aim will be to show the topological ordering of diagrams induced 
by the large-A'c power counting scheme. For clarity, it is convenient to use 
't Hooft double-line notation. 
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Figure 2.1: Double line notation for quarks, anti-quarks and gluons. 



muuw 





Figure 2.2: Feynman diagrams with its counterparts in the double line notation 
introduced by 't Hooft. 



The Double-Line Notation and Planarity 

In perturbative QED there is only one coupling constant which shows 
up to couple ferniions and antifermions. That is why Feynman diagrams 
are so useful to organize calculations in powers of the coupling constant: 
you only need to count the number of vertices. In QCD, however, the 
ordinary coupling constant gg is not really a free parameter, because in view 
of the renormalization group, it is absorbed into defining the scale of masses. 
Therefore, we have to change our strategy. In Yang-Mills theory at iVc —)■ oo, 
the mesons and glue states are free, stable, and non-interacting. Meson 
decay amplitudes are of order l/\/77^, and meson- meson elastic scattering 
amplitudes are of order l/Nc- These elastic amplitudes are given, as in 
Regge phenomenology, by a sum of tree diagrams involving the exchange, 
not of quarks and gluons, but of physical mesons. 

In the 1/Nc expansion, we need to keep track not only of the vertices 
of the theory (we will show later on that the coupling constant at large- A'c 
is color-dependent) but also of color flow inside the diagram. We would 
like to have a pictorial approach to be able to determine in an easy way 
the scaling of physical amplitudes with the color factor. The double-line 
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notation introduced by 't Hooft makes it transparent to extract the color 
scaling of a given amplitude altering only slightly the Feynman diagrams we 
are used to. Consider the following pictorial recipe: represent every quark 
line by a color line, right-faced arrow meaning color flow, left-faced arrow 
meaning anticolor flow, as is shown in figure [2TT1 

This does not change the Feynman picture much. When it comes to 
gluons, however, we have to interpret their color indices as a system of color 
and anticolor line. Figure 12.21 shows some examples of converting Feynman 
diagrams to double-line diagrams. 

From the gluon loop contribution to the gluon propagator, figure 12.21 
it is easy to see that even after the color quantum numbers of the initial 
and final states are specified, there are still N possibilities for the quantum 
number of the intermediate state gluons. As a result, this diagram receives 
a combinatoric factor of N^. On the other hand, in the same figure 12.21 
there is also a factor of coupling at each of the two interaction vertices. If 
we want the one-loop gluon vacuum polarization to have a smooth limit 
for large- A'c, we must choose the coupling constant to be g/\/T^, where 
g is to be held fixed as Nc becomes large. With this choice, introduced 
by Witten \30\ I31j . the two vertex factors of gj^fWc in figure [22] combine 
with the combinatoric factor of N^ to give a smooth large-A^c behavior: 
{g/y/N^c)"^ X Nc = g'^, independent of N^. Thus, the vanishing for large- A^c 
of the coupling constants cancels the divergence of the combinatoric factor, 
to produce a smooth large- iVc limit (as expected). 

Due to the normalization for the coupling constant we choose, in order 
to survive as Nc ^ oo, a. Feynman diagram must have combinatoric factors 
large enough to compensate for the vertex factors. Only a certain class of 
Feynman diagrams, the so-called planar diagrams, have the combinatoric 
factors large enough to just cancel the vertex factors. The other diagrams 
vanish for large- A'^c. The large- A'c limit is therefore given by the sum of the 
planar diagrams. 

As we have mentioned above, figure 12.21 has a combinatoric factor of 
Nc- In the right side of this figure the gluon loop is been redrawn in the 
double line notation. Following the arrows, the color lines at the edge of the 
diagram are contracted with those of the initial and final states. However, 
at the center of the figure, there is a closed color line that is contracted 
only with itself. The color running around this loop is unspecified even 
when the initial and final states are given, and the sum over the quantum 
number of this loop gives a factor of Nc- This is the combinatoric factor of 
Nc associated with the upper left diagram of the figure 12.21 

Using the double line notation, it is not difficult to determine whether a 
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Figure 2.3: Two-loop contribution to the gluon propagator, in the ordinary notation 
and in the double line notation, respectively. 



given diagram survives in the large- A'c limit. For example, in figure [231 there 
is drawn, both in the ordinary notation and in the double line notation, 
a two-loop contribution to the gluon propagator. This diagram has four 
interaction vertices, each contributing a factor of l/-v/iVc, but it has two 
closed color loops that are self-contracted, each contributing a factor of Nc- 
Altogether, the diagram is of order {l/^/N'c)^ x A'^^ = 1 and so survives in 
the large- iVc limit. 

Therefore, exist a relation between the vertices and the loops to indicate 
when a diagram survives or not. A counterexample of a diagram that do not 
survives is drawn in figure 12.41 both in ordinary and double line notation. 
This diagram has six interacting vertices, but only one large and tangled 
closed loop. The diagram is, therefore, of order (l/^/iVc)^ x Nc = 1/A^ and 
vanishes like 1/N^ as Nc becomes large. The basic difference between this 
diagram and the previous one is this diagram is nonplanar, it is impossible 
to draw this diagram on the plane without line crossing. The diagrams of 
figures [2^2] and [231 are, by contrast, planar, they can be drawn on the plane. 

With this example, we arrive to the first "selection rule" in the large- Ac 
limit, that is nonplanar diagrams are suppressed. 

There is also a second " selection rule" that reflects the fact that for large- 
Nc there are A^ gluon states but only Ac quark states, so that diagrams 
with internal quark lines have fewer possible intermediate states and smaller 
combinatoric factors. 

Consider the one quark loop contribution to the gluon propagator. This 
diagram is drawn in figure 12.51 both in the ordinary and the double line 
notation. Because the quark propagator corresponds to a single color line, 
not two, the closed color line present in figure [2^2] is absent in figure [231 As 
a result, figure [231 has no large combinatoric factor, and its only dependence 
on Nc comes from factors of 1/ ^/Wl at each of the two vertices. So figure 
12.51 vanishes like 1/Ac for large-A'c- 

In sum, there are two selection rules for Feynman diagrams in the large- 
A^c limit; 



Nonplanar diagrams are suppressed by factors of 1/A, 
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Figure 2.4-' A nonplanar example diagram as a contribution to the gluon propagator 
in the ordinary and double line notation 




Figure 2.5: One quark loop contribution to the gluon propagator, in the ordinary 
and in the double line notations. 



• Internal quark loops are suppressed by factors of l/N^. 

The leading diagrams for large- A'^c are the planar diagrams with a mini- 
mum number of quark loops. 

Properties of mesons in the Large- A'^c limit 

Before describing mesons for large- A'^c, it is necessary to make an impor- 
tant assumption. QCD is a confining theory at Nc = 3, and we will assume 
that the confinement persists also at large- A'^c. Combining the assumption 
of color confinement and the knowledge from large- A'c and planar diagrams, 
the large- A'^c theory has the following properties: 

• Mesons for large-A'c are free, stable, and noninteracting. Mesons 
masses have smooth limits for large-A'c, and the number of meson 
states is infinite. 

• Meson decay amplitudes are of order l/^/Wl] meson- meson elastic scat- 
tering amplitudes are of order 1/Nc, and are given by a sum of tree 
diagrams involving the exchange of physical mesons. 
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Figure 2. 6: A representative Feynman diagram contributing to the two-point func- 
tion together with its counterpart in double line notation. The cut illustrates the 
fact that in planar diagrams only one-particle intermediate states are allowed. 



• Zweig's rule is exact at large- A'c; singlet-octet mixing and mixing of 
mesons with glue states are suppressed, so that mesons come in nonets; 
and mesons for large- A'^c are pure qq states. 

The first point to establish is that the operator J(x), as a generic quark 
bilinear, acting on the vacuum, creates, in the large- A'^c limit, only one-meson 
states. It is equivalent to claim that the only singularities of the two-point 
function of J are one-meson poles. In other words, to lowest order in 1/Nc, 

{Ak)J{-k))=Y^^^ (2.62) 

n ^ 

Here, rrin is the mass of the n meson, and a„ = (0| J|n). To show this 
property, one can cut the leading contribution to the two-point function 
of J{x), an see that the only intermediate states that appears are one- 
meson intermediate states, as it is shown upper right of figure 12. 6t where the 
intermediate states (the two gluons) go with one quark and one antiquark. 
In a confining theory, the qq pair are always bound together into a meson. 

Following the previous properties, we can write the two-point functions 
with an one-meson intermediate state as it is shown in figure 12. 7i Notice 
that in the case of three-point correlation functions, the diagrams, figure 
\'2.S\ may be of two types. To determine the amplitude for two-body decays 
such as A ^ BC, we will began with the fact that three point function is 
of order Nc, because in free field theory it is given by the one-loop diagram 
in figure [231 and we know that the more elaborate planar diagram have the 
same dependence on Nc that the free field theory has. But this amplitude 
has a term which is of the form {0\J\m)^Tmmm- Since {0\J\m) is of order 
\/^! '^mmm must be of order l/^/Wl. By similar arguments, we can see 

-(k—2) 

that a local vertex with k mesons is of order \/Nc 

It is also possible to extend this analysis to include glue states. Ampli- 
tudes with arbitrary number of mesons and glue states are given, to lowest 
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Figure 2. 7: Two-point correlation function in the Large-Nc limit 
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Figure 2. 8: Three-point correlation function in the Large-Nc limit represented as a 
summation of two possible ways to interact. 



order in 1/Nc, by sums of tree diagrams. In these diagrams, the general 

local vertex with k mesons and / glue states is of order Nc ^ 

In conclusion, with the two selection rules, it is easy to understand the 
meson phenomenology explained at the beginning of this section in terms of 
the Large- A^c limit. 



2.4.2 Chiral Perturbation Theory in the Large-A^c hmit 

Once we have seen how to understand QCD in the Large- A'^c limit, one 
can proceed to study its low-energy representation. Actually, the U{Nf)l x 
U{Np)ji chiral symmetry of Large- A'^c QCD is spontaneously broken to a di- 
agonal U{Np)v vector symmetry, resulting in (pseudo-) Goldstone bosonrl^' 



""The axial anomaly is l/Nc- See section \%5 




Nr. 



Figure 2.9: Three-point correlation function loop in the Large-Nc limit as a free 
field theory. 
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This form for the breaking can be proven in the large- A'^c hmit |35) . The 
low-energy interactions of the pseudo-Goldstone bosons of QCD, the vr, K, 
rj, and now including rj' (we will se why soon), can be described, as in section 
12.2. It in terms of the effective Lagrangian, the Large- A'c Chiral Lagrangian 
|381l39j. which can be computed by evaluating the QCD functional integral 
with sources for the pseudo-Goldstone bosons (the source terms are actually 
fermion bilinears). In the large- A'c limit we have seen that the leading order 
diagrams that contribute to correlation functions of fermion bilinears are of 
order Nc, and contain a single quark loop. This implies that the leading 
order terms in the chiral Lagrangian are of order Nc- These leading order 
terms can be written as a single flavor trace, since the outgoing quark flavor 
at one vertex is the incoming flavor at the next vertex. Similarly, diagrams 
with two quark loops have two flavor traces, and are of order unity, and in 
general, those with r quark loops have r traces, and are of order N^~^- This 
simple rule will help us on the following lagrangian description. 

The chiral Lagrangian in the large- A'c limit is written in terms of a 
unitary matrix 

[/ = g2ffl//o^ (2.63) 

where /o ~ 92 MeV is the pion decay constant in the chiral limit, and 



^S+^+^ -" ^" ^ 



V2 



° V _L V' 1^0 



K- K^ -^ + ^ 



(2.64) 



is the matrix of pseudo-Goldstone bosons. The rj' has been included since it 
is related to the octet (now nonet) pseudo-Goldstone bosons in the large- A'c 
limit, by Zweig's rule (from SU{3) to U{3), the dual group es A'^^ = 9). The 
order p"^ terms in the chiral Lagrangian (in presence of external fields as in 
sect. 2.2.1) are 

£(2) = i^i^Df^uD^U-^ + BoiM^U + MU-^)), (2.65) 

where Ai is the quark mass matrix in the QCD Lagrangian. The first 
term is order Nc, since /o oc vNc (as we have said before). The second term 
in Eq. ()2.65p also has a single trace and is of order Nc, so Bq is of order unity. 
The U field can be expanded in powers of H/ fo, thus each additional meson 
field has a factor of l/fo oc 1/vNc, which gives the required 1/vNc sup- 
pression for mesons derived earlier. The effective Lagrangian Eq. (j2.65p has 
then an overall factor of Nc- Graphs computed using the chiral Lagrangian 
have a 1/Nc suppression for each loop. 
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Li(Mp) Value Order 

2Li - L2 

La 

Lg 

Lj 

L2 

L3 

L5 

Ls 

L9 

Lio 

Table 2.2: Experimental values for the coefficients of the order p'^ terms in 
the chiral Lagrangian at Large- A'c- Values are from ref. [4H I42j 



0.6 ±0.5 


1 


0.3 ±0.5 


1 


0.2 ±0.3 


1 


0.4 ±0.15 


1 


1.3 ±0.7 


Nc 


4.4 ±2.5 


Nc 


1.4 ±0.5 


Nc 


0.9 ±0.3 


Nc 


6.9 ±0.2 


Nc 


5.2 ±0.3 


Nc 



The order 0{p^) terms in the chiral Lagrangian were presented in Eq. ()2.38p |17j . 
It turns out that one finds the A^c-counting rules for the LECs: 



0{Nc) : Li,L2,L3,L^,Ls,Lg,LiQ 

0(1): 2Li-L2,L4,L6,^7 

These are the A'^c-counting rules given in ref. pj'J, with the exception of 
L7, which is taken from ref. [4Dj. In ref. ^7\, L7 was argued to be of order 
N^. The experimental values for the Lj's in the Large- A'^c limit are given 
in Table 12.21 The terms of order Nc are systematically larger than those of 
order unity. 

Higher derivative terms in the chiral Lagrangian are suppressed by pow- 
ers of the chiral symmetry breaking scale A-^ ~ 1 GeV J19l I43j . In the 
large-A'^c limit A^ is of order unity and so stays at around 1 GeV. Loop 
graphs in the chiral Lagrangian are proportional to l/(47r/^)^ and are of 
order l/Nc- Thus in the large- A^c limit, the chiral Lagrangian can be used 
at tree level, and loop effects are suppressed by powers of 1/Nc- 

2.5 Resonances in the Large- A^c limit 

When we discussed chiral symmetry and its implementation to study the 
low energy dynamics of the pion octet fields, we estimated the radius of 
convergence of the chiral expansion to be A^ ~ IGeV. This scale lies close 
to the first resonance multiplet and then a natural extension of the theory 
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would exploit chiral symmetry by incorporating these higher energy exci- 
tations states in a Lagrangian. In such a way, this Lagrangian would then 
allow a prediction of the LECs of the strong interactions in terms of the 
masses and decay constants of the included resonances. 

The building of an effective Lagrangian out of Goldstone bosons and 
resonance excitation fields in a chiral-invariant way can be achieved once 
the resonance fields are embodied with a chiral representation. This is the 
avenue followed by Resonance Chiral Theory (RChT) which is a description 
of the Goldstone-resonance interactions in a chiral invariant framework \42\ 
HH I45j . Alternatively to the chiral counting, it uses the 1/Nc expansion of 
QCD in the limit of a large number of colors as a guideline to organize the 
perturbative expansion. At leading order just tree-level diagrams contribute 
while loop diagrams yield higher order effects. Integrating out the heavy 
resonance states leaves at low energies the corresponding chiral invariant 
effective theory, ChPT. Many works have investigated various aspects of 
RChT: Green functions (061 |171[l8lll5lll9l[50l|5ll[52l[53l[5l[55l [56]), 
apphcations to phenomenology ([MlEZlESlESEniEIlEaES])- 

In RChT the pseudo-Goldstones enter through the exponential realiza- 
tion U = exp(i0/v2/). As the standard effective field theory momentum 
expansion is not valid in the presence of heavy resonance states, RChT takes 
the formal large- A^c limit expansion as a guiding principle |30t [31] . The la- 
grangian can be then organized according to the number of resonance fields 
in the interaction terms, 

^RChT = ^GB + ^Ri + ^RiRj + ^RiRjR^ + ••■ ■, (2.66) 

where Ri stands for the resonance multiplets and the first term in the r.h.s. 
of the equation contains the operators without resonance fields, Eq. (|2.20p . 
Let us consider a chiral-invariant Lagrangian £/j, describing the cou- 
plings of resonance nonet multiplets V-^'^ {1 ), ^^''(1"'""'"), 5,(0"'""'") and 
Pi(0~"^) to the Goldstone bosons ^[M], and the kinetic terms [HI [12], 

^^ = E {^ (V^r/+..> + ^ {Vr^.n^) + 0= (<V-..)(2.67) 
+ Cd^ (SiUf'Uf,) + Cm, {SiX+) + idm, {PiX~) 
where 
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with F^'^^ the field-strength tensors of the l^ and r'^ flavor fields defined 
in Eq. ^(TWif and 

X± = u^XU^ ± ux^u- 

The resonance couplings Fy^, Gy,, Fa,, Cd^, Crm and drm are of O {^/N^). 

The lightest resonances have an important impact on the low-energy 
dynamics of the pseudoscalar bosons. Below the resonance mass scale, the 
singularity associated with the pole of a resonance propagator is replaced by 
the corresponding momentum expansion. Therefore, the exchange of virtual 
resonances generates derivative Goldstone couplings proportional to powers 
of 1/M^. At lowest order in derivatives, this gives the large-A^c predictions 
for the 0{p^) couplings of chiral perturbation theory |44^ I64j . 

All these couplings are of 0{Nc), in agreement with the counting indi- 
cated in Table [221 while for the couplings of 0(1) we get 2Li — L2 = L4 = 
Lg = L7 = (as we will see experimentally, except for Lj, in table [23]) . 

Owing to the U{1)a anomaly, the r/i field is massive and it is often 
integrated out from the low-energy chiral theory. In that case, the SU{3)l (8) 
SU{3)ji chiral coupling L7 gets a contribution from r/i exchange (TTHMlEl] 
is: 

2.5.1 Integrating out Resonances 

Once the Lagrangian C-ji is chosen, the natural step is to integrate out the 
resonances. We are left with the original chiral lagrangian, the coefficients 
being functions of hadronic parameters (masses and couplings), to be related 
to the Lj's through a matching procedure. The results are as follows, now 
with couplings of 0(1) and 0{Nc) together [il]ll2]: 
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r<2 „2 :2 
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^2=l4+ + + + + 
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VI 



S ^^"P 



(2.69) 



Upon comparison with experiment one gets the values shown in Table 



This Table 12.31 shows a remarkable agreement, strongly supporting res- 
onance saturation. Still, there are two flaws in the argument: first the 
value /i* = Mp is, though very suggestive, completely arbitrary. The right 
way to proceed would be to determine /x* from Cn at one loop level. The 
renormalized couplings LKij,) depend on the arbitrary scale of dimensional 
regularization fi. This scale dependence is of course canceled by that of the 
loop amplitude, in any measurable quantity. However, contrary to Cqcd in 
Eq. (j2.22p . Cr has no obvious expansion parameter. One way out is rely on 
the fact that Cr is meant to be an EFT of QCD and then we can resort to 
Large- A'c methods to tackle the problem. 

The second flaw is related to the values used for the masses My, Ma-, Ms, 
Ms, and Mp to obtain the estimations shown in Table [2131 To describe the 
Lagrangian Eq. ()2.66p . one uses instead the Lagrangian Cr, Eq. ()2.67p . where 
one has done the assumption that at low energies the main contribution 
coming from the resonant sector is mainly dominated by the lowest lying 
states. However, we already know that in the large-A'^c limit of QCD, the 
Green's functions are described by a infinite number of states. The the 
precedent Lagrangian has also to be described with a infinite number of 
states. In practice, however, one relies on Eq. (j2.67p as an approximate 



Numerical values for masses and couplings are given in Refs. [41[|42] , 
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Li{M,) 


experimental values (xlO ^) 


V 


A 


S 


Si 


Vi 


total 


Li 


0.7 ±0.3 


0.6 





-0.2 


0.2 





0.6 


L2 


1.3 ±0.7 


1.2 














1.2 


L3 


-4.4 ±2.5 


-3.6 





0.6 








-3.0 


La 


-0.3 ±0.5 








-0.5 


0.5 





0.0 


H 


1.4 ±0.5 








1.4 








1.4 


U 


-0.2 ±0.3 








-0.3 


0.3 





0.0 


Lj 


-0.4 ±0.15 














-0.3 


-0.3 


Ls 


0.9 ±0.3 








0.9 








0.9 


L9 


6.9 ±0.2 


6.9 














6.9 


Lio 


-5.2 ±0.3 


-10.0 


4.0 











-6.0 



Table 2.3: Values of low-energy constants appearing in Eq. (|2.38|) 



Lagrangian which only includes one resonance per channel. We will see in 
the next chapter than in that situation, i.e., using a finite number of states 
instead of an infinite one leads to a rational approximation which poles and 
residues that may have nothing to do with the physical ones. Then, using 
the values for the lowest meson states from the Particle Data Book to obtain 
the estimations shown in Table [231 may induce a systematical source of error 
which is not been properly considered up to now. 



2.5.2 Minimal Hadronic Approximation 

We have seen in sec. 12.41 that the operator J(x), as a generic quark bilinear 
creates, in the Large- A'^c limit, only one narrow meson state. In other words, 
largeA'^c QCD is able to constrain the analytical form of QCD Green func- 
tions to be meromorphic functions, Eq. ()2.62p . In order to reproduce the 
parton model logarithm, the sum of states in Eq. (|2.62p extends to infinity. 
Analytically, the infinite number of poles and residues of the Green function 
could be known if the solution to large- A'^c QCD is known. However, as we 
have said, this is still to be determined and then Eq. ()2.62p is not much 
useful phenomenologically speaking. 

The so called Minimal Hadronic Approximation (MHA) is, in the con- 
text of large- A^c QCD, a much useful phenomenological approximation. Is, 
in itself, a rational approximation which allows us to make the matching 
conditions a useful tool where our knowledge of a certain Green's function 
is centered in its chiral expansion and also the expansion given by its OPE. 
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In that sense, one writes down Eq. ()2.62p with a finite number of states, i.e. 

n=l 

The MHA [651 EHl EH EHl EH] to Large- A^c QCD is an interpolating func- 
tion between these two regimes. These matching conditions imply the fact 
that we have to use two different languages to describe the same physics: 
while at short distances we use quark and gluon as a variables, at long dis- 
tances we use meson fields. These entail that one has to match an expansion 
in powers of Og coming from short distances to an expansion in powers of 
meson momenta at long distances. 

In recent years, a large amount of work has been dedicated to studying 
the consequences of these ideas [70l[7Il[72l[73l[56l[711[Za[76l[50l[771ll9l[l5l 

EE]. 

It was then realized that all the successful results from a chiral resonant 
lagrangian could be encompassed at once as an approximation to large- 
Ac QCD consisting in keeping only a finite (as opposed to the original 
infinite) set of resonances in Green's functions. The main advantage of that 
approximation is, after ensuring that fulfills all the constrains from low and 
high energies imposed, one can go improving the result by adding more and 
more terms in Eq. (j2.7Up expecting a smooth convergence. 
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Chapter 3 

Pade Theory 



3.1 Introduction to Pade Theory 

In the study of mathematical analysis, the relation between a certain func- 
tion and the Taylor coefficients of its expansion is still a profound mathemat- 
ical question. Classically, the answer is that if the Taylor series expansion 
converges absolutely, then defines uniquely the value of a function which is 
differentiable an arbitrary number of times. In practice, however, two prob- 
lems may appear. The first is how to compute the Taylor coefficients of a 
certain function, and the second what would be the range of applicability of 
that Taylor expansion. 

Moreover, when perturbation methods are used to solve a problem, the 
answer usually emerges as an infinite series. If the perturbation series con- 
verge rapidly, summing the few calculated terms gives a good approximation 
to the exact solution. However, it is more common for the series to converge 
slowly, if it converges at all. A regular perturbation series could converge 
slowly when the modulus of its parameter of expansion e is less than the 
radius of convergence. Singular perturbation series diverge for all values of 
e 7^ 0, and even if the series is asymptotic, the value of e may be too large 
to obtain much useful information. 

While discovering that a perturbation series diverges is really discourag- 
ing, it is possible to assign a sensible meaning to the sum of the series and 
even to use the first few terms to approximate this sum. 

Also, a divergent series indicates the presence of singularities and then 
shows the inability of a polynomial to approximate the function adequately. 
Therefore, a summation algorithm is needed which requires as input only a 
finite number of terms of that series. 

On top of these kind of methods, Euler summation and Borel summation 
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are examples of summing an infinite series even when diverges and are been 
used extensively in many fields in physics and biological sciences. 

However, both Euler and Borel summation methods have the difficulty 
to all the terms of the divergent series must be known exactly before the sum 
be done. These cases are, nevertheless, not always realistic enough since, as 
we said, sometimes one can only compute few terms of that series. 

A well-known summation method having this property is called Pade 
summation. 

3.1.1 What is a Pade Approximant? 

The present work will present this method and its properties, and will be 
applied to several physical cases where we will learn aspects of resonance 
saturation, unitarization process, spectral functions, etc. Before going into 
a detailed mathematical explanation we think that it would be appropriate 
to show first the mean advantages of the Pade Approximants (PAs in the 
following) method with a nice and simple example. Later on, we will explain 
in a more rigorous way the main properties of Pade Approximants, especially 
those related with the convergence properties of our approximations. 
Lets consider, then, the following function F{x): 



F(x) = \ = > a„x" = l + + + 0(x^). (3.1 

^ ^ V 1 + x ^^ 2 8 16 128 V ; V ; 

n=0 

As an amusement, lets suppose that the true function F{x) is not known 
but only its Taylor expansion. However, we would like to compute the value 
liuix^oo F{x) (for simplicity, called F(oo)). Since the function F{x) is not 
known (if was, lim^^-^oo F{x) = v2), we will compute it through its Taylor 
expansion. However, since the radius of convergence of that function is 
r = 1/2, for any value of |x| > 1/2 this expansion fails. In this particular 
example we can use a special trick to transform the series into one which will 
let us estimate the value F(oo). Suppose the change of variables x = jr^^ 
then 



n=0 

(3.2) 
For the new variable w, the limit x — )• 00 is translated into w — )• 1/2 and 
now the new Taylor series representation, Eq. (j3.2p . converges aX, w = 1/2. 
The first few successive approximation to -F(oo), namely 'Ylin=Q^nW^ ■, for 
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k 





1 


2 


3 


4 


ELo bnW- 


1 


1.25 


1.34375 


1.38281 


1.39990 



Table 3.1: Results for the lowest partial sums of Eq. (|3.2 



k = 0, 1, 2, 3, 4, • • • are shown in Table 13.11 evaluated at w 
sequence of values are certainly converging to y/2 ~ 1.4241. 



1/2. The 



In fact, the partial sums X]n=o ^"^" ^^^ ^^ written in terms of the 



original variable x following w 
find 



x/{l + 2x). For the first /c = 0, 1, 2, • • • , we 



1, 



1 + (5/2)x 1 + (9/2)x + (43/8)2;2 



(3.3) 



l + 2x ' (l + 2x)2 

which are rational fractions in terms of the original variable x. 

A particular kind of such a rational approximation is the Pade Approxi- 
mant. The mean idea is to match the Taylor series expansion into a rational 
function so that it would tend to a finite limit as x tends to infinity as in 
the case of Eq. (j3.3p . The simplest PA to the function F{x) is: 



R{x) 



1 + bix 



(3.4) 



The three unknown parameters are fixed in the way that the Taylor 
expansion of R{x) be the same Taylor expansion of F{x) in Eq. (|3.ip . We 
find: 



R{x) 



1 + i^ 
l + |rc 



1 + 



5x^ 25x'^ 
"S""^ 32 



125x^ 
128 



+ Oix^). 



(3.5) 



As shown in Eq. (j3.5p . if one expands R{x), finds that the first three 
coefficients are exactly the same that those of Eq. (j3.ip by construction 
but the other coefficients are different. Evaluating now lim^^-i^oo Rix) = 1.4 
which is a better determination that any of the approximations of Table [HTTl 

The next such approximation, which needs 5 parameters, is: 



13, 



41,^2 



i H "^x -\- jqX 



- ^-^ I = 1.41379. 



(3.6) 



Further such approximations converge quite well and one sees an im- 
provement of the result while going further in the approximant sequence. In 
fact, using 7,9, and 11 Taylor coefficients we obtain: 
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1.414201183, 1.414213198 and 1.414213552, 

respectively, the last one is off by only 10^*^. 

We can also form the same type of approximation to the expansion but 
now in powers of the variable w: 

-I w 1 3«) I uP_ 

1 _!^ ±. 1 4 ^ 16 ... ^oyN 

' 1 _ 3h' 1 5w , 5w^' ^'^■' ' 

^ 4 J- ~ — + TT 

Evaluating the approximants at u) = 1/2 we get the results 1, 1.4, 41/29, • • 
identical values obtained using the series expansion in x. 

This invariance principal is a general and important property of Fade 
Approximants and is the basis of their ability to sum the x series in our 
example and give excellent results, even at x — )• oo. 

The successive approximations we have built increase monotonically. Al- 
though this property is not a general one, it can be proven to hold in a wide 
variety of cases. For these cases, one can prove that certain Fade Approxi- 
mants form converging upper and lower bounds. 

Rigorous definition of a Fade Approximant 

Since we already have introduced the notion of a Fade Approximant, lets 
now define them in a more rigorous mathematical way. Let a function f{z) 
have an expansion around the origin of the complex plane of the form 

oo 

n=0 

with a certain radius of convergence r. One defines a Fade Approximant 
to f{z), denoted by Pjf{z), as a ratio of two polynomials Qm{z)-,Rn{z]^ 
of order M and A^ (respectively) in the variable z, with a contact of order 
M + N with the expansion of f{z) around z = 0. This means that, when 
expanding P^ {z) around z = 0, one reproduces exactly the first M-|- A^-|- 1 
coefficients of the expansion for f{z), Eq. 



P^'iz) = f^ « /o+/i z+h z' + ...+fM+N z*^+^+0(z^+*^+i) . (3.9) 

At finite z, the rational function P^(z) constitutes a resummation of 
the series ()3.8p . Is in that sense that sometimes that method is also called 
a Fade summation. 



^Without loss of generality we define, as it is usually done, -Rjv(O) — 1. 
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Of special interest for the present work will be the case when M = N+J, 
for a fixed J. The corresponding PAs P^ "*" (z) belong to what is called the 
near-diagonal sequence for J 7^ 0, with the case J = being the diagonal 
sequence. 

As seen with the example, the mean advantages of that method are that 
the full power series representation of a function need not to be known to 
construct a Pade Approximant -just the first M + A^ + 1 terms. We will see 
that the series can be divergent and still be able to construct a PA sequence 
that approximates quite well the original function. Moreover, sometimes 
they work rather well even beyond their proven range of applicability. 

3.1.2 Examples 

To point out that these advantages are more general that the case we 
have already explained, we would like to show a few well known examples 
of convergence for Pade Approximants in both convergence and divergence 
series representation. 

The convergence series case 

The function e^ has a Taylor series '^'^=qz" /n\, and then can be easily 
transformed into a sequence of Pade Approximant. For all complex z, both 
the Taylor series approximant and the Pade Approximant approach e^ as 
n, M, N — )• 00. The rates of convergence are well known in both cases and 
it can be shown |79j that the relative error done by a Pade Approximant is 
2" times smaller than the error in the Taylor approximant, as n — )• c«. 

Another example can be done using the function l/r(z). In that case, 
although Pade summation enhances the convergence of the Taylor series, the 
effect is not dramatic. The Pade sequence requires about half the number of 
terms that the Taylor series requires to achieve 1% accuracy. Since l/r(z) 
has no singularities, the benefits of Pade summation are marginal. 

The divergence series case 

In these two previous examples, one can not see a compelling reason to 
use Pade summation because the Taylor series already converge for all z and 
the improvement of convergence is not astounding. However, the real power 
of Pade summation arises when applying the method to divergent series. 

As a first example of that case, lets examine the function f{z) = log(l + 
z), which is a Stieltjes function. A function f{z) is called Stieltjes function 
if it can be expressed as a Stieltjes integral according to 
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Here, ^p{t) is a measure on < t < oo and has for all m > finite and 
positive moments Hm defined by 

POO 

^im = / r#(t). (3.11) 

Jo 
The formal series expansion for f{z)^ which do not need to be convergent, 

oo 

/(z) = ^(-l)>™z™ (3.12) 

m,=0 

is called a Stieltjes series if its coefficients ^m are moments of a measure 
■0(i) on < t < oo. According to Eq. (|3.1ip : 

QQ /"OO 

f{z) = Y,i-^rz'^ r#(t). (3.13) 

m=0 •^'^ 

An example of a Stieltjes series is, then, the series for f{z) = log(l + z): 
Log(l + z) = Y, ^ \, (3.14) 

m=0 

For \z\ < 1, the power series converges absolutely, for |z| = 1 it converges 
conditionally, and for \z\ > 1 it diverges. However, as long as the |Arg(2;)| < 
TT, the divergent series can at least in principle be summed. 

The series shown in Eq. p.l4p follows from the integral representation 

■'°'^<i + ^> = ^rrf5- <'^''* 

We only have to set 'ip{t) for < t < 1 and 'ip{t) = 1 for 1 < t < oo. The 
moments fim of this positive measure ip{t) are given by 

/i^ = / rdV(i) = / t'^dt = — . (3.16) 

Jo Jo m + L 

If the positive measure ip{t) corresponding to the Stieltjes function f{z) 
could be determined directly from the Stieltjes moments fim, the value of 
f{z) could at least in principle be computed via the integral representation 
Eq. ()3.10p . even if the series diverges. 

A necessary and sufficient condition which guarantee that a Fade Ap- 
proximant are able to sum a divergent Stieltjes series is that all Hankel 
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determinants D{m, n) formed from the moments fim have to be strictly pos- 
itive, [HD], i.e., 



D{m,n) 



f^m t-^m+l ■ ■ ■ fJ"m+n 

fJ-m+l l^m+2 ■ ■ ■ l^m+n+l 



> . (3.17) 



t^m+n l-l"m+n+l ■ ■ ■ l^m+2n 

Unfortunately, the practical application of this condition is by no means 
simple, in particular if only the numerical values of a finite number of Stielt- 
jes moments Ho^fJ-i, • • • fJ-n are available. However, there is a comparatively 
simple sufficient condition, the so-called Carleman condition |80) . that says 
that if the moments Hm satisfy 



Y^fi^=oo, (3.18) 



m=0 

then the PA sequence P^ ^ converge for J > — 1 to the value of the cor- 
responding Stieltjes function as A^ — )■ oo. It was shown in Ref. [ST] that 
the Carleman condition is satisfied if the moments /Xm do not grow faster 
that C™'^^(2m)! as A^ — )• oo, with C being a suitable positive constant. The 
Carleman condition implies, for example, that the series Eq. ()3.14p can be 
approximated by the Pade sequence P^~^ for J > —1 for all z in the cut 
plane |Arg(2:)| < it, [82j. 

The reason why the Pade Approximants provide a better description 
that the Taylor series for the function f{z) = ln(l + z) is that the poles 
of P^~^ {z) have the ability to approximate the effect of its branch-cut 
singularity. The poles of P^^ (z) and Pj^~ (z) all lie on the real negative 
axis and, as A^ — )• oo, the poles become dense and in some sense resemble 
the branch cut [79| [82l [83] . 

Among other examples of application of Pade summation one can find 
the well know asymptotic series case, for example the Stirling series, which 
can be rapidly approximated by a Pade summation. Or a simple and 
extremely efficient computation of a asymptotic series is the case of the 
parabolic cylinder function, D^(x), which for x > 1 is the solution of the 
parabolic cylinder equation y" + {u + 1/2 + l/4x^)y = 0. 

3.1.3 Convergence theorems 

In despite of these nice examples, since many functions may be asymp- 
totic to the same divergent series, which of these functions, if any, will the 
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Pade sequence select as its limit? Moreover, if the limit function is multi- 
valued, which branch of the function will be singled out as that given by the 
limit of the Pade sequence? There is no general theory of Pade summation 
for arbitrary series yet. Nevertheless, the convergence theory of Pade Ap- 
proximants for the special class of Stieltjes series, meromorphic functionqj, 
continued fractions, hyper geometric functions, Bessel functions is relatively 
completqj. For Stieltjes series, one can partially answer these questions. 
Since a Stieltjes series is the prototype of an asymptotic series 

oo 

^(-l)"n!x" (3.19) 

n=0 

which is a formal solution to a certain differential equation, one can not add 
up all the terms of the divergent series because the sum does not exist. By 
"summing" we mean finding a function to which the series is asymptotic. 
To sum the series ()3.19p we invoke the integral identity n! = L e~^t^dt and 
the interchange X]^o(~^^)" ~ TTH' -^°^ l^^l ^ -'^• 

Y,{-'^Tri\x'^ ^Y.^-xT e-'edt^ dte-'Y.^-xtT (3.20) 

n=0 n=0 -^O -^^ n=0 

The last piece can be solved by splitting the integral in two parts: 

/ dte-*V(-a;t)"= / dte"* V(-xi)" + / die"* V(-xi)". (3.21) 

•^0 n=0 -^0 n=0 "^^ n=0 

Using X]^o(~^^)" ~ T+St ^^^ ^'^^^ integral becomes /g jj^dt. The second 
integral becomes J^ T+^'^^ using the change of variables —xt = 1/y and 
the sum X^.^oy~" ~ «^- '^^^ resulting integral, summing both parts, 
y{x) = Jq°° iTTH^* ^^ ^ Stieltjes integral as Eq. (j3.1Up and satisfies exactly 
the differential equation x'^y" + (1 + 3x)y' + y = 0. 

Finally, integrating y{x) by parts iteratively, one recovers the asymptotic 
series expansion valid as x — )• 0^, Eq. ()3.19p . 

A more interesting case is the so called Generalized Stieltjes integral, 
given by 



^A function is said to be meromorphic when its singularities are only isolated poles. 
^In the next section we will briefly schematize few convergence theorems that will we 
used along these pages. 
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where the weight function p{t) is nonnegative for t > and approaches 
zero so rapidly as t — t- oo that the moment integrals 

oo 

epit)dt (3.23) 



exist for all positive integers n. 

The Pade sequence of a Stieltjes series representation has some remark- 
able convergence properties when its weight function is nonnegative. One 
can show that, for x > 0: 

• The diagonal Pade sequence P^ (x) decreases monotonically as N in- 
creases. 

• The Pade sequence -P^~ (x) increases monotonically as N increases. 

• The sequence P^ (x) has a lower bound, while the sequence Pf^f~ (x) 
has an upper bound. 

These properties imply that limAr_>.oo P^i^) ^^'^ hmiv-!.oo Pj^^ (x) exist and 

also that 

limAT^oo P^~^{x) < limAf^oo Pn{x)- 

Even more, all Stieltjes functions F{x) with the same finite series repre- 
sentation used in the Pade summation satisfy 

lim Piy~\x) < F{x) < lim P^ (x) (3.24) 

These previous inequalities can be generalized as follows. For any J > 
— 1 the Pade sequence P^~^ (x), if generated from a Stieltjes function, is 
monotone increasing when J is odd and monotone decreasing when J is 



even. 

Moreover, since Stieltjes functions are analytic in the cut plane |Arg (x)| < 
TT and since the Pade approximants P^{x) and P^ ~ (x) converge to the orig- 
inal function, all the poles of P^ (x) and Ppf~ (x) must lie on the negative 
real axis. Finally, if a series is a Stieltjes series, then the Pade sequences 
P^{x) and P]\i~ (x) converge on the cut— x plane and the limit functions 
are Stieltjes functions. 

The convergence properties of the PAs to a given function are much 
more difficult than those of normal power series and this is an active field of 
research in Applied Mathematics. As we have alreadu said, those which con- 
cern continued functions, hyperbolic functions and meromorphic functions 
are rather well-known and, specially this last case, will be of particular inter- 
est for this work. There are also several well known theorems that will we of 
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our interest, among them Montessus de Ballore and the Nuttal-Pomerenke 
theorem. The fohowing section presents these theorems without so much 
detail but accurately, since these theorems will we used along these pages 
and we will refer it several times. 

Montessus de Ballore'' s theorem applies to a function f{z) which is mero- 
morphic with its M poles within a circle. In practical terms, provided M is 
known, Montessus' theorem asserts uniform convergence of the sequence of 
Pm approximants within the circle. The complete statement of Montessus' 
theorem [M] when the function has simple poles is the followingfj: 

• Let f{z) be a function which is meromorphic in the disk \z\ < R, with 
precisely M simple poles at distinct points zi,Z2,- " ■, ^M, where 

< \zi\ < \z2\ < ■■■ < \zm\ < R- 

Then 

lim Pt = f{z) (3.25) 

uniformly on any compact subset of 

TiM = {z, \z\ < R,z^ z„i = l,2,--- , M} . 

Later on, R. de Montessus de Ballore extended the theorem for the pos- 
sibility of multiple poles instead of simple poles: 

• Let f{z) be a function which is meromorphic in the disk \z\ < R, with 
m poles at distinct points zi,Z2,- ■ ■ z^ with 

0<\zi\<\z2\<--- <\Zm\ <R- (3.26) 

Let the poles at z^ have multiplicity fk, and let the total multiplicity XlfcLi ^fc - 
M precisely. Then 

lim Ptj = fiz) (3.27) 

uniformly on any compact subset of 

TIm = {z, \z\ < R,z^ Zk,k = l,2,--- , m}. 

The main advantage of that theorem is that guarantees uniform conver- 
gence of the approximation within a circle defined by \Zfn\ < R- No spurious 
poles coming from the approximation will be found inside this disk. Also, 
the position of the M poles will be correctly determined by taking the limit 



An easy proof of the theorem can be found in [80], chapter 6.2. 
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limi_s.oo in Eq. (|3.27p . given the number of poles and multiplicities M is 
known. An interesting application of that theorem can be found in Ref. |85j . 

In despite of these convergence statements, it is difficult to go beyond 
that since there are no proven rate of convergence for the Pj^ sequence to 
the original function, or upper and lower boundary constrains as it happens 
when one constructs the sequences P^ and -P/y^ for Stieltjes functions. 

When the function is Stieltjes and meromorphic one has to choose which 
theorem would be better to use. It will depend on the underlying motivation 
to use Montessus' theorem or the theorem for Stieltjes functions that we have 
already described. The main difference between both convergence theorems 
is that the Montessus' theorem guarantees convergence uniformly inside a 
certain disk while the other one provides a boundary constrains for the 
approximation . 

When the function is meromorphic but is not Stieltjes, the Montessus' 
theorem says nothing about upper and lower boundary constrains. Then, in 
practice, one constructs a Fade sequence and estimates the accuracy of the 
method by looking into the difference between consecutive PAs, Refs.jSUj 

An extension of the Montessus de Ballore theorem is the Nuttall-Pommerenke's 
Theorem [89] (commonly known as Pommerenke's Theorem) which asserts 
that the sequence of (near) diagonal PA's to a meromorphic function is 
convergent everywhere in any compact set of the complex plane except, per- 
haps, in a set of zero area. This set obviously includes the set of poles 
where the original function f{z) is clearly ill-defined but there may be some 
other extraneous poles as well. For a given compact region in the complex 
plane, the Nuttall-Pommerenke's theorem of convergence requires that, ei- 
ther these extraneous poles move very far away from the origen as the order 
of the Fade Approximant increases, or they pair up with a close-by zero be- 
coming what is called a defect or a Froissart doublevQ in the mathematical 
jargon [90]. These are to be considered artifacts of the approximation. Near 
the location of these extraneous poles the FA approximation clearly breaks 
down but, away from these poles, the approximation is safe. 

To extend the Montessus' theorem one has to consider that the circle 



^At thirties, M. Froissart observed that the Pade approximants to a power series per- 
turbed by random noise are characterized by some zeros and poles unusually located in 
the vicinity of the unit circle (which was the range of convergence of the particular series 
he used). Each such zero was accompanied by a pole at the distance proportional to the 
scale of the noise, forming the so-called Froissart doublets. Since then, several studies are 
been done and now this defect phenomena is better understood. Actually, for Stieltjes 
functions the location Zd of these defects have to obey \zd\ > R where R is the convergence 
radius of its Taylor expansion around z — 0. 
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of convergence of a P^j to a certain meromorphic function f{z) can grow 
provided exist M' > M, M' been the poles of the extended subset 3a/'- 
Then, the first step (a weaker form of Montessus' theorem) is a theorem 
apphcable when the degree of the denominator is known to be grater than 
or equal to (instead of precisely equal to) the number of poles within the 
circle of convergence of the P^ . 

Let F{z) be analytic at the origin and also in a given disk \z\ < R, except 
for m pols, counting multiplicity. Consider a sequence of P^ to f{z) with M 
fixed, M > m, and L — )• oo. Suppose that arbitrarily small, positive e and 
5 are given. Then Lq exists such that \f{z) — P^\ < e for any L > Lq and 
for all \z\ < R except for z £ ei, where e^ is a set of points in the z-plane of 
measure less than 6. 

An interesting corollary is deduced at this point: with the hypothesis of 
the theorem, the more general PA sequence P^/ satisfies \f{z) — Pj^l \ < e 
for any k > k^ and for all \z\ < R excepting z S e^, e^ a measure less than 
6, provided 

• Lk/Mk — )• oo as fc — )• oo {Mk ^ 0), and 

• Mk>M for all k > /cq. 

An interesting case consequence of that corollary is when M^ = Lk, 
called the diagonal Fade sequence. 

This last comment is a weak form of both what is known to be true 

and what is expected to be true about convergence in measure of the Fade 

^L 

M- 



sequence Pj^j. Nonetheless, it provides a bass for further development. 



First, the diagonal sequence may be replaced by the sequence P^/ , k = 
1,2,... provided that, for any A in the range < A < 1 however small, 
A < jf- < X"^. Frovided that inequality holds and L^ + Mk — >• oo, this 
weaker constraint is sufficient to allow convergence in measure. 

Second, f{z) need not be meromorphic, but may also have a count- 
able number of isolated essential singularities. This means that for example 
f{z) = exp(— (1 — z)~^) and g{z) = exp(2:r(z)) are allowable functions, but 
not function whose singularities have a limit point in the finite z-plane. 

All these steps leads us to the Fommerenke's theorem, which reads: 

• Let f{z) be a function which is analytic at the origin and analytic in the 
entire z— plane except for a countable number of isolated poles and essential 
singularities Then the sequence of approximants P^ with L/M = A(A 7^ 



The inclusion of essential singularities was a notable step forward in the extension of 
the Nuttall's theorem [91] and in the global convergence theory of Pade Approximants. 
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0, A 7^ oo) satisfies: 



lim P{;f = f{z) (3.28) 

on any compact set of the z— plane except for a set of points of zero 
measure. 

An interesting corollary of this theorem [80j says that it can be gener- 
alized for Pade sequences of the type P^~^ , for a fixed k, and then also 
convergence for these sequences is found when N ^ oo. 

The main advantage of the application of this theorem in front of Montes- 
sus' theorem is now one does not need to know in advance the number of 
poles of the original function inside a certain disk. The disadvantage, how- 



ever, is when you want to construct a sequence of -Pjy''' > ^° increase one 
order on the sequence, one needs, at least, two different inputs (one for the 
numerator and one for the denominator of the PA) while following Montes- 
sus' theorem, increase one order in the sequence only demands one new 
input (since the denominator of the PA is fixed to a certain degree M). 

As N ^ oo, the Pommerenke's theorem ensures convergence of the se- 
quence of PAs to the original meromorphic function, in any compact set in 
the complex z plane except at a finite number of poles. Of course, where 
there is convergence, the PA may be considered an approximate resumma- 
tion of the Taylor series around the origin. On the other hand, the set of 
points where there is no convergence certainly includes the position of the 
poles since not even the original function is defined there, but, as we said, 
there may appear other artificial poles which have no counterpart in the 
original function (when, for example, k < 0). One would naively think that 
the presence of these artificial poles would cause a major distortion and com- 
pletely spoil the rational approximation. However, one can actually show, 
as the theorem states, that as the order of the Pade increases, i.e. as N 
grows, these artificial poles either move to infinity in the complex plane and 
decouple or they get "almost-canceled" by the appearance of nearby zeros, 
d la Froissart. Although, in general, this cancelation is not complete, it is 
efficient enough to make the region of distortion of the artificial pole only of 
zero measure. This is why and how the Pade Approximation works. For an 
explicit example where all these properties come to play in the context of a 
Regge-inspired model, we refer to Ref. [86]. 
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3.1.4 Generalization and Extensions of Pade Approximants 

Pade Types and Partial Pade Approximants 

The last theorems implies somehow certain knowledge of the position of 
the poles of the original function. Certainly, as L — )• oo, the M poles of Eqs. 
()3.25p . ()3.27p and ()3.28p will tend to reproduce the original ones. When the 
position of the poles in the original function is known (for example, for a 
Green's Function, at least the lowest lying states), it is interesting to devise 
a rational approximation which has this information already built in. The 
corresponding approximants are called Partial Pade Approximants (PPAs) 
in the mathematical literature [92j and are given by a rational function 

" -(z) = — ^^^^Vt' (3-29) 



jM / X _ Qm{z) 

^'^^"^ - Rm{z) Tk{z) 



where Qm{z)-,Rn{z) and Tk{z) are polynomials of order M,N and K (re- 
spectively) in the variable z. The polynomial Tk{z) is defined by having 
K zeros precisely at the location of the lowest lying poles of the original 
functioiu i.e. 

Tk{z) = {z + zi) {z + Z2) ... {z + zk) . (3.30) 

As before the polynomial Rn{z) is chosen so that i?Ar(0) = 1 and, together 
with Qm{z), they are defined so that the expansion of F^ j^{z) matches 
exactly the first M + A^ + 1 terms in the expansion of the original function 
around z = 0, i.e. : 

KM^) ^fo + fiz + hz' + ... + fM+N ^^'+^ + 0(z^+*^+i) . (3.31) 

At infinity, ¥^ ^{z) obviously falls off like l/z'^^^^^'^ . Exactly as it happens 
in the case of PAs, also the PPAs will have defects for a general meromorphic 
function, the Froissart doblets. 

Finally, another rational approximant defined in mathematics is the so- 
called Pade Type Approximant (PTA) [92] T^{z) : 

TfJiz) = ^ , (3.32) 

where T]\f{z) is also given by the polynomial (|3.30p . now with A^ preassigned 
zeros at the corresponding position of the poles of the original function, f{z). 
The polynomial Qm{z) is defined so that the expansion of the PTA around 



''For simplicity, we will assume that all the poles are simple. 
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z = agrees with that of the original function up to and including terms of 
order M + 1, i.e. 

r^{z) ^fo + fiz + f2Z^ + ... + fM z'' + 0(z*^+i) . (3.33) 

At large values of z, one has that T^{z) falls off like l/z^"^"'^. Clearly the 
PTAs are a particular case of the PPAs, i.e. T^{z) = ¥q^^{z) and coincide 
with what has been called the Hadronic Approximation to large- A'c QCD in 
the literature [651 [Ml ESI EH EZ] . 

Since in that last case all the poles are preassigned precisely in the loca- 
tion of the lowest lying poles of the original function, Froissart doblets will 
not we found. 

We may not confuse Pade Types with the Montessus' theorem state- 
ments. The Pade Type sequence converges to the original function when 
both M, A^ — )• c«. In the Montessus' theorem instead, the denominator was 
fixed to a certain value and only the degree L of the numerator was enlarge 
to reach convergence (Eq. (|3.27|) ). 

N-point Fade Approximants 

The Pade method that we have already introduced could be called One- 
point Pade method because the approximants are constructed by matching 
them with a power series about a particular point zq. However, the function 
in question may have been investigated in the vicinity of two or more points. 
One may wish to incorporate information from all these expansions in a 
single sequence of Pade approximants. Suppose /(z) has the asymptotic 
expansions 

CX) 

f{z) = J2 an{z - zor, Z ^ Zo, (3.34) 

n=0 

oo 

f{z) = Y,bn{z-z,r, z^zi, (3.35) 

n=0 

in the neighborhoods of the distinct points zq and zi, respectively. A Two- 
point Pade approximanllfl to f{z) is a rational function F{z) = Rn{z)/Sm{z) 
where 5^/(0) = 1. Rn{z) and Sm{z) are polynomials of degree N and M, 
respectively, whose (A + M + 1) arbitrary coefficients are chosen to make 
the first J terms (0 < J < N + M + 1) of the Taylor series expansion of F{z) 
about Zq agree with Eq. ()3.34p and the first K terms of the Taylor series 



We are basically interested in Two-Point Fades since they are the closest version 
of that method commonly used in the physical literature and called Minimal Hadronic 
Approximation. 
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expansion oi F{z) about zi agree with Eq. (|3.35p . where K+J = N + M + \. 
This idea can be extended to a A^— point Pade when the rational function 
fits various points, not necessary distinct. Here we illustrate the method by 
one Two-point Pade example [79j. Define the function f{z) as 

This function is the unique solution to the differential equation 2zf'{z) = 
~(1 + 2z)f{z) + 1. This solution has power series expansions about both 
z = and z — ?• oo. The series expansion are: 



^ C-4)"n' 

/(z) = J^a„z", an = ^^^ I •, , |z|^0(.3.37) 

n=0 



oo 



f{z)^}_^bnz~^, 6„ = 2(2n-2)!4-"/(n-l)!, 60 = 0, |z| ^ oc(.3.38) 

n=0 

For definiteness we discuss only the diagonal Pade sequence with N = M 
and use as input J = N + 1 and K = N. For this particular function, since 
60 = 0, all the PAs constructed with the expansion z — )■ 00 will have the 
same behavior as 1/z for z — )• 00 as the function f{z) has. 

Two-point Pade Approximants are always more constraint than One- 
point PAs. One can compare the diagonal sequence One-point Pade ap- 
proximants obtained from the expansions at z = and 2; — )• 00 with the 
Two-point Pade from both expansions. In particular, the lowest One-point 
PA constructed with the expansion z = reads 



1 _ 4z 

i + ¥' 



(3.39) 



5 
while the lowest One-point PA constructed with the expansion z — )• 00 reads: 

Both PAs needs three inputs. To construct the lowest Two-point PA we 
will also use three inputs in such a way that the PA will mucfo one order 
from the expansion z — >• 00 and two orders from z = 0. That PA reads: 

(3.41) 



^ ^ 3 



^All the PAs we are constructing are diagonal and then we need an odd number of 
inputs. In that situation, when dealing with a Two-point PA one can always choose 
K < J or K > J but never K = J. 
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PA 


z = 0.1 


z= 1 


z= 10 


z = 100 


FA^-o Ea.(l3.39p 


1.00 


0.97 


-6.30 


-124 


PA^^oo -h^q.dMOj) 


-1.34 


1.86 


0.99 


1.00 


2point-PA Eq.dMU) 


1.00 


1.12 


2.46 


2.94 



Table 3.2: Comparison between a One-point PA constructed using the Tay- 
lor expansion (first raw), One-point PA using the expansion at z — )• c« (sec- 
ond raw), and a Two-point PA (third raw). All the numbers are normalized 
to the function /(z). 



PA 


z = 0.1 


z = 1 


z= 10 


z= 100 


PA,-o Eq.(l3i2j) 


1.00 


1.00 


3.46 


92.8 


PA.^oo Eq.(l3i3D 


-2.78 


1.12 


1.00 


1.00 


2point-PA Eq.(|3.44p 


1.00 


1.02 


1.23 


1.04 



Table 3.3: Comparison between a One-point PA constructed using the Tay- 
lor expansion (first raw). One-point PA using the expansion at z — )• oo (sec- 
ond raw), and a Two-point PA (third raw). All the numbers are normalized 
to the function f{z). 

The Two-point approximant gives a more uniform approximation to /(z) 
while in the neighborhood of each point the corresponding One-point Pade 
is slightly more accurate as can be see numerically in the table [3^2] (we show 
the values for the ratio PA//(2;)). 

We can easily improve these results by going to the next PA in the 
corresponding sequence. In the three cases we will need 5 inputs and the 
One-point PA constructed using expansion 2 = will be: 



1 



2z_ _|_ 322^ 
9 "*" 945 
1 I 4£ I 42£ 
-^ "•" 9 "•" 63 



while the One-point PA constructed using the expansion around z 
be: 

-5 I 22 
3 "•" 3 



1 



4z + ^ 



In its turn, the Two-point PA will be: 



1+22 
^ ^ 15 



1_L 42 , 42^ 

-^ ^ 5 ^ 15 



(3.42) 
oo will 

(3.43) 
(3.44) 



In that second case, the numerical results are shown in table 13.31 
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For small z the Two-point Fade is significantly more accurate that the 
ordinary One-point Fade about z — t- oo. For large z, the Two-point Fade is 
significantly more accurate that the One-point Fade from z = 0. One has 
to consider that in those cases a Two-point Fade of degree N is comparable 
in accuracy than the One-point Fade of degree iV/2. 
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Chapter 4 

Pade Theory and 
Meromorphic functions 



In the second chapter we have discussed the main properties of the Pade 
Theory. We have summarized the main convergence theorems we will use 
in these pages and we have shown few examples of the procedure to follow. 
Now comes the turn of the physical case. 



4.1 Rational Approximations: generalities 

In the physical case the original function f{z) will be a Green's function 
G{Q'^) of the momentum variable Q^. In general, the structure of that kind 
of function is not know. However, in QCD in the large Nc limit this Green's 
function is meromorphic with all its poles located on the negative real axis in 
the complex Q^ plane. These poles are identified with the meson masses. On 
the other hand, the region to be approximated by the Pade Approximants 
will be that of euclidean values for the momentum, i.e. Q^ > 0. The 
expansion of G{Q'^) for Q^ large and positive coincides with the Operator 
Product Expansion (section [2. 3p . 

In general a meromorphic function does not obey any positivity con- 
straints and, as we will see, this has as a consequence that some of the 
poles and residues of the PAs may become comp/eaiJ- This clearly precludes 
any possibility that these poles and residues may have anything to do with 



^On the contrary, remember the case of a Stieljes functions which obeys positivity 
constrains. In that case the poles and residues of the diagonal and paradiagonal PA 
sequences are purely real and with the same sign as those of the original function |93l I82| . 
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the physical meson masses and decay constants. However, and this is very 
important to reaUze, this does not spoil the validity of the rational approx- 
imation provided the poles, complex or not, are not in the region of Q^ 
one is interested in. It is to be considered rather as the price to pay for 
using a rational function, which has only a finite number of poles, as an 
approximation to a meromorphic function with an infinite set of poles. 

Given the meromorphic function G{Q'^) = Q^Yly-Ai—Q'^) in the Q^ 
complex plane with an analytic expansion around the origin, as in Eq. (|4.7|) . 
it is possible to construct a Pade Approximant, P^ {Q"^), such that its ex- 
pansion in powers of Q^ matches that of the original function up to, and 
including, the term of 0{Q'^^^^^'). Since the function falls off at large Q^ 
as Q~^ up to logarithms (see Eq. ()4.5|) ). we choosci A^ = M -|- 2 in order 
to optimize the matching of the rational approximant at large Q^ to this 
behavioo We emphasize, however, that this choice does not affect the prop- 
erties of convergence of Pade Approximants, as described in section 13.1.41 
Later on we will show that this procedure is the most efficient way to deal 
with what we are interested in. However, this is not a mathematical theo- 
rem and it depends on the analytical properties of the function and also the 
motivation of doing PAs. To study Taylor coefficients, poles and residues 
this strategy is certainly the best option to follow however you have to treat 
with the Proissart doublets of the approximation. If one is interested in 
take under control the analytic properties of the PAs under a certain disk 
on the Q^ variable, one should proceed for example following Montessus^ 
theorem, where the dominator of the approximant is fixed in a certain way 
(sec. I3.1.3P and then keep it under control. On the contrary, since we are 
interested in Taylor and OPE expansions for G{Q'^) = Q^IIv-a{—Q'^), we 
will use Pommerenke^s theorem (sec. l3.1iB1) in this chapter. 

Pommerenke's theorem (Eq. (|3.28|) ) is important in this case because it 
teaches us useful information about the qualitative behavior of how Pades 
approximate meromorphic functionqfj. Regretfully, when asking more quan- 
titative questions such as the rate of convergence, which is the first step 
towards an estimate of the error, such a theorem is only of limited practical 
importance. In practice, one can take a more useful approach towards an 
estimate of the error by studying the behavior of a set of successive rational 



■^This choice is not arbitrary as we will see in section H31 

^Due to the presence of logarithms in Ea. (|4.5p . however, this matching cannot be 

perfect. 

^Remember that Pommerenke's theorem can be applied to any meromorphic function 



which need not to be Stieltjes. 
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approximants, as we will show later on. 



4.2 Testing rational approximations: a model 



4.2.1 The model 

A particular model should help on the understanding of all the PA prop- 
erties presented in the previous chapter. To reassemble a physical case we 
shall consider a correlation function of a vector and axial- vector currents, 
which is particularly sensitive to properties of Spontaneous Chiral Symmetry 
Breaking, sec. 12.21 namely the two-point function 

^^'J'iQ) = i j d'xe^^^J^^^H^)Jl ^'^(0)) , (4.1) 

with Jy{x) = d{x)^^u{x) and J^(x) = d{x)'y^j^u{x). 

In the chiral limit where the light quark masses are set to zerqj, this 
two-point function only depends on one invariant function (Q^ = —q'^ > 
for q'^ spacelike) which is 

Kl'^iQ) = {%<lu - g^uq^) ny,A(g') . (4.2) 

As it is known, the difference H-v-A^ff') defined as follows, satisfies the 
unsubtracted dispersion relation given bjo 



1 f^ dt 1 

ny-A(g') = 77(ny('?') - nA(g')) = llm / - luiIiy^Ait) . 

I A-s>oo Jo t — q"^ — l€ TT 

(4.3) 
and is usually called IIlr{Q'^), the self-energy function, and its proper- 
ties are rather well known in the limit of vanishing light quark masses. First, 
Ilijl{Q'^) vanishes order by order in perturbation theory and is an order pa- 
rameter of the spontaneous breakdown of chiral symmetry for all values of 
the momentum transfer. It also governs the electromagnetic tt~^ — tt^ mass 
difference [95] 

^1+\eM = -V72 / dQ' Q'^LRiQ') ■ (4.4) 

^ o/tt Jo 



^See section [2^ for details. 

^The upper cutofT which is needed to render the dispersive integrals mathematically 

well defined can be sent to infinity provided it respects chiral symmetry [94) . 
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This integral converges in the ultraviolet region as shown by Shifman, 
Vainshtein and Zakharov [27j using Wilson's [28] Operator Product Expan- 
sIoelI 



lim -Q^Ulr{Q^) = ^ [Svr^ (^ + ©(a^ log{Q^))) {^^ 



(4.5) 
where Og is the QCD coupling constant and (?/''0) is the quark condensate. 

Unlike the expansion around the origin of energies, the existence of non- 
vanishing anomalous dimensions, even in the large— A'c limit, gives rise to 
the logQ^ terms and, unlike ()4.7p . renders the expansion around infinity in 
(|4.5p not analytic. 

Witten [961 197] has furthermore shown that 

- Q^IIlr{Q^) > for 0<g2<oo, (4.6) 

which in particular ensures the positivity of the integral in Eq. ()4.4p and thus 
the stability of the QCD vacuum with respect to small perturbations induced 
by electromagnetic interactions. This positivity is also seen experimentally 

M- 

The low Q^ behavior of this self-energy function [981 199] , is governed by 
chiral perturbation theory 

- Q'niH(Q') = In + ^LwQ^ + SCsrQ^ + O(Q^) , (4.7) 

where Liq is the only LEG |17j, sec. 12.21 with an axial content in the 
0{p'^) low energy effective chiral Lagrangian, Eq. ()2.38p . and Csr is the 
corresponding LEG at 0(p^) pnjin]. 

In the large-A'c limit, the spectral function associated with IIlr{Q'^) 
consists of the difference of an infinite number of narrow vector states and an 
infinite number of narrow axial-vector states, together with the Goldstone 
pole of the pion [TOOlfTOT] : 



^ oo 

-lmny(t) = 2F^6{t-M^) + 2Y^F^{n)6{t-M^{n)) , 

n=0 
^ oo 

- Imn^(t) = 2F^6{t) + 2 y Fl{n)5{t - Ml{n)) . (4.8) 

n=0 

Here Fp, Mp are the electromagnetic decay constant and mass of the 



p meson and Fv^A{n) are the electromagnetic decay constants of the n' 



th 



^To obtain that expression one can follow the procedure shown in section E3] and using 
the Weinberg Sum Rules (WSR), App. [H 
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resonance in the vector (resp. axial) channels, while My^Ain) are the cor- 
responding masses. Fq is the pion decay constant in the chiral limit. The 
dependence on the resonance excitation number n is the following: 

Fyj^{n) = F^ = constant , Myj^{n) = niyj^ + n A^ , (4.9) 

in accord with known properties of the large- A'c limit of QCD |30[ [3T] as 
well as alleged properties of the associated Regge theory [1021 1103] I104J . 
The combination (j4.3p thus reads: 



n=0 



nLR(g') ^2 + _q2 + M2 + ^ \ -g2 + M^{n) -q^ + M^ (n) 



(4.10) 



or 



(4.11) 

where Ff = ffMf and — g^ = Q^. This last expression will we useful when 
studying OPE and Taylor expansions, Eq. (|4.5|) and Eq. (|4.7|) respectively. 
This two-point function can be expressed in terms of the Digamma func- 
tion ijj{z) = ^logr(z) as [M] 

(4.12) 
To resemble the case of QCD, we will demand that the usual parton- 
model logarithm is reproduced in both vector and axial- vector channels and 
that the difference ()4.3p has an OPE which starts at dimension six ()4.5p (or 
at dimension four considering Q^^lr{Q^))- A set of parameters satisfying 
these conditions is given bjo 

Fq = 85.8 MeV , Fp = 133.884 MeV , F = 143.758 MeV , (4.13) 

Mp = 0.767 GeV, uia = 1.182 GeV, my = 1.494 GeV , A = 1.277 GeV , 

This set of parameters has been chosen to resemble those of the real world, 
while keeping the model at a manageable level. For instance, the values of 
Fp and Mp in Eq. ()4.13p are chosen so that the function I\.lr in (|4.12p has 



*These numbers have been rounded off for ttie purpose of presentation. Some of the 
exercises which will follow require much more precision than the one shown here. 
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vanishing l/Q^ and l/Q^ in the OPE at large Q^ > 0, the Weinberg Sum 
Rules, as in real QCD (see App. |B]for more details). Also, the model ad- 
mits the introduction of finite widths (which is a l/A'c effect) in the manner 
described in Ref. |100j . after which the spectral function looks reasonably 
similar to the experimental spectral function. This comparison can be found 
in Fig. 5 of Ref. |101j . But this model is also interesting for a very differ- 
ent reason. In Ref. [lUlj several attempts were made at determining the 
coefficients of the OPE by using the methods which have become common 
practice in the literature. Among those we may list: Finite Energy Sum 
Rules [lOHl dUSl UnZl, with pinched weights [Ml |109], Laplace sum rules 
[110| and MHA [651 ESI EZ]- As it turned out, when these methods were 
tested on the model, none of them was able to produce very accurate re- 
sults. We think that this makes the model very interesting as a way to assess 
systematic errors |lllj . 

Defining the expansion of the Green's function ()4.1ip in Q^ = —q^ 
around Q^ = 0, Q^ — ?• oo as 

Q2niR(-Q2)~^C'2fcg2'= , with A: = 0,±l,±2,±3,... (4.14) 

k 

one obtains that the coefficients accompanying inverse powers of momentum, 
akin to the OPE at large Q^ > 0, are given by {p = 1, 2, 3, ... with k = 1 —p): 



C2k — —Fn Spl + 



(_l)p+i 



p2j^2p-2 _ }_p2^2p-2 f o (^^\ _ R (^^ 



(4.15) 



where Bp{x) are the Bernoulli polynomials [112J. As stated above, Fp and 
Mp are defined by the condition that the above expression Eq. ()4.15p van- 
ishes for k = 0,-1 (in agreement with Eq. ()4.5p above), enforcing that 
<5^ IIlr{—Q'^) ~ Q~^ at large momentum as in QCD, Eq. (|4.5p . We em- 
phasize that the above coefficients of the OPE ()4.15p can not be calculated 
by a naive expansion at large Q^ of the Green's function ()4.10p . In other 
words, physical masses and decay constants do not satisfy the naive WSR 



On the other hand, for the coefficients accompanying nonnegative pow- 
ers of momentum, akin to the chiral expansion at small Q'^, one has {k = 
1,2,3,...): 



m,. \ ,fi,_u / m 



C0--F0 , C2k-{-i) Ml^-jk^Yy,j2k{^' [-7^)-^ VA2 



(4.16) 
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Co 


C2 


C4 


c. 


C_4 


C_6 


C-8 


-7.362 


21.01 


-43.92 


81.81 


-2.592 


1.674 


-0.577 



Table 4.1: Values of the coefficients C2k from the high- and low-Q^ expan- 
sions of Q2 IIlr{-Q'^) in Eq. (HT^ in units of lO'^ GeV"^-"^^ . Notice that 
C-2 = and Cq = —Fq (the pion decay constant in the chiral limit), see 
text for details. 

where ip^ '{z) = d il){z)/dz . In Table WA\ we collect the values for 
the first few of these coefficients C2k foi' both expansions Eq. ()4.15p and 
Eq. gH}. 

Following the definitions in section 13. H let us start with the construction 
of the rational approximants to the function Q^ ^lr{—Q^) using as input 
the coefficients of its chiral expansion (j4.16p . Since our original function 
Eq. (j4.12p falls off at large Q^ as Q"^, this is a constraint we will impose on 
all our rational approximants by selecting the appropriate difference in the 
order of the polynomials in the numerator and denominator. 



4.2.2 Construction of Pade Approximants 

The simplest PA satisfying the right falloff at large momentum is P2 (Q^)) 
so we will begin with this case. In order to simplify the results, and unless 
explicitly stated otherwise, we will assume that dimensionful quantities are 
expressed in units of GeV to the appropriate power. Fixing the three un- 
knowns with the first three coefficients of the chiral expansion, Eq. (j4.16p . 
for the function defined in Eq. (j4.12p (i.e. 6*0,2,4)) one gets the following 
rational function: 



P'^iQ' 



R 



(Q2 + ^ii)(Q2 + 4) 



'R 



3.379X 10"^ , ZR = 0.6550+i 0.1732 

(4.17) 

We can hardly overemphasize the striking appearance of a pair of complex- 
conjugate poles on the Minkowski side of the complex Q^ plane. Obviously, 
this means that these poles cannot be interpreted in any way as the meson 
states appearing in the physical spectrum ()4.8l I4.10p . In spite of this, if 
one expands ()4.17p for large values of Q^ > 0, one finds C_4 = — r^ = 
—3.379 X 10^^ which is not such a bad approximation for this coefficient of 
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Im(q") 
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10 
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Figure 4^.1: Location of the poles of the Fade Approximants Fi2{—q)(blue), 
P32 (— q^)('(7reenj and F^2{~1^)(''~^d.) ^"^ ^^^ complex q^ plane. When the order of the 
FA is increased, the overall shape of the figure does not change but the two branches 
of complex poles move toward the right, away from the origin. 



the OPE, see Table 14.11 Even better is the prediction of the fourth term in 
the chiral expansion, which is Cq = 79.58 x 10^'^. 

This agreement is not a numerical coincidence and the approximation 
can be systematically improved if more terms of the chiral expansion are 
known. In order to exemplify this, we have amused ourselves by constructing 
the high-orders PAs P^^iQ'^), P^iQ"^) and P^^iQ"^) (see FigSl]), always 
keeping in mind the appropriate falloff. The particular case of -P52 (Q^) 
correctly determines the values for C_4__6,-8 with 52, 48 and 45 decimal 
figures (respectively). In the case of C206) which is the first predictable 
term from the chiral expansion for this Pade Approximant, the accuracy 
reaches some staggering 192 decimal figures. This is all in agreement with 
Pommerenke's theorem [89], sec. 13.1.31 

As it happens for the PA ()4.17p , also higher-order PAs may develop some 
artificial poles. In particular. Figure HTJ shows the locations of the 12-(blue), 
32-(green) and 52-(red) poles of the PAs P^^iQ^), PEiQ"^) and P^^iQ"^) in 
the complex q^ plane, respectively. As can be seen in that figure, Pi2{Q'^) 
has 6 poles which are real and close to the lowest part of the spectrum of the 
model, (|4.10ti4.13|) . and the other 6 poles are complex conjugate pares. At 
the same time, -P32 (Q^) ^as 16 real poles and P^2{Q'^) ^as 26. Notice that 
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Figure 4-2: Location of the poles (dots) and zeros (squares) of the Fade Approx- 
imant F^2i~Q^) *'^ ^'^c complex q^ plane. We recall that Q^ = —q^. Notice how 
zeros and poles approximately coincide in the region which is farthest away from 
the origin. 

the lowest real poles of the three PAs shown in Fig l4.1l are overlapped and 
that the complex conjugate pares move away (to infinity) when increasing 
the order of the Pade approximant. 

Figure HT2] also shows the locations of the 52 poles of -P52 (Q^) but now 
including the location of its zeros (blue squares). In that higher case, the 
first 26 poles are purely real and the rest are complex-conjugate pairs. A 
detailed numerical analysis reveals that the poles and residues reproduce 
very well the value of the meson masses and decay constants for the lowest 
part of the physical spectrum of the model given in Eqs. (|4.1Uti4.l5|) . but 
the agreement deteriorates very quickly as one gets farther away from the 
origin, eventually becoming the complex numbers seen in Fig. 14.21 It is by 
creating these analytic defects that rational functions can effectively mimic 
with a finite number of poles the infinite tower of poles present in the original 
function (fiJ2]) . 

For instance the values of the first pole and residue in -P52 (Q^) reproduce 
those of the p in Eq. ()4.13p within 193 astonishing decimal places. However, 
in the case of the 26th pole, which is the last one still purely real, its location 
agrees with the physical mass only with 3 decimal figures. This is not to be 
considered as a success, however, because after the previous accuracy, this 
is quite a dramatic drop. In fact, the residue associated with this 26*^ pole 
comes out to be 29 times the true value. The lesson we would like to draw 
from this exercise should be clear: the determination of decay constants and 
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masses extracted as the residues and poles of a PA deteriorate very quickly 
as one moves away from the origin. There is no reason why the last poles 
and residues in the PA are to be anywhere near their physical counterparts 
and their identification with the particle's mass and decay constant should 
be considered unreliable. Clearly, this particularly affects low-order PAs, for 
example that of Eq. ()4.17p . 

A very good accuracy can also be obtained in the determination of global 
euclidean observables such as integrals of the Green's function over the in- 
terval < Q^ < oo. Notice that the region where one approximates the true 
function is far away from the artificial poles in the PA. For instance, one 
may consider the value for the integral 

In = (-1) / dQ2 Q^Ulb{Q'^) = 4.78719 x 10-^ (4.18) 

Jo 

which, up to a constant, would yield the electromagetic pion mass difference 
in the chiral hmit [TT3] in the model (jil^ . The PA -Pja (Q^) reproduces 
the value for this integral with more than 42 decimal figures. This suggests 
that one may use the integral ()4.18p as a further input to construct a PA. 

For example if we fix the three unknowns in the PA P2 (Q^) by matching 
the first two terms from the chiral expansion but now we complete it with 
the pion mass difference (j4.18p instead of a third term (C4) from the chiral 
expansion as we did in (j4.17p . the approximant results to be 

_ ^2 
P^{Q^) = —^- ^ , with rl = 2.898X 10-3 , zr = 0.5618+i 0.2795 . 

(4.19) 
This determines C_4 = -2.898 x 10"^ and C4 = -41.26 x 10"^, which 
shows that using the pion mass difference is not a bad idea (see table 14.21 
for comparison). Notice how the position of the artificial pole has changed 
with respect to ()4.17p . We remark that this last procedure, although reason- 
able from the phenomenological point of view, strictly speaking lies outside 
the standard mathematical theory of rational approximants [SU\ I92j . and 
while we will still call -P2 (Q^) ^ Pade Approximant, this is not a strict PA 
anymore. In section 14.41 we will see how crucial is go beyond the classical 
Pade procedure and we able to control the impact of the modification we 
are introducing here. 



4.2.3 Pade Types and Partial Pades: a model example 

In real life, the number of available terms from the chiral expansion for the 
construction of a PA is very limited. Since the masses and decay constants of 
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P^AQ') 


ZR 


C4(xl0-3) 


C6(xl0-3) 


C_4(xlO-3) 


/.(XIO'3) 


P^iQ') 

p^m 


0.66 + i0.17 
0.56 + i0.28 


input 
-41.26 (6%) 


79.58 (3%) 
64.40 (21%) 


-3.379 (30%) 
-2.898 (12%) 


5.04 (5.3%) 
input 


Model Eq.(|4jU) 


Mp,mA 


-43.92 


81.81 


-2.593 


4.78 



Table 4.2: Summary of predictions for P2 (Q^) ^^id P'iiQ'^) compared with 
the real values written in Eqs. (|4.im4.lB|) . In parenthesis we show the 
relative error defined as P'-edictedvaiue-reaivaiue Qf g^^^ prediction. 

the first few vector and axial-vector resonances are known, one may envisage 
the construction of a rational approximant having some of its poles at the 
prescribed values given by the known masses of these resonances. If all the 
poles in the approximant are prescribed this way (as in the MHA), we have 
a PTA. On the contrary, when some of the poles are prescribed but some 
are also left free, then we have a PPA (see section r3.1.4|) . 

Pade Type Approximants (PTAs) 

Assuming that the first masses are known, let us proceed to constructing 
the PTAs, Eq. ([3:32]) . The lowest such PTA is T§(Q2)^ which contains 
two poles at the physical masses of the p and the first A in the tower, the 
lowest one in each tower, Eq. (|4.13|) . Fixing the residue through the chiral 
expansion to be Co 



-Fq , one obtains 



tO(Q2) 



F^APpMl 



(Q2 + M2)(Q2 + Ml) 



(4.20) 



Even though it has the same number of inputs (Co and the two masses), this 
rational approximant does not do such a good job as the PAs Eq. (|4.17p 
or Eq. (|4.19p . For instance, C_4 is 2.3 times larger than the true value in 
Table WA\ As we have already stated, one way to intuitively understand this 
result is the following. The OPE is an expansion at Q^ — t- 00 and therefore 
knows about the whole spectrum because no resonance is heavy enough with 
respect to Q^ to become negligible in the expansion, i.e., the infinite tower 
of resonances does not decouple in the OPE. Chopping an infinite set of 
poles down to a finite set may be a good approximation, but only at the 
expense of some changes. These changes amount to the appearance of poles 
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and residues in the PA which the original function does not have. This is 
how the PA in Eq. ()4.17p manages to approximate the true function Eq. 
()4.12p . However, by construction, the PTA Eq. (|4.20p does not ahow the 
presence of any artificial pole because, unlike in a PA, all its poles are fixed 
at the physical values. Consequently, it only has its residues as a means to 
compensate for the infinite tower of poles present in the true function and, 
hence, does a poorer job than the PA Eq. (|4.17p . particularly in determining 
a large-Q^ observable like C_4. Indeed, the role played by the residues in 
the approximation can be appreciated by comparing the true values of the 
decay constants to those extracted from Eq. ()4.20p . Although the one of 
the p is within 30% of the true value, that of the A is off by 100%. 

A different matter is the prediction of low-energy observables such as, 
e.g., the chiral coefficients, C2k for /c > 0. In this case heavy resonances make 
a small contribution and this means that the infinite tower of resonances 
does decouplelj Truncating the infinite tower down to a finite set of poles is 
not such a severe simplification in this case, which helps understand why a 
PTA may do a good job predicting unknown chiral coefficients. Indeed, Eq. 
(j4.20p reproduces the value of C2 within an accuracy of 15%, growing to 
22% in the case of C4. A global observable like I-,^ averages the low and the 
high Q^ behaviors and ends up differing from the true value Eq. ()4.18p by 
35%. This gives some confidence that observables which are integrals over 
Euclidean momentum may be reasonably estimated with MHA as, e.g., in 
the Bk calculation of Ref. [HI fTT5lfTT6] . 

Improving on the PTA of Eq. (|4.20p by adding the first resonance mass 
from the vector tower produces the following approximant 

^pw^2^_ a + bQ^ -tV, / "= -13.5x10-3, 

^3W } (Q2 + M2)(Q2 + M2)(Q2 + M2 ) ' ^'^"^ \ b = +1.33 X 10-4 , 

(4.21) 
where the values of the chiral coefficients Co and C2 have been used to 
determine the parameters a and b. The prediction for C4 is much better now 
(only 2% off), in agreement with our previous comments. The prediction 
for C_4 is still very bad, becoming now 19 times smaller than the exact 
value. Table 133] shows, for the lowest PTAs, the predictions for the Taylor 
coefficients C2,C4,Cq, the first OPE coefficient C_4, the prediction for the 
integral of Eq. (j4.18p and also the prediction for the lowest residue, i.e., Fp, 



This is because the residues F in the Green's function Eq. (|4.12|l stay constant as 
the masses grow. This behavior does not hold in the case of the scalar and pseudoscalar 
two-point functions [112) . 
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KiQ^) 


C2 X 103 


d X 10^ 


Cq X 10^ 


C_4 X 10^ 


U X 10^ 


Fp X 10^ 


T0(Q2) 

Ti(g2) 


17.8 (15%) 
input 


-34.0 (22%) 
-43.2 (2%) 


60.5 (26%) 
79.5 (3%) 


-6.0 (133%) 
-0.1 (95%) 


6.5 (35%) 
4.2 (12%) 


112.8 (16%) 
131.0 (2%) 


Eq. dlllD 


21.01 


-43.92 


81.81 


-2.593 


4.78 


133.8 



Table 4.3: Summary of predictions for T2((3^) and T\{Q'^) compared with 
the real values of Eqs. (|4.im4.lB|) . In parenthesis we show the relative error 
of each prediction. 



Eq. ()4.13p . Nevertheless, it eventually gets much better if PTAs of very high 
order are constructed. For example, as shown in the Table [331 not only the 
OPE coefficients are improved hierarchically by increasing the order of the 
PTAs, but also the Taylor coefficients and the integral of Eq. (j4.18p improve. 
For instance, we have found C_4 = —2.58 x 10"^ for the approximant Tg 
with 9 poles. 

However, another matter is the prediction of the residues. For instance, 
the prediction for the decay constant of the state with mass My in Eq. (|4.2ip 
is smaller than the exact value in the model in Eq. ()4.13p by a factor of 2. In 
general, we have seen that the residues of the poles always deteriorate very 
quickly so that the residue corresponding to the pole which is at the greatest 
distance from the origin is nowhere near the exact value. We again explicitly 
checked this up to the approximant Tg, in which case the decay constant 
for this pole is almost 5 times smaller than the exact value. The conclusion, 
therefore, is that PTAs are able to approximate the exact function only at 
the expense of changing the residues of the poles from their physical values 
(this fact can be check in the last column of Table 14. 4p . Identifying residues 
with physical decay constants may be completely wrong in a PTA for the 
poles which are farthest away from the origin. 



Partial Fade Approximants (PPAs) 

As an intermediate approach between PAs and PTAs, there are the PPAs, 
Eq. p.29p . where some poles are fixed at their physical values while some 
others are left free. From a theoretical point of view one should expect 
PPAs to have both advantages from PAs and PTs. The simplest of such 
rational approximants is Pj i(Q^) (see section [3.1.4l for notation). Fixing its 
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n^Q') 


C_4 X 103 


U X 10^ 


Fp X 10^ 


Fa X 10^ 


Flast X 103 


Ti(Q2) 
Tf(Q2) 


-1.96 (24%) 
-2.48 (4.2%) 
-2.58 (0.6%) 


4.73 (1.2%) 
4.782 (0.12%) 
4.787 (0.02%) 


133.87 (0.01%) 
133.884 (??%) 
133.884 (??%) 


142.6 (0.8%) 
143.7 (0.01%) 
143.76 (??%) 


51.6 (64%) 

39.7 (72%) 
30.1 (80%) 


Ea.()4.11|) 


-2.593 


4.787 


133.884 


143.758 


143.758 



Table 4.4: Summary of predictions for tl{Q'^), Tf(Q2) and T^(Q2) com- 
pared with the real values of Eqs. (|4. 11114. 13]) . We have included here the 
prediction of the second residue which corresponds to ruA ■ The last column 
is the prediction of the last residue of each PT which corresponds to the 
fifth, seventh and ninth masses of the spectrum for the respective PTs. In 
parenthesis we show the relative error of each prediction. (??%) means too 
close agreement that can not be shown in the Table. 



3 unknowns with M^, Cq and C2, one obtains 



AAQ' 



_ J. 



R 



(Q2 + M2)(Q2 + ^^) 



, with r| = 3.75x10^3 



ZR = 0.8665 . 

(4.22) 

As can be seen, the mass (squared) of the first A resonance is predicted to 
be at zr which is sensibly smaller than the true value in Eq. (|4.13p . 

At that point a comment is in order because, intriguingly enough, this 
is also what happens in the real case of QCD |4H W2\ I117J . The Partial 
Pade P5 i(Q^) predicts the location of the uia at 0.931GeV when its value 
in the model is 1.182GeV, Eq. ()4.13p . The same kind of discrepancy was, 
actually, found in [42j. In that reference, RChT (see sec l2.5p is used to 
predict the value of the ruA when saturating the spectrum of the theory 
to one resonance per channel (only Mp and mA) and fixing the value of 
Mp to the physical value. Technically, this is a Partial Pade since ruA is 
left free and low-energy properties are used to fix the unknowns. Then, 
to fix its value, the WSR were used and a particular relation between Mp 



and mA was found, i.e, ttt-a = MpJ\ — f^/Fy = 0.968GeV. At that time 
they used /,r = 93.3MeV and Fy = 154MeV. Nowadays, with the values [5] 
we will find niA = 0.962GeV which has to be compared with the physical 
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rriA = 1.230 it 0.040GeV. The authors agreed that the value predicted was 
a good determination of the physical rriA. 

A similar example is found in Ref. jllTj where the MHA is used and 
the value of Mp is again fixed to the physical value [5]. After imposing 
several OPE and low-energy constrains, the value of m^ is obtaineco, niA = 
(938.7ibl.4)MeV. Since the set of equations used by these authors is different 
that the set used in Ref. |42j . the final value for rriA is also different, but still 
smaller than the physical one. 

Another similar case can we found in Ref. |118j . They used also a Partial 
Pade where the Mp was fixed to the physical value but to predict tua they 
used Fa coming from the radiative decay oi — )• tt^ Eq. ()4.5ip . After imposing 
several OPE constraints, they found niA = 998 it 49MeV which is, again, 
smaller than the physical value. 

Since all three determinations are done in the framework of large- A^c 
QCD, one could interpret the results as a large- A^c ruA instead of a real niA 
[118) . From the Pade Theory point of view, however, one could understand 
this discrepancy in another way without needing to say that this prediction 
is a large- A'c effect. As can we seen in Table 14.51 the Euclidean quantities 
like C4 or Cq are rather insensible to the value of niA used, even if it is a 
complex number. Then, in the other way, imposing C4 or Cg could lead to 
a different values of niA depending on the set of equations used. We can not 
exclude the conclusion that these references addressed but it may also be 
that the location of the second pole zr in Eq. ()4.22p has nothing to do with 
the physical value of tua ■ 

The rational function Eq. (fi:22|) predicts C_4 = -r| = -3.75 x 10~^ 
which is a better determination than that of the PTA Eq. (j4.20p with the 
same number of inputs, and C4 = —45.52 x 10~^ which is not bad either. 
Concerning the pion mass difference, one gets /•„- = 5.22 x 10~^. However, 
as compared to the PAs Eq. (|iT7l) or Eq. (jil^ . the PPA Eq. ()i:22|l does 
not represent a clear improvement. 

In order to improve on accuracy of the PPA, we can try to go to the 
Partial Pade P2i(<5^) or its counterpart F\2{Q'^)- The first one needs 
Mp,Co,C2 and C4, and the second one needs Mp,mA,Co and C2. The 
results are shown in Table H31 Since the predictions for the decay constant 
and for the pion mass difference for both cases are very accurate, one may 



^'^In this reference the authors also studied what would happen if one includes one extra 
vector state in the approximation, a MHA+V. They found that the impact of including 
higher states in the determination of tua is irrelevant. We will also discuss this result in 
sec. 11231 
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try to use the mass and decay constant of the first resonance, M^ and Fp, in 
addition to 1,^ and the chiral coefficients Co, C2 and construcl|ij the P2 liQ'^): 
which can be written as: 



F^M^ a - F'JMJ Q2 (^ ^743 ^ ^^q- 



3 



2,1^^ ) Q2 + M2 ^ (Q2 + ^^) (Q2 + ^*) ' I ^^ = 1.24 + i 0.34 . 

(4.23) 
This PPA, upon expansion at large and small Q^, determines C_4 = — 2.47x 
10~^ and C4 = —44.0 x 10^'^ to be compared with the corresponding coeffi- 
cient in Table 14. li The accuracy obtained is better than that of Eq. (|4.19|) , 
but this is probably to be expected since Eq. (|4.23p has more inputs. 

Based on the previous numerical experiments done on the model in Eqs. 
(|4.12l I4.13P (and many others), we now summarize the following conclu- 
sions. Although, in principle, the PAs have the advantage of reaching the 
best precision by carefully adjusting the polynomial in the denominator to 
have some effective poles which simulate the infinite tower present in Eq. 
(j4.12p . they have the disadvantage that some of the terms in the low-Q^ 
expansion are required precisely to construct this denominator. This ham- 
pers the construction of high-order PAs and consequently limits the possible 
accuracy. 

When the locations of the first poles in the true function are known, 
there is the possibility to construct PTAs (with all the poles fixed at the 
true values) and PPA (with some of the poles fixed and some left free). 
As we have seen, although the PTA may approximate low-Q^ properties of 
the true function reasonably well, the large-Q^ properties tend to be much 
worse, at least as long as they are not of unrealistically high order. The 
PPAs, on the other hand, interpolate smoothly between the PAs (only free 
poles) and the PTAs (no free poles). Depending on the case, one may choose 
one or several of these rational approximants. However, common to all the 
rational approximants constructed is the fact that the residues and/or poles 
which are farthest away from the origin are in general unrelated to their 
physical counterparts. 



4.2.4 Comment on Complex poles and their consequences 

Artificial poles and analytic defects are transient in nature, i.e. they 
appear and disappear from a point in the complex plane when the order of 



^^This is not, strictly speaking, a PPA since we use here /^ and Fp instead of Taylor 
coefficients. We use the tilde to identify this approximant. Although we do not lie on a 
strict mathematical theorem the results are quite good, see Table B31 
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niA 


C4 X 10^ 


Ce X 10^ 


C_4 X 10^ 


U X 10^ 


Fp X 10^ 


p?,l(Q') 


0.867 


45.5(3.6%) 


88.7(8.4%) 


-3.75(45%) 


5.22(9%) 


151.4(28%) 


pil(Q') 


1.3 + i0.3 


input 


input 


-2.32(10%) 


4.74(0.9%) 


133.6(0.2%) 


P},2(Q') 


input 


input 


81.9(0.2%) 


-2.03(22%) 


4.68(0.5%) 


134.4(0.4%) 


nM) 


1.2 + i0.3 


44.0(0.2%) 


82.1(0.3%) 


-2.47(5%) 


input 


input 


Eq.(|4.11D 


1.182 


43.9 


81.81 


-2.59 


4.78 


133.88 



Table 4.5: Summary of predictions for P?i(Q2)^ IP2,i(<5^)' ^i^iiQ"^) and 
¥\ 2(Q^) compared with the real values of Eqs. (|4. 11114. 13]) . All the results 
are shown in absolute value. In parenthesis we show the relative error of 
each prediction. 

the Pade approximant is changed, that is the typical pattern for a defect or a 
Proissart doublet. Sometimes their presence mean certain noise in the Taylor 
expansion used the built them, sec. 13.1.31 On the contrary, the typical sign 
that a pole in a Pade is associated with a truly physical pole is its stability 
under these changes in the order of the Pade. Of course, when the order in 
the Pade increases there have to be new poles by definition, and it is natural 
to expect that some of them will be defects. Pade Approximants place some 
effective poles and residues in the complex Q^ plane in order to mimic the 
behavior of the true Green's function, but it can mimic the function only 
away from the poles, e.g. in the Euclidean region. Obviously, PAs cannot 
converge at the poles, in agreement with Pommerenke's theorem [89], since 
not even the true function is well defined there. The point is that what may 
look like a small correction in the Euclidean region may turn out to be a 
large number in the Minknowski region. 

To exemplify this in simple terms, let us consider a very small parameter 
e and imagine that a given Pade P(Q^) produces the rational approximant 
to the true Green's function G{Q'^) given by 



G(Q2) ~ P(Q2) = ^(g2) ^ 



g2 + M2 ' 



(4.24) 



where R{Q^) is the part of the Pade which is independent of e. Although 
for (5^ > there is a sense in which the last term is a small correction 
precisely because of the smallness of e, for Q^ < Q this is no longer true 
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because of the pole at Q^ = — M^. This pole is in general a defect and may 
not represent any physical mass. In fact, associated with this pole, there is 
a very close-by zero of the Fade P{Q'^) at Q^ = -M^ - e i?(-M2)"\ as 
can be immediately checked in (j4.24p . This is another way of saying that 
a defect is characterized by having an abnormally small residue and is the 
origin of the pairs of zeros and poles in the y-shaped branches of Fig. 14.21 
Therefore, not only are defects unavoidable but one could say they are even 
necessary for a Fade Approximant to approximate a meromorphic function 
with an infinite set of poles. 

Similarly to masses, also decay constants may be unreliable. To see this, 
imagine now that our Fade is given by 

again for a very small e. As before, the term proportional to e may be 
considered a small correction for Q^ > 0. However, at the pole Q^ = — M^ 
the decay constant becomes F + e"^ which, for e small, may represent a huge 
correction. When the poles are preassigned at the physical masses, like in 
the case of FTAs, it is the value of the residues that compensates for the 
fact that the rational approximant lacks the infinite tower of resonances. As 
we saw before, the residues of the poles in the Fade which lie farthest away 
from the origin are the ones which get the largest distortion relative to their 
physical counterparts. 

These comments could be related with the results found in Ref. |117] . In 
that reference the properties of the function nifl(Q^) in the Large- A'^c limit 
are also studied assuming the MHA. MHA has a spectrum consisting of a 
pion state, a vector state and an axial vector state. Fart of the motivation of 
that study was to understand a discrepancies between various phenomeno- 
logical analyses Refs. [m [Ml 1201 [Ml [Ml [Ml [M] using the same 
experimental data on hadronic r-decays published by the ALEFH |122j and 
OFAL [123j collaborations at LEF. Table 1 of Ref. [11 7j shows the conflicting 
results for the chiral condensates Eq. ()4.5p which modulate the asymptotic 
behavior of the n/,ij(Q^) function at large Q^ valuea^n. To perform their 
analysis, the authors of Ref. |117j study two successive approximations to 
nLi?(Q^) in the large-A^c limit. After imposing some good high-energy be- 
havior in several Green's functions and form factors throug a fit procedure, 
Ref. [117] obtains, keeping only one vector state V and one axial- vector state 



^■^With our analysis we can also provide some insights in this discussion but we left it 
to the end of sec. 14.3.11 
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A, the precise values My = 775.9 ± 0.5 MeV and Ma = 938.7 ± 1.4 MeV. 
Actually, they performed a Partial Pade P^ ^ since they fixed the My to 
the physical value and led Ma free. Going one step forward, they improved 
the approximation adding an extra higher vector state, a MHA+y. They 
again performed a Partial Pade, but now a P} 2 since they fixed the My to 
the physical value and led Ma and Myi to be free. However, in this last 
case, imposing My = (775.9 ± 0.5) MeV, they found Ma = (939.4 ± 1.1) 
MeV and M[, = (1258.2 ± 2.5) MeV. The value of Ma is rather stable. 
More interesting, they found a zero in the numerator of the approximation 
at 1262.9 lb 2.5 MeV. Pole and residue are so close that they are almost 
canceled. Therefore, when they expand at large values of Q^, for both MHA 
and MHA+y they found very similar results showing that the inclusion 
of a higher vector state does not improve the approximation. They argued 
that this situation is consistent with the fact that the MHA approximation 
seems to have already the hulk of the full large-Nc information. In other 
words, adding an extra V'-pole appears to be compensated, at a very good 
approximation, by the position of the nearby zero. 

Prom our point of view, we suspect that this compensation could be 
understood as a Froissart doublet instead of a physical pole, and then the 
value My = (1258.2 it 2.5) would not correspond to the physical one. We 
have also explored a Partial Pade sequence. Table 14.51 and we have seen 
that letting the poles be complex helps improving on the predictions for the 
OPE and Taylor coefficients. 



4.3 Selecting the best approximant sequence 

At the beginning of this chapter, at sec. 14.11 we argued that we would use 
PA of the form P^f_^_2 because the fall off the function Q'^HiniQ'^) has a fall 
off of Q"^. We also said that this procedure is not a theorem and could be 
a source of a systematic error on our predictions. Now we are in position to 
discuss the quality of this choice with the help of several examples. 

In our case we have the advantage of certain knowledge of the behavior 
at Q^ — )• 00 that the function has. However, from a more realistic point of 
view, we may know just few coefficients of its Taylor expansion and nothing 
more. Neither poles, nor residues, and even nor OPE coefficients may be 
known. In that situation, a priori, one can not decide what is the best 
kind of PA sequence to approximate the function. In our model Eq. ()4.1ip . 
since we know a priori certain global information about the function, we 
can impose extra constraints to our PAs. But now, we would like to test 
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whether assuming -Pjv/+2 ^^ correct or not. For that reason, we will perforin 
an exercise. Lets suppose that only the first Taylor coefficients of our model 
Eq. ()4.1ip are known, i.e. only Cq + C2Q^ + • • • + C2q{Q'^Y^ are known from 
Eq. (j4.16p . With all of these coefficients, we can construct several PAs, such 
as Pf,Pi, Pi, P!,Pi, Pi, Pf, Pi, Pg- Then, with ah of them we can try to 
predict the position of the first pole and residue of the original function 
and also the first not used Taylor coefficient, namely the C22- Defining the 
relative error e as e = f^^'^^rfedm^t.e-rea; m;«e ^ ^-^^ ^^^^^^^ ^^^ ^-^^ prediction 

of the first pole position, the first residue and the C22 are shown in Fig. 14.31 

In all these pole, residue and C22 predictions the best estimation comes 
from the same Pg . This PA has an extra particularity, is the only PA from 
that group with the same fall off behavior as Q~'^ as our model Eq. ()4.1ip 
has, and then, we can also estimate the value of the first nonvanishing OPE 
coefficient. In that case, we obtain a prediction for C-4 with a relative error 
£ = 4-10-1 

Also, as already shown in section 14.21 the analytical structure of the Pg 
follows the same pattern as its counterparts shown in Fig. 14. 1[ As can be 
seen in Fig l4.4l the lowest poles are close to the original once while when 
one goes away from the origin, its location worsens becoming eventually 
complex. 

PAs are not only useful to predict Taylor coefficients and the lowest 
masses and residues but also tell us global properties of the original function 
without the need to impose any extra constrain. In such a way, we obtain 
not only a prediction for the first OPE coefficient but we predict the right 
fall off behavior of our model. 

To show that this characteristic behavior of the PAs is a general fact and 
not model dependent, we have also studied a couple of different models with 
different properties. However, in all of them we find the same conclusion. 
The first of theses models is close to which we have already studied, that 
means, is also a meromorphic function with a particular falloff when Q^ — )• 
00. Lets define the function G{Q'^) as: 



7?2 I?2 7712 00 ^ p2 Tp2 

^\^ ) ci^^ rft ^ A/r2 n2_L A/r2 "^Z^ 



Q2 ' Q2 + ^2 Q2 + ^J2^ ' z_^ I Q2 + Af2 (n) Q2 + M\{n 

(4.26) 
Again, we impose the dependence on the resonance excitation number n 
as in Eq. (H^]): 

p2^^(n) = p2 = constant , M^_^(n) = m\j^ + n A^ . (4.27) 
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Figure 4.3: Determination of the position of the first pole (top), the firs 
residue (center) and the Taylor coefficient C22 (bottom) from several PAs. 
See the text for details. 
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Figure 4.4: Location of the poles (red) and residues (green) in the complex 
g^-plane for the Pi. 



Thanks to Eq. ()4.27p . the infinite sum in Eq. (|4.26p is summable. Also, 
since we are not now concern about the physical properties of the model, 
just mathematics, we do not need to reassemble QCD. We choose a set of 
parameters satisfying this time the fall off constraint that the OPE starts 
at dimension 10: 



Fo = 85.8 MeV , Fp = 97.7 MeV , Fao = 110.8 MeV , F = 143.758 MeV(,4.28) 



Mp = 0.385 GeV, 



Mao = 1-268 GeV, my = 1.494 GeV , A = 1.277 GeV . 



For instance, Fp, Mp, Fag and Ma^ are chosen so that the function G{Q'^) 
in Eq. I^^7IU\\ has vanishing 1/Q^ 1/Q^, 1/Q^, 1/Q^ in the OPE at large 
Q^. Defining the expansion of our G{Q^) function around Q^ = 0, Q^ — t- oo 
as 



Q2 G(-Q2 



E 



921 Q 



2k 



with / = 0,±1,±2,±3, , 



(4.29) 



one obtains the values for the coefficients of both Taylor expansion 
(/ > 0) and OPE expansion {I < 0). In Table HT6l we collect few of these 
coefficients for the function G{Q'^). 

Following the definitions presented in sec. 14. H we can proceed with the 
construction of rational approximants to the function G(Q^). Again we 
will play the same game than before, i.e., assuming that we only know 
the lowest Taylor coefficient of G(Q^) and nothing more, i.e, go + 52Q^ + 
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90 


52 


54 


56 


5-8 


5-10 


5-12 


-7.362 


61.18 


-429.6 


2901.7 


-4.95 


22.46 


-61.87 



Table 4.6: Values of the coefficients §21 from the high- and low-Q^ expansions 
of Q^ G{-Q^) in Eq. (lOSD in units of lO'^ GeF^-s/. Notice that 5.2 = 
g_4 = 5-6 = and go = —Fq, see the text for details. 



■■■ + 52o(Q^)^'^) Eq- (|4.29p . With all these coefficients, we can construct 
the PAs Pf, Pi, Pi, Pi, Pi, Pi, P7,Pi, Pi- Once all these PAs are know, by 
reexpanding them we can predict g22, the first Taylor coefficient not used in 
their construction, and also we can try to predict the first pole and residue 
of G(Q^) defined in Eq. (|4.28p . Fig. 14.51 shows the results for the case of the 
model Eq. (j4.12p , using also the definition of relative error already presented. 
The preferred approximant is now the Pf, which gives the best prediction 
for the first pole and residue position, and also for (722 • Also, is the only PA 
among them that has the right fall off as Q~^ as the function G{Q'^) has, 
giving an extra prediction for the first nonvanishing OPE coefficient with a 
relative error eg_g = 1%. 

Again, its poles and residues. Fig. 14.61 are located in the complex Q^ 
plane in a similar way than the cases shown in Figs. 14.11 14.21 and 14.41 

At that point, one can think that this selection-rule-behavior shown by 
PAs comes from the fact that the function we are approximating is mero- 
morphic. To this end, we would like to consider also an example of a Stieltjes 
function and study again which PAs is the best approximant when only few 
of the Taylor coefficients are known and no more information about the 
function is available. Let us consider the following function defined by the 

sum [innlfTnT] 



H{Q') = y: 



1 



(4.30) 



^ n{Q'^ + n) 

as a very simple toy model for the once-subtracted vacuum polarization 
function in the large- A'^c limit, |93j . Obviously its spectrum consists of masses 
located at M^ = 1, 2, 3, ..., and one may use the mass of the lowest resonance 
as the unit of energy. An exact result for H[Q^) can be readily obtained by 
summing Eq. g3Q]), [93]: 



HiQ') = ^{HQ' 



l) + 7} 



(4.31) 
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Figure 4.5: Determination of the position of the first pole (top), the firs 
residue (center) and the Taylor coefficient 322 (bottom) from several PAs, 
approximating the function G{Q'^). See the text for details. 
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Figure 4.6: Location of the poles (red) and residues (green) in the complex 
q^-plane for the Pj. 



where the function il^{Q'^ + 1) is the Digamma function, defined as 

^P{z) = —\ogT{z) , 



(4.32) 



with T{z) the Euler Gamma function and 7 ~ 0.577216 being the Euler- 
Mascheroni constant. The function H{Q'^) has the low-energy expansion 



p=0 p=0 



(4.33) 



where CiP + 2) = Xlnii ^ ^ ^ is the Riemann's C function. On the other 
hand, the large Q^ expansion is given by (|Arg((5^)| < it) 



HiQ') 



7 + logQ^ 



+ 



E 



B 



2p 



2Q4 Z^ 2p Q4P+2 

p=l 



(4.34) 



(_l)p+i ( P^^ C(2p) are the Bernoulli numbers for p 



where i?2p = (— 1)p+-^ ('2npp C(2p) are the Bernoulli numbers for p = 1, 2, 3, .... 
While the chiral expansion has a finite radius of convergence \Q'^\ < 1, 
the large Q^ expansion has zero radius of convergence (i.e. it is asymptotic) 
due to the factorial growth of the Bernoulli numbers. Notice also the pres- 
ence of the logQ^ behavior in ()4.34p . That is the main different between 
this example and the previous one. We can not think about a certain falloff 
when Q'^ — )• 00 of the original function since H{Q^) diverges. However we 
can again do the same game with PAs. We just take the first Taylor coeffi- 
cients of Eq. (|4.33p . ho + /i2<3^ + • • • + h2o{Q^y^, and construct the lowest 
PAs Pf,Pi,P^,Pi,Pi,Pi,P^,Pi,P^ . We can, then, predict the position 
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of the first pole and residue, and by reexpanding around (^ = 0, predict 
/i22- The results are shown in Fig. 14.81 The best approximant in this case is 
the diagonal P^ which also manages to be de one with the closest behavior 
as a plateau as the ratio Q~^ \og{Q~^) has. Of course, we will be never able 
to predict the very first coefficient of the OPE expansion since this is not 
a coefficient but a logarithm. However, PAs somehow organizes the infor- 
mation coming from the Taylor expansion to describe properly the function 
even for a high values of Q^. That can be easily tested numerically since is 
a direct application of the convergence theorem for Stieltjes functions, i.e., 
limQ2^Q.2 \H{Q^) - P^\Q^)\ < e, for any large Q*^ |Arg(Q*2)| > ^^ ^nd 
small e. 

Trying to better understand the particular case of model Eq. (|4.34p . 
we have also studied a different set of PAs, now restricting ourselves to 
one less Taylor coefficient, /120, giving the opportunity to our set to have 
a subdiagonal PA since the previous exercise had not. We can construct, 
then, the set Pf, P^-, P^,PI,P^,Pq,Pj, Pg and again predict the first pole, 
the first residue and the first predictable Taylor coefficient. The results 
are shown in Table HTTl In that last case, the best approximant is the P^, 
which predicts the mean piece of the first coefficient in the OPE, Eq. ()4.34p . 
Figs. 14.71 14.81 tell us that best PA lies half way as a subdiagonal PA, half 
way a diagonal PA, and if one can construct both PA sequences, one will 
see that the original function can be constrain with an upper and a lower 
limit, as pointed out by the convergence theorem for Stieltjes functions, Eq. 
iKMi in sec. ED 

An interesting feature shared by both PAs is that its analytical structure. 
Since both are approximating a Stieltjes function, their poles and residues 
must be all reals, according again to the convergence theorems of sec. 13.11 
That is shown in Fig. 14.91 At the top, the lowest poles (red points) and 
residues (green points) of P5 and at the bottom, the lowest poles (brown 
points) and residues (sky-blue points) of P5 . All of their poles are actually 
purely reals. 

The convergence theorem for Stieltjes functions guarantees convergence 
for all the Pade -P^"*" {J > —1) sequences. As we have said, for J > — 1 
all the poles of our PAs have to lie on the same real axes where the brunch 
cut is located. Now we have the opportunity to see what would happen if 
we use J < —1. If we look to the analytical properties of PQ,Pj,Pg and 
pf, P^,Pg we see that they have several complex conjugate poles. Actually, 
there is also a way to count the number of complex conjugate poles given a 
certain value for J. If a certain PA has degree Pj^, and M = N + k, the 
PAs will have, at least, A^ + 1 real poles and {k — l)/2 complex conjugate 
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pares when possible, i.e., for odd values of k. 



4.3.1 Beyond Pade Approximants from Taylor Expansion 

Throughout the last section we have carefully studied Pade Approximants 
constructed from a Taylor expansion. Now, two questions may be addressed 
in order to compleat our analysis. The first one is what would be the role 
of other Pade Approximants, such as Pade-Type Approximants, compared 
with the Pade Approximants? The second is what may happen if instead 
of the Taylor expansion defined when Q^ — )• we use the expansion when 
Q^ — >• oo to construct Pade Approximants? In the present sub-section we 
will try to answer both question. 



Pade Approximants vs Fade Types 

The first answer would be better understood if one compares PTs with 
PAs defined to be of the same degree, i.e, P^^ and T^ . Using again the 
model of Eq. ()4.10p and using only the first 11 coefficients of its Taylor 
expansion, Eq. (j4.7p . i.e., Cq + C2Q^ + • • • + C2q{Q'^Y^ ■, we construct the 
following set of PAs: Pf,Pi,P^,Pi,Pi,Pi,P^,Pi,P^. Now imagine that 
we construct a similar set of PTs with the same Taylor expansion. Now, all 
these PTs have their denominator fixed in advance. This will lead to the 
following set of PTs: T?, T|, T^, T|, T|, T^, Tf, T^, T^. We have now exactly 
the same number of input parameters, 11, but as a mixture within Taylor 
coefficients and pols of the function. For example, T^ has only one pole fixed 
in advance, the M^ and 10 Taylor coefficients. On the contrary, Tg has 9 
correctly located poles and only one Taylor coefficient. When all the PTs 
are constructed, by reexpanding we can predict the Taylor coefficient C22 
and compare with the results found in Fig. 14. 3i The comparison is shown 
in Fig. 14.101 where the red points are the predictions for the PTs and the 
orange points are the predictions for the PAs. 

From the point of view of PAs the best approximant is the Pg which also 
gives the best prediction for the first pole and residue. It is also the only 
PA from the selected set able to give a prediction for the first nonvanishing 
OPE coefficient. However, from the PTs point of view things are different 
since the best PT is not the Tg as one would expect but T4 which has a 
divergent behavior for large Q^. This PT also gives the best prediction for 
the first residue (remember that all the poles are fixed in advance), while 
is hard to see difference between T4 and T3 or T5. In any case, Tg differs 
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Figure 4.7: Determination of the position of the first pole (top), the firs 
residue (center) and the Taylor coefficient /i2o (bottom) from several PAs, 
approximating the function H{Q'^). See the text for details. 
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Figure 4.8: Determination of the position of the first pole (top), the firs 
residue (center) and the Taylor coefficient /i22 (bottom) from several PAs, 
approximating the function H{Q'^). See the text for details. 
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Figure 4.9: Analytical structure of the P^ (top) and ^5^ (bottom), both 
approaching HiQ"^), Eq. (gJO]). 
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point and Pj^ when is an orange point. 

numericahy one order of magnitude from T| in the prediction of C22 or for 
the first residue. 

With the same among of input, PAs play a better role on determining 
Taylor coefficients and residues that their counterparts the PTs. That hap- 
pens because PAs have an extra degree of freedom that PTs do not have. 
They can have complex conjugate poles to mimic the properties of the model 
when PTs are constrained to have real poles. The advantage of the first are 
clear since the difference among them are of several order of magnitudes. 
Of course, concerning masses PTs are un a different league since masses are 
inputs in that case. 

This comment agrees with what was found in the Ref . p5J . In that refer- 
ence, the authors show that many short-distance constraints (which means 
imposing to the approximant the right fall off when Q^ — )• 00) can be easily 
incorporated in a general model based on resonance saturation. At the same 
time, however, they pointed out that these approximations cannot reconcile 
all short-distance constraints due to a general conflict between those con- 
strains coming from Green functions, and those coming from form factors 
and cross-sections. The authors also show that if one wants to remedied this 
situation without spoiling the matching into the OPE, one can not proceed 
using a single or any finite number of resonances per channel type of ap- 
proximations. In that sense, an earlier example where single resonance does 
not allow to reproduce all short-distance constraints can be found in Ref. 

ESI. 
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Using the OPE to construct Pade Approximants 

Concerning the second question we addressed ourselves, in 1977, A. A. 
Migdal [124j suggested PAs as a method to extract the spectrum of large- Ac 
QCD from the leading term in the OPE of the (VV) correlator, i.e., from 
the parton model logarithm. However, nowadays this proposal should be 
considered unsatisfactory for a number of reasons [1251 fT26j , the most simple 
of them being that different spectra may lead to the same parton model 
logarithm [1271 1128) . In fact, the full OPE series is expected to be only an 
asymptotic expansion at Q^ — t- oo (i.e. with zero radius of convergence) , and 
PAs constructed from this type of expansions cannot in general reproduce 
the position of the physical poles |129j . Migdal's approach has been recently 
adopted in some models exploiting the so-called AdS/QCD correspondence 
[130) and, consequently, the same criticism also applies to them. And now 
we can address this criticism based on Pade Theory just constructing PAs 
from the 1/Q^ expansion at infinity. 

We will show how the PAs constructed from this expansion (akin to the 
OPE) do not in general reproduce even the first resonances in the spectrum 
of the model Eq. ()4.12p . unlike those constructed from the chiral expansion. 
Recalling the definition of the OPE for the model Eq. ()4.12p . given in Eq. 
()4.5p . with the corresponding coefficients, Eq. ()4.15p . it is straightforward 
to construct a PA in l/Q^ around infinity, i.e. by matching powers of the 
OPE in l/Q^. Since the function Q^Ii^LRi—Q^) behaves like a constant for 
Q^ — )• 0, we will consider diagonal Pade Approximants, i.e., of the form 
P^{l/Q'^), in order to reproduce this behavior. Figure 1^.111 shows the po- 
sition of the poles and zeros of the PA P|q (— 1/g^) in the complex q'^ plane 
(recall that Q^ = —q^)- As it is clear from this plot, the positions of the 
poles have nothing to do with the physical masses in the model, given by 
Eqs. ()4.9M.13p . even for the lightest states. This is to be contrasted with 
what happens with the PA constructed from the chiral expansion around 
Q^, which is shown in Fig. 14. 2i The difference between the two behav- 
iors is due to the fact that, while the chiral expansion has a finite radius 
of convergence, the radius of convergence of the OPE vanishes because this 
expansion is asymptotic. In spite of this, one can see that the PA P^q (1/Q^) 
is an excellent approximation to the function Q'^^lr{—Q'^) in the euclidean 
region. 

On interesting feature of that expansion is that one can try to predict 
chiral coefficients. 

The simplest PA with only one resonance per channel constructed using 
the OPE expansion is the PA: 
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Figure 4.11: Location of the poles (red-dots) and zeros (blue-squares) of the 
Fade Approximant P^^{-l/q2), constructed from the OFF of Q'^Ii-LRi-Q^) 
defined in Eqs. ()4.16l and I4.15p . in the complex q'^ plane. We recall that 
Q^ = —q^. The poles are all complex-conjugate pairs. 



piiVQ' 



-F^ + 






(4.35) 



For the four unknown parameters, we can use the first four OFF coeffi- 
cients, i.e, the following set of equations |131j : 

(4.36) 





7?2 


-Fl = 


7?2 
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-Flm\ 
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Flm\ - 
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F^ml - Flm\ 



(4.37) 



where eg and eg are the OFF coefficients of order 0{Q^^) and 0{Q~^) 
respectively, Eq. (j4.5p . Once the four undetermined parameters of Eq. ()4.3.ip 
are known, Fq. (|4.36p . by reexpanding at Q^ — )■ one can estimate the LECs 
Fq and Lio, one obtain this nice relation among OFF coefficients and LEG 
coefficients: 



L 



10 



4e 



(4.38) 



which determines uniquely the sign of the condensate of order 0{Q ^). 
One can go improving this result by constructing higher and higher FAs 
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Figure 4.12: Predictions for Fq (left) and Lio (right) from PAs constructed 
using the OPE expansion. 



from the OPE. Fig. 14.121 shows the prediction for Fq and Lio for all the 
PAs up to the P^q{1/Q^). The sign of Lio never change, the prediction 
improves by increasing the order of the PA, and the location of the poles 
and residues are completely far away from their physical counterparts. Fig. 
14.111 In sec. 14.2.41 we have commented on Ref. [117j and the discrepancy 
on the values of the condensates of order 0{Q~^) and 0{Q~^) found in the 
literature. After our analysis, we can give a new insight on that discussion 
since we found a relation between the Lio and the OPE coefficient of order 
0(Q~^), eg. Nevertheless, we can not say anything about the relative sign 
of 0{Q~^) versus 0{Q~^) since in our result eg always appears squared. In 
any case, since we have seen that Lio has the same sign as eg, that fact may 
contradict the predictions for the OPE coefficient of order 0{Q^^) obtained 
in Refs. |lU6t \10b\ I1U8| which found different sign for it. 



89 



4.4 The QCD case 

After all the mathematical discussion about Pade Theory applied to the 
model Eq. ()4.12p . let us now discuss the real case of large- A^c QCD in the 
chiral limit. In contrast to the case of the previous models, any analysis now 
is limited by two obvious facts. First, any input value will have an error 
(from experiment and because of the chiral and large- A^c limits), and this 
error will propagate through the rational approximant. And second, it is not 
possible to go to high orders in the construction of rational approximants 
due to the rather sparse set of input data. In spite of these difficulties one 
may feel encouraged by the phenomenological fact that resonance saturation 
approximates meson physics rather well and the quite well understanding of 
Pade Theory. 

The simplest PA to the function Q^Hlr{—Q^) with the right fall-off as 
Q-4 at large Q^ is P^iQ^): 

By simple re-expansion around Q^ = it is then possible to predict 
an estimate for the term of 0((5^) in (j4.7p . In the QCD case, the values 
of the three unknowns a, A and B may be fixed by requiring that this PA 
reproduces the correct values for Fq, Licl£| and I-,^ o given by 

Fq = 0.086 ± 0.001 GeV , 
5m^ = 4.5936 ± 0.0005 MeV =^ U = (5.480 ± 0.006) x lO^^GeV^ , 
Lio(0.5 GeV) < Lio < Lio(l.l GeV) =^ L^ = (-5.13 ± 0.6) x 10"^ . (4.40) 

The low-energy constant Liq is related to the chiral coefficient C2, in the 
notation of Eq. ()4.14p . by C2 = — 4Lio. Since Liq does not run in the large- 
Nc limit, it is not clear at what scale to evaluate Lio{n) [6O I 11331 1134^ [ST]. 
In Eq. (|4.40p we have varied ^ in the range 0.5 GeV < fi < 1.1 GeV as a 
way to estimate 1/Ac systematic effects. The central value corresponds to 
the result for Liq (Mp) found in Ref. |119j . The other results in ()4.40p are 
extracted from Refs [17[ [5] . 



^•^ Since rris decouples from F^ in the large- A^'c limit, the value of Fq is estimated in Eq. 

(|4.4Up by extracting the chiral corrections from _F,r using SU{2) x SU{2) chiral perturbation 

theory, but doubling the error as compared to Ref. 1132] . 

^■^ Recall that F is, up to a constant, the electromagnetic pion mass difference SniTr |117j 



and is defined in terms of TIlr in Eq. H4.18|) 
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Obviously, the PA ()4.39p can also be rewritten as 

in terms of two poles zy^A- In order to discuss the nature of these poles, we 
will define the dimensionless parameter C, by the combination 

C = -4Lio ^ = 2.06 ±0.25 , (4.42) 

where the values in ()4.40p above have been used in the last step. Imposing 
the constraints (j4.40p on the PA (|4.4ip one finds two types of solutions 
depending on the value of Q: for C > 2 the two poles zy^A are real, whereas 
for C < 2 the two poles are complex. At (" = 2, the two solutions coincide. 
To see this, let us write the set of equations satisfied by the PA (j4.4ip as: 



F'o - 


^2 


zyZA 


4Lio = 


\zy ZA/ 


U = 


2 ZyZA . ZA 

= ^o' , log, • 

ZA - Zy Zy 



(4.43) 

The first of these equations can be used to determine the value of the residue 
r^ in terms of zyZA- In order to analyze the other two, let us first assume 
that both poles zy^A are real. In this case, they also have to be positive or 
else the integral I-^ will not exist because it runs over all positive values of 
Q^. Let us now make the change of variables 

zy = R{l-x) , ZA = R{l + x) . (4.44) 

The condition zy^A > translates into R > 0, \x\ < 1. In terms of these new 
variables, the second and third equations in ()4.43p can be combined into 

C=-log\±^, (4.45) 

X 1 — X 

where the definition ()4.42p for (" has been used. With the help of the identity 
log(l+x/l— x) = 2 th~ X (for |x| < 1), one can finally rewrite this expression 

"^ 2 

C = - th~^x , (x real) (4.46) 

X 

which is an equation with a solution for x only if C > 2. Once this value of 
X is found, the value of R can always be obtained from one of the last two 
equations (|4.43p and this determines the two real poles zy^A from ()4.44p . 
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On the other hand, when ^ < 2, Eq. ()4.46p does not have a solution. 
However, according to ()4.42p . ( can also be smaller than 2. In order to 
study this case, we may use the identity th~^(z y) = i tan~-'^(y) to rewrite 
the above equation ()4.46p in terms of the variable x = i y {y real) as 

2 

C = - tan^^y , (y real). (4.47) 

y 

One now finds that this equation has a solution for y when C ^ 2. In this 
case the poles of the PA ()4.39p are complex-conjugate to each other and 
can be obtained as zv,a = R{^ ^ i y)- These poles, obviously, cannot be 
associated with any resonance mass and this is why this solution has been 
discarded in all resonance saturation schemes up to now. However, from 
the point of view of the rational approximant ()4.39p there is nothing wrong 
with this complex solution, as the approximant is real and well behaved. 
From the lessons learned in the previous section with the model, there is 
no reason to discard this solution since, as we saw, rational approximants 
may use complex poles to produce accurate approximations. Therefore, we 
propose to use both the complex as well as the real solution for the poles 
zv,A, at least insofar as the value for ( ^ 2. In this case we obtain, using 
the values given in Eqs. ()4.40p . 

(C > 2) , r^ = -(4.1 ± 0.5) X 10~^ , zy = (0.77 ± 0.10)^ , za = (0.96 ± O.^l^jlS) 
(C < 2) , r^ = -(3.9 ± 0.1) X 10^3 , zv = z*^ = (0.81 ± 0.04)^ + i (0.25 ± Q.^ap) 

in units of GeV for r^, and GeV for zy^A- The two solutions in Eqs. 
()4.48|4.49p have been separated for illustrative purposes only. It is clear 
that they are continuously connected through the boundary at (^ = 2, at 
which value the two poles coincide and zy = za — (0.85)^. The errors 
quoted are the result of scanning the spread of values in ()4.40p through the 
equations ()4.43p . 

With both set of values in Eqs, ()4.48l I4.49p . one can get to a prediction 
for the chiral and OPE coefficients by expansion in Q^ and l/Q"^, respec- 
tively. These expansions of the PA can be done entirely in the Euclidean 
region Q^ > 0, away from the position of the poles zv,a, whether real or 
complex. Recalling the notation in Eq. (|4.14p . the above P2 (Q^) produces 
the coefficients for these expansions collected in Table 14.71 The values for 
the OPE coefficients C_4^_6,-8 in this table are compatible with those of 
Ref. |117] . after multiplying by a factor of two in order to agree with the 
normalization used by these authors. However, the spectrum in our case is 
different because of the complex solution in (|4.49p . As we saw in the pre- 
vious section with a model, this again shows that Euclidean properties of a 
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C4 


Ce 


Cg 
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C-6 


C-8 


-Fo' 


-4Lio 


-43 ± 13 


81 ±53 


-145 ± 120 


-4.1 ±0.5 


6±2 


-7±6 



Table 4.7: Values of the coefficients C2k in the high- and low-Q^ expansions 
of g2 ni/j(-Q2) in Eq. (H3iD in units of lO'^ GeF^-sfc, j^g^all that 
C_2 = 0. 

given Green's function, such as the OPE and chiral expansions, or integrals 
over Q^ > are safer to approximate with a rational approximant than 
Minkowskian quantities, such as resonance masses and decay constants. 



4.4.1 Improving Cgy 

Once the rational approximant is known, upon reexpansion around Q^ = 
0, higher order unknown coefficients of the chiral expansion may be pre- 
dicted, as shown in the left side of Table HT71 lEB], [88], [87]. If the rational 
approximant is a better description of the original function than the partial 
sums of the chiral expansion, one may expect this prediction to be reliable. 
In a real QCD case, although the LEG Liq is pretty well known [n\ I119J . 
this is not so for the following coefficient Csj- It is therefore important to 
obtain a new determination of this LEG with its associated error. 

Using P2 (Q^) '^6 obtained an estimate of the O(Q^) term in the ex- 
pansion (|4.7|) , or C4 in Table 14. 7[ which translates into the value Cs7 = 
(5.4 ± 1.6) X 10~^ GeV~^, C4 = — 8C87. In the present section, we would 
like to reassess this particular value with a more complete analysis. The 
same can be done for higher order coefficients, e.g. 0{Q^), but with less 
precision as higher the coefficient. 

Improving the precision of our prediction means going further among the 
sequence of approximants. In order to be able to construct that sequence of 
rational approximants it is of course crucial to have enough number of inputs. 
Since PAs are constructed from the coefficients of the Taylor expansion (j4.7p 
one immediately faces an obvious difficulty. Since what one wishes is an 
estimate of Cs7, only the two coefficients Fq and Liq may be used. With 
these two coefficients as input, the only PA vanishing at large Q^ is Pf and 
gives a prediction for Cgr = (7.1 ± 0.5) x 10~^ GeV~^. Actually, another 
possibility is the Pq but this case is just the Taylor expansion itself and, of 
course, will give the prediction Cgr = 0. Since we know from section 14.31 
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that the best approximant has to fahoff as Q~^, and P^ only falls off as Q^^ 
which is too slow as compared to Eq. ()4.5p . is necessary to consider more 
general rational approximants than the standard PAs. 

The model confirms that one may estimate the unknown LECs with 
PTAs and PPAs where, in the first case, the physical masses were chosen 
in increasing order, i.e. Mi < M2 < M3... For instance, with the PTA 
Tj||,2(Q^) we could see that one has a good prediction for the term of 
0{Q'^^'^'^^^') in the low-Q^ expansion, which is the first one not used as 
input, with a precision which improves as the order of the approximant, 
M, increases. Furthermore, the accuracy obtained for the unknown coeffi- 
cients of the Taylor expansion is very hierarchical: the accuracy obtained 
for the term ©(Q^CM+i)) ig better than that for the term of ©(Q^C^z+s))^ 
and that better than for the term ©(Q^' •*), with a quick deterioration 
for higher-order terms. The case of PAs follows the same pattern. As to the 
description of the spectrum, we found that PAs also reproduced the values 
for the residues and masses in a hierarchical way: while the first masses and 
residues are well reproduced, the prediction quickly worsens so that the last 
pole and residue of the PA has no resemblance whatsoever with its physical 
counterpart. The same is true for the residues of a PTA (since the masses 
are fixed to be the physical ones by construction). 

In working with PTs we have seen with the help of the model in section 
14.2.31 that using the residue of the heaviest resonance as input is not a 
good strategy (see, for example. Tables 113] and 231) • Consequently, in the 
following PTs we will not use the residue of the last (and even next to last) 
resonance as input. 

Based on the above, one can envisage the following strategy for getting 
a sequence of estimates for the 0{Q'^) LEG Cgy. Assuming that the vector 
and axial-vector meson masses stay approximately the same in the large- A'^c 
and chiral limits, one can use their values extracted from the PDG book 
[5] to construct several PTAs. We think that this assumption is reasonable 
for both limits. First, for the chiral limit, this is because the up and down 
quark masses are very small, [5]. Second, for the large- A'c limit, there is 
a non negligible amount of phenomenological evidence in favor of the p 
meson being a qq state |135j . Besides, the success in the spectroscopy of 
the quenched lattice results for the lightest vector mesons is also suggestive 
that 1/Nc corrections may not be very large |136)^i . Therefore, we will use 



^^Be that as it may, whether the assumption is correct or not will ultimately be judged 
by the final results obtained. 
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for the masses 

rup = 0.7759 ± 0.0005 , rup, = 1.459 ± 0.011 , nip,, = 1.720 ± 0.020, m^/// = 1.880 ± 0.030 

TJia, = 1.230 ± 0.040 , rria'^ = 1.647 ± 0.022, (4.50) 

where all the numbers have been expressed in GeV. 

For instance, with only Fq and the masses of the p and oi, one can 
construct the PTA T^{Q^) and predict the value for Lio = (-4.32 ±0.02) x 
10~^, which is not bad when compared, e.g., with Eq. ()4.40p . The small 
error comes from the small error of -Fq since the impact of the error of 
the masses is rather negligible. The next term in the expansion gives the 
following value for Cgy = (4.00 it 0.09) x 10^^ GeV~^ which is similar to 
that obtained in section H^ with the Pade P2 • However, since this value for 
Cs7 comes from the second unknown term in the expansion of T2((5^) rather 
than the first (we predicted also Lio), it is quoted here only for illustrative 
purposes and will not be included in our final estimate. Adding Liq and 
the p' mass to the previous set of inputs one can then construct TI{Q'^), 
which produces Csy = (5.13 ± 0.26) x lO^^ GeV-^. The PTA T| can be 
constructed if one also uses the pion mass difference Eq. (j4.4p and m„/ , 
yielding in this case Cs7 = (5.24 ±0.33) x lO'^ GeV-^. We find the stability 
of these predictions quite reassuring. 

A comment on the quoted error estimates is in order. These quoted 
errors are the result of the propagation of errors from the input via the 
montecarlo method |137j . As such, they do not reflect the intrinsic system- 
atic error due to the approximation itself which will be estimated, at the 
end, as the spread of values obtained with the sequence of different approx- 
imants. On the other hand, the propagation of the error from the input via 
the montecarlo method consists in the following. Taking each input in Eq. 
(|4.40p and Eq. (|4.50p . we have constructed a sample of data with a gaussian 
probability distribution yielding as the average and standard deviation pre- 
cisely the corresponding input value and its quoted error, respectively. For 
each member of this sample, the rational approximant is then constructed 
and, upon reexpansion, the LEG is obtained. The distribution of the differ- 
ent values for Cgy so obtained happens to be also gaussian to a very good 
approximation. Therefore it will have an average value X and a standard 
deviation Y which are then used to quote the result for Cg? as X ±Y. 

To be able to construct further rational approximants one needs an extra 
assumption. Although, as we have emphasized above, the residues of the 
heaviest poles in a rational approximant do not come out anywhere close to 
the corresponding physical decay constants, this is not true for the lightest 
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t;^ 


inputs 


CS7 


TO 


/o ; rup, ma 


4.00 ±0.09 


n 


fo, Lio ; nip, ma, m^i 


5.24 ±0.33 


^2 (a) 


/o, Lio, 6M.^ ; mp, ma, m^i, m^i 


5.25 ±0.10 


Tp2(5) 


fo, Lio, Fp ; mp, ma, m^,, m^, 


6.00 ±0.15 


r|p3(a) 


fo, Lio, Fp, 6M^ ; mp, ma, m^, m^, m^i 


5.78 ±0.21 


r^Sib) 


fo, Lio, Fp, Fa ; mp, ma, m^,, m^, m^,, 


6.24 ±0.38 


n 


fo, Lio, Fp, Fa, SM^ ; mp, ma, m^', m^>, m^", m^/// 


6.03 ±0.38 



Table 4.8: Set of inputs used for the construction of the different Fade Type 
Approxiniants in the text. The last column is the summary of predictions 
for C87 explained in the text. 



ones. In particular, section [4.2.31 shows that the value of the residue for the 
first pole in a PTA could reproduce the exact value in the model with very 
good precision if the order of the PTA is high enough and, more importantly, 
it is improving as the order of the PTA grows. Consequently, if we are willing 
to use the decay constant Fp, and perhaps also the Fa^, one can go for the 
construction of higher PTAs. These two residues can be gotten from the 
decays p — >• e+e" and oi — )■ 7r7, respectively, and their values are |117j 



0.156 ±0.001 F, 



ai 



0.123 ±0.024 



(4.51) 



in GeV units. 

For instance, using Fq,Lio,5M^ and Fp, as well as the five masses 
iTT-p, m-ai , nipj,m,ai and mipii, one can construct the PTA T^. Upon expanding 
this approximant, one obtains the value Cgy = (5.78 ± 0.21) x 10^'^ GeV~^. 
Alternatively, one can also use FQ,Lio,Fp and only the first four masses to 
construct a T| approximant, which is different from the other T| considered 
above. The value obtained for Cgr, i.e. Cgy = (6.00 ±0.15) x 10"^ GeV'^, is 
nevertheless very similar, which again brings confidence on the prediction. 

In this way we have constructed a variety of rational approximants which 
we have listed on Table 14.81 in increasing order of the degree in the denom- 
inator, together with the set of inputs used and the prediction for Csy- We 
have gone all the way up until the Tg, with the six masses listed on (j4.50p . 
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Figure l3.13l shows the prediction for the LEG Cst from the corresponding 
rational approximant shown on the abscissa, upon expansion around Q^ = 0. 
We also included the previous result obtained with the PA P2 R-^f. [86] , [88] , 
but with the montecarlo method for the treatment of errors. As one can 
see, the stability of the result is quite striking. After averaging over all 
these points, we obtain as our final result in the large- A'^c limit, 

^87 = (5.7 ± 0.5) X 10"^ GeV-^ (4.52) 

The error in Eq. ()4.52p is mainly dominated by the error on the input for 
Lio in Eq. (|4.40p and is rather insensitive to the errors on the other inputs. 
For instance, one could increase the error on Fq to 5 MeV in Eq. (|4.4U|) . or 
the error on rup to 50 MeV in Eq. (|4.50|) . or the error on (^M^ to 0.5 MeV 
in Eq. ()4.40p , without falling out of the error band given in Eq. ()4.52p . 

For comparison, we also show in Fig. 14.131 the result of several other 
estimates for this LEG. Reference [99] (shown as 'A') uses the residues in 
Eq. ()4.5ip and the p and a\ physical masses to construct, in effect, what 
we could call the PTA T\ to Q^'^y-A- The difference between this result 
and ours stems from the fact that this rational approximant falls off like a 
constant at large Q^, unlike Eq. (j4.5p . Also, as we have already emphasized, 
the use of the physical decay constant Fa-^ (j4.5ip in a rational approximant 
which has the a\ as the heaviest pole is a potential source of error. 

Reference [5B] also obtains an estimate for Cgy (shown as 'B') based on 
the construction of a rational approximant which effectively coincides with 
the PTA T^ but using the physical value of F^ = 92.4 MeV [5j instead of 
the value of Fq in Eq. (|4.40p . Had they used Fq, the result would have been 
lower, and would have agreed with the value we mentioned in the paragraph 
right after Eq. (|4.50p . Therefore, our comments on the Pade Type Y\ found 
in that paragraph also apply to this determination in [56]. 

One can also get still another estimate for Cgy from the PTA Y\ in [56] by 
assuming that the a\ mass in the large- A'c limit is not approximated by the 
physical value in Eq. (|4.50p . but by a value which comes from the radiative 
pion decay saturated with the p and the a\. This value turns out to be 
vfiax ~ 998 MeV |118j . This lower number for the a\ mass is the reason for 
a higher value for Cgz than that obtained in [56], and is shown as 'G' in Fig. 
14.131 However, there is no compelling reason to associate this different mass 
of the a\ with the large- A'^c limit as we have discussed in sec l4.2.3l In fact, 
our results show how similar values for Cgy can be obtained with the physical 
masses of the mesons used for the poles. Moreover, one of the advantages 
of our method is that one can get a rough idea about the systematic error 
involved by looking at the dispersion of the values obtained. 
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Figure 4.13: Prediction for Cgy in the large- A'^c limit from the PA Pg ™ Refs. 
[86], [87], and the different PTAs discussed in the text and appearing in Table 
1. For comparison we also show the estimate from Refs. [991 [56l |TT8] , which 
we label 'A', 'B' and 'C (resp.). 



Recently, two more estimates for Cgr appear in the literature. The first, 
Ref. [63] , is based on Resonance Chiral Theory [42 ^ W\] and obtains a esti- 
mate Cgy = (3.9 lb 1.4) • 10~^ by studying the Next-to-leading order term in 
the large- A'^c limit of the function Q^IlLRiQ^)- The second one, Ref. [138j . 
based on r-decay analysis, gives Cg^ = (4.89 ± 0.19) • 10"^ trough a direct 
fit to the experimental data using a sum rule approach. 

Of course, in the large- A'^c limit CgT does not run with scale whereas in 
the world at Nc = 3 it does. This is an additional source of systematic 
error in the result ()4.52p . However, phenomenological evidence as well as 
theoretical prejudice [44] suggests that a reasonable guess for this systematic 
error may be obtained by varying the scale in C87(/u) between the range 
0.5 GeV ^ /i ^ 1 GeV (compare with Liq in (j4.40p l. Using the running 
obtained in Ref. |139j . this error turns out to be ~ 30 per cent, right in the 
ballpark expected for a typical 1/Nc effect. This systematic error should be 
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PartialPades 

Figure 4.14: Prediction for Cs7 in the large- A'^c limit from the different PPAs 
discussed in the text. 
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added to our large- A'^c result in Eq. ()4.52p in order to obtain an estimate for 
CsrifJ- ~ 0.7) in the real world. In this case, all the different results in Fig. 
14.131 can be encompassed by this error. 

We would like to finish this section by recalling that PAs and PTAs are, 
in a way, two extreme versions of a rational approximant. While in the 
latter all poles are fixed at the physical masses, in the former the poles are 
left free, and they are obtained by demanding that the expansion around 
Q^ = reproduces that of the original function to the highest possible or- 
der. Besides these two rational approximants, there are also the Partial 
Pade Approximants [86l I92j (introduced in section |4. 2. 3p which, from a cer- 
tain point of view, lie half way between PAs and PTAs. As the poles of 
its denominator are only partially preassigned, there is no reason why, in 
general, the poles of a Partial Pade should come out to be purely reaJlj, 
unlike those of a PTA, which are of course real by construction. We have 
constructed seven of these Partial Pades, see Figure [3. 141 with a polynomial 
in the denominator up to fifth order in Q^. In some of the cases the poles 
were actually complex, as it was also the case of the PA P2, Eqs. (J4.17I4.19|) . 
However, the results obtained for Cgr are almost identical to those in Fig. 
14.131 although with errors which are somewhat larger. This feature rein- 
forces the stability of the result shown in Fig. 14.131 and gives us reassurance 
about the reliability of our result. Finally, we would like to mention that 
predicting 0{p^) LECs may also be another straightforward application of 
this method. 



^^ Although, when complex, they always come in complex conjugate pairs. This just 
means that, in general, the poles of a rational approximant are not necessarily physical. 
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Chapter 5 



Pade Theory applied to 
Meromorphic functions of 
Stieltjes type 



5.1 Introduction 

In the present chapter we will steer things in a different direction. Up to 
now, we have focused our attention on meromorphic (but not Stieltjes) func- 
tions, but as we have already said in the introduction Pade Approximants 
can be successfully applied also to a Stieltjes functions. In the following 
chapter, we will present two different examples of that kind of application. 
The first one, related with the Linear Sigma Model and the prediction of 
a resonance mass (all the details can be found in [82], |140j . [85]). and the 
second one related to the vacuum polarization function of a heavy quark 
(see [83] for further details). 

One of the mean advantages of PAs with regard to Taylor expansion is 
its range of applicability (see the examples of sec. 13.11 and the underlying 
motivation of chapter 3) . In chapters 2 we saw that PA can deal with mero- 
morphic functions when trying to predict LECs. What happens, however, 
when instead of a set of poles, your function has a brunch cut? When the 
function is Stieltjes Eq. (j3.10p . i.e., its spectral function is positive defined 
(see de discussion after Eq. (|3.10p ). the convergence of the PAs are guar- 
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anteed in a region out of the branch cut. When its spectral function is not 
positive defined (either meromorphic) , we do not know of any mathematical 
result assuring the convergence of the PA sequence. Nevertheless, in few of 
these cases, PAs are surprisingly useful even when the convergence of PA 
are not properly guaranteed for any theorem. There is also a particular 
function that its analyticity is defined by a cut and a finite set of poles, the 
so-called meromorphic functions of Stlietjes type, for which also convergence 
theorems are known and can be applied |92[ I14H I142| . In the first section 
of this chapter we will go into these kind of functions and in section 15.51 we 
will concentrate our efforts on an example of a Stieltjes function. 



5.2 The Linear Sigma Model case 

The physical motivation for these studies is based on the need of further 
ingredients when your Effective Field Theory (EFT) is unable to describe 
the resonance region. One example of that situation occurs in Chiral Per- 
turbation Theory. As an EFT is a very useful tool for the description of 
low-energy physics [m [151 ES] but its range of validity is just a few hundred 
MeV, breaking down as one approaches new states not included in the EFT. 
Several works, then, have tried to extend the range of validity of the EFT 
by means of the unitarization of the low-energy amplitude. Unitarization 
methods have been applied extensively to Quantum Chromodynamics and 
xPT |143j , where the issue of the scalar resonances is particularly interesting. 
However, their range of applicability is wider; for example they have also 
been applied to W^Vl^-scattering |144t 1145] . Among the most usual ones, 
Pade Approximants were already used before the nineties [146t I147J and 



more recently one of its variants, the Inverse Amplitude Method (lAM) [1481 
I149^ I150t I15H I152J , although some remarks and criticisms on the reliability 
of these and other unitarization methods have been raised [147t I153t I154J . 
In the present section we study up to what point one can rely on them to 
describe the resonant properties of the theory. By means of a couple of 
models (the Linear Sigma Model and a vector resonance model), we show 
that these unitarization procedures may lead to improper determinations of 
the resonance pole position (masses and widths) \l5b\ 1156] . 

The starting point of our analysis is, therefore, a model where the prop- 
erties of the resonances are known. Then we derive and unitarize the corre- 
sponding low-energy amplitude. The predictions for the masses and widths 
obtained from the PA sequence P^~^ , {J > — ijj quickly converge for 



^While one can use J > —1 an obtain similar results, we choose the diagonal case to 
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Stieltjes functions, as A^ is increased. We show that also for a meromorphic 
function of Stieltjes type (such as the first line of Eq. ()5.14p the convergence 
is also obtained thanks to a theorem of convergence [92^ I157t I14H I142J . 

We also show that the same convergence patterns does not occur for any 
PA sequence, i.e., the restriction J > — 1 should be fulfilled. Actually, when 
applied to a Stieltjes function it is known that all the sequences P^^ , {J > 
— 1), have guaranteed convergence (see sec Jci.l.ljj) . For the case J < —1, your 
sequence is not anymore protected. In particular, could happen your PAs 
have undesirable poles. Sometimes these poles will be complex-conjugate 
poles in the form of defects or will tend to infinity when increasing the order 
of the PA. Eventually, they could lie on the real axis of the complex plane 
where the function has no branch cut. One of this cases is when J = 1 — N, 
the sequence P^ (referred to as lAM in some works |147l 11481 11491 11501 I15H 
11521 1154] ). In the present section we will study a Stieljes function and also 
a meromorphic function of Stieltjes-type with both P^ and P^ sequences. 
The results are compared and found to be quite different from those of the 
original model. For simplicity the chiral limit will be assumed all along the 
section, but this will not alter the main results. 



5.3 One-loop Linear Sigma Model and its unita- 
rization 

5.3.1 Unitarization of the xPT amplitude 

To proceed we will first compute the P^ sequence. To show that certainly 
the Inverse Amplitude Method [I AM) is a P^ sequence we will first define 
lAM. 

The lAM provides an amplitude that is unitary not only at the per- 
turbative level but exactly. Considering the partial waves projection 

tj{s) = -^J dcose Pj{cos0)T{s,t,u)^, (5.1) 

where T{s, t, u) are the isospin amplitudes and 9 the scattering angle in 
the vrvr center-of-mass rest frame. Pj are the Legendre polynomial. In the 
elastic limit one has for s > 

Imi(s) = |t(s)p, (5.2) 



exemplify the procedure. 
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where the indices IJ are assumed, i.e. t = tj. 

This relation can be reexpressed as a relation for the inverse amplitude: 



Imtis)"^ 



(5.3) 



Thus, the imaginary part of t{s)~^ becomes completely determined and one 
only needs to specify the real part Ret~^. The I AM relies then on a low- 
energy matching to x PT (since lAM has been applied in the literature to 
that EFT) with 



VVt ~*(2) [l-*(4)A(2) + •••] , 



(5.4) 



in order to fix the unknown part of the amplitude. For the first partial 
waves tj{s), with IJ = 00, 11, 20, one finds the following 0{p^) amplitudes, 

to(s)(2) 
*l(s)(2) 

*o(s)(2) = - ooli^a ' (5-5) 



167rF2 ' 

s 

967rF2 ' 

s 
~327rF2 ' 



and at 0{p'^ 



tl{s){i) = to(s)(2) X ^ 

-2s 



*o('S)(4) = io(«)(2) X 



3M2 



1- 

— s 

1 



5 



2647r2 
1 + 



''^'-M + ^-sj + f l+o'^^ 



(In — In 
M2 M2 



26 



+ 0(5^ 



5 /9, — s 11, s 163 \ „, 2n /^ ^n 



where we have used the relation 2gF^ = M^ |26J. 
Thus, at 0{p^), one has the unitarized amplitude, 

t( 



Ham 



K2) 



1 



H4) 
*(2) 



(5.7) 



This expression is sometimes also known as a P/ Fade Approximant of the 
partial-wave amplitude. The JAM has been also extended up to 0{p^) by 
means of what is sometimes named as a P2 approximant |152|, I154| : 

t{2) 



hkM 



«(4) 



(4) 



*(2) *(2) V*(2) 
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2 • 



and the P^ Fade at 0{p^ 



^lAM 



t(2) 



I Hi) 

t(2) 



«(6) 



7^+7^ 



H2) 



'(4) 



■(2) 



t(2) 



2tif,\t 



(6)^4) 



(2) 



However, we want to remark that i(>s)j^]y[ is not a PA in the variable s: 
It is not a rational approximant since it also contains the logarithms from 
the pion loops, as one can see in the specific expressions of the correspond- 
ing amplitudes Eq. ()5.6|i4 Thus, strictly speaking no theoretical argument 
ensures the recovery of the physical amplitude. Only in the tree-level limit 
^(•^)lAM becomes a PA. In any case, we will see that both the whole and 
the tree-level lAM amplitudes are unable to reproduce the original partial 
waves in the resonance region. 

To proceed given the O(p^) and 0{p^) amplitudes, ti2) ^'^'^ ^(4)) Eqs. 



(|5.5p and ()5.6p . is possible to extract the poles of the corresponding t(s)iAM 
for the LSM, satisfying 1 = t(s)(4)/t(s)(2) at s = Sp. For the firsts partial 
waves tj{s), one finds (see Ref. |.82j): 
IJ=00 



_6_ 

IT' 



-Mi 



1 + 



9 /193 



2647r2 V 3 



251n' 



_6_ 
IT 



ISivr + 0{g^ 



IJ=11 



IJ=20 



:Mi 



-Mi 



1 + 



1 + 



487r2 



26 



+ ivr + 0{g' 



g /163 ^3 llfvr 

— ^ h 5 In - H 

24^2 \^ 24 2 4 



+ 0{g' 



, (5.^ 



(5.9) 



. (5.10) 



These are the poles that appear in the unphysical Riemann sheet as one 
approaches from upper half of the first Riemann sheet. There is also a 
conjugate pole at s* if one approaches the real s~axis from below. 

The first thing to be noticed is that poles appear in the 7J = 11 and 
20 channels even for small values of g, contrary to what one expects in the 
LSM, where no meson with these quantum numbers exists. Furthermore, 
these "states" are not resonances, as they are located on the left-hand side 
of the complex s-plane, out of the physical Riemann sheet, and carrying a 
negative squared mass. 



See Ref. [H] for farther details about the calculations of these amplitudes. 
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As for the IJ = 00 channel, one finds a resonance with pole mass and 
width, 



0{g' 



MJ 121 levr ' -^- '■ (^-1^) 

These expressions have to be compared with the original ones in the 
LSM, which are given in Ref. |82j : 
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I. (5.12) 



The I AM predictions for M^ and MpTp result, respectively, 40% and 
% smaller than the original ones in the LSM, shown in Eq. ()5.12p . The 
JAM poles remain badly located even in the weakly interacting limit, so 
this failure cannot be attributed to non-perturbative effects. In the limit 
when 5 — )■ and Mg- is kept fixed one finds that the poles predicted in 
all the different channels fall down to the real s-axis. This points out the 
low reliability of this particular method in order to recover the hadronic 
properties of the theory from its effective low-energy description. 

We are left with just tree-level amplitudes and the expressions become 
greatly simplified. Due to the smoothness of this limit, it will be assumed in 
the next analysis of higher order Fade Approximants P^ and in the study 
of the vector model in Section EI 



5.3.2 Higher order Pade Approximants for tree-level ampli- 
tudes 

In this subsection we consider higher order Pade Approximants to the 
partial wave amplitudes, with the hope that this will provide some insight 
on the nature of the unitarization process discussed above. We will see that 
the PA sequence associated with the lAM does not converge properly, and 
that diagonal and paradiagonal sequences, sec J3.1l are much more suitable 
for this purpose. 
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Tree-level PAs in the LSM 



In order to be able to handle the amplitude at higher orders, we will 
consider the vrvr scattering at tree-level. This is equivalent to working in 
the limit g <^ 1 and keeping just the first non-trivial contribution in the 
g expansion. Thus, the vrvr-scattering is determined in the LSM by the 
function 

Ms,t,n) = -^^j^^-, (5.13) 



F2 M2 - S ' 
By means of the partial wave projection in Eq. (|5.ip . this provides 
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(5.14) 



The In [l -|- sjM'^ logarithms come from the partial-wave projection of the 
tree-level exchanges of resonances in the crossed channel. They have abso- 
lutely nothing to do with the logarithms of the xPT amplitudes in Eq. ()5.8p , 
which come from the vrvr loops. 

From a mathematical point of view both i\ and tg are Stieltjes functions 
(see Eq l3.10l and below) with a brunch cut in —Ml > s > — cxd. In that 
situation, the convergence of a Fade sequence -P/y^ , J > — 1 is guaranteed 
(see sec l3.1.3| . However, Iq is what is called a meromorphic function of 
Stieltjes type [5 ^ll57tll41|,ll42] . that is a function that can be decomposed 
in two parts, a Stieltjes function and a rational function. This last piece 
must have a finite number of poles all located out of the branch cut of the 
Stieltjes part. In our case, ig has a simple pole at s = Ml and a brunch 
cut -oo < s < -Ml. Besides, Refs. [92l |1571 EU US] demonstrate a 
theorem of convergence for Fade Approximants for diagonal Fade sequences 
[Ha 157] and also for P^^-^ , J > -1, [HIl [92] that wih we used in this 
section. 

The relation between Stieltjes functions and its corresponding Fade Ap- 
proximant is already explained in chapter [31 Further details will be studied 
in section 15. 5[ The most interesting part here is the study of t9. who has 



different analytical properties than its counterparts Eq. ()5.14p . From now on 
we will basically focuss on this particular function. 
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At low energies the amplitude becomes 
A{s,t,u) 



s 
F2 



s s^ 



(5.15) 



so the partial waves are given by: 
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(5.16) 



The comparison between the low-energy expansions and the whole result 
provides a first insight of the piece of information that is lost in the unita- 
rization procedure. At high energies, the partial waves contain poles on the 
right-hand side of the s-plane, related to s-channel resonance exchanges, 
and a left-hand cut, related to the crossed-channel resonance exchanges. At 
low energies, both kinds of exchanges contribute equally to the low energy 
couplings, so the crossed resonance exchanges shift the JAM poles from their 
physical value. Although t and u channels are not so relevant in the region 
close to the resonance pole, at low energies they are as important as the 
s-channel. 

The simplest Pade Approximant, Pi, gives the prediction 
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(5.17) 



which agrees with the one-loop calculation from Eqs. ()5.8p - (|5.10p if one 
remains at leading order in g. 



5.3.3 Higher order PAs in the LSM 

Now we will proceed to study higher order PAs. In particular, we will 
employ and compare the P^ and P^ sequences for the study of the tttt 
partial wave scattering amplitudes. We will also comment on PAs of the 
Pat^ , e.g Ppj~ and Pm~^ ■ In the next lines we will focus our attention 



N 



N 
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on the IJ = 00 partial wave, but analogous results are found for the other 
channels. 

Former works pointed out that the PAs and other unitarizations fail to 
incorporate the crossed channel resonance exchanges [15811159] , Nonetheless, 
we will see that as A^ grows, the poles of the sequence P^ actually tend to 
mimic not only the s-channel poles but also the left-hand cut contribution 
from diagrams with resonances in the t and u channels. 

Another advantage of the sequence P^ is that, as already explained in 
sec. 14.31 is the best PA sequence one can use when trying to predict global 
properties of the original function. Fig. 15.11 shows the determination of the 
position of the M^ pole, its residue and the Taylor coefficient an from several 
PAs, all constructed using only the first Taylor coefficients from Eq. ()5.16p . 
Exactly following sec. I4.3l we take the t^{s) w ais^ + a2S^ + • • • aios^^ and we 
study all the possible PA constructed using only these coefficients. These 
lets to a set of nine PA: Pf, P^^P],, Pf , P^, Pi, P|, /f , Pg. The best PA is 
the diagonal one, who gives the best prediction for the three parameters 
we are studying. Also, this diagonal PA gives an extra insight about the 
behavior of to(s) for s — )• oo, Eq. (|5.18|) . Of course, a rational function will 
never we able to predict a logarithm function but, the diagonal PA finds the 
compromise between the 1/s piece and the log(l/s) piece: 



■' ''°'^'^ {l/s)' + Oi{l/sf) (5.18) 



V S / 



S 

The results shown in these plots reassess the conclusions extracted from 
seciOl 

Our results for the sequence P^ are summarized in Fig. 15.21 No conver- 
gence is found with this sequence. In the case of N odd. Fig. I5.2i a. shows 
that the P^ pole closest to M^ does not approach this value even for very 
large N, always remaining a 30% below. The analytical structure of the 
original amplitude (s-channel sigma pole plus left-hand cut) is never recov- 
ered since the P^ PAs always set the poles in the circular pattern shown in 
Fig. 15.21 b. This suggests that the use of further P^ approximants to extend 
the lAM is not the optimal way to proceed, even if we had an accurate 
knowledge of the low-energy expansion up to very high orders as is now the 
case. 

Alternatively, the use of sequences such as P^'^'^ (e.g. P^~^, P/v , P/v^^ • 
seems to be a better strategy. In the following we analyze the sequence Pjy , 
as it ensures the appropriate behavior at high energies \t{s)\ < 1. The P^ 
pole closest to M^ is shown in Fig. I5.3i a. One finds a quick convergence of 



^Imposing unitarity to our PA sequence is a constrains that lies on a physical motivation 
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Figure 5.1: Determination of the position of the M^ pole (top), the first 
residue (center) and the Taylor coefficient S22 (bottom) from several PAs, 
approximating the function to(s), Eq. ()5.16p . See the text for details. 
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Figure 5.2: Top: Position of the nearest pole to M^ for the first PAs of the 
form P^ with N odd (for even N ah the poles are complex). Bottom: Poles 
of the Pqi in the complex plane. 



the sequence: P/ reproduces the sigma pole a 40% off but P| disagrees by 
less than 1%, Pg by less than 0.1%, etc. Notice that already Pi provides a 
much better description than Pq^, although one includes far more low-energy 
information in the latter. All this points out the sizable discrepancy of the 
first element of the sequence (P/) with respect to the original amplitude. It 
also indicates that the P^ PAs do not produces a serious improvement. On 
the contrary, the P^ sequence provide a far more efficient strategy with a 
quick convergence. 

Likewise, Fig. 15. 21 b shows how the P^ PAs are unable to recover the 



but as far as concern only with PAs, it is not a necessary assumption to be able to establish 
a converge sequence as explained in sec l3.ll 
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Figure 5.3: Top: Location of the closest pole to M^ for the first Pj^ PAs. The red 
line showns the exact value s/M^ = 1. Bottom: Poles of P|q . 



analytical structure of the original amplitude, whereas the P^ sequence, 
besides providing the isolated pole of the sigma, tends to reproduce the left- 
hand cut as A'' increases, as expected since t^ is a Stieltjes-type function. 



The poles of P^o ^^® plotted in Fig. 15.31 Although a PA is a rational function 
without cuts, these are mimicked by placing poles where the cuts should lie. 
The P|q has one isolated pole near M^ (with an accuracy of 10^^") and 
nineteen poles over the real axis at Sp < —M^, i.e. on the left-hand cut of 
the original function. As N is increased, the number of poles lying on the 
branch cut increases too. 

A remarkable feature found for the first P^ approximants {Pi, P2, -P3) 
is that they obey exact unitarity, as it happened with the lAM sequence 
pi 
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In several situations the PAs set all poles over the left-hand cut position 
and one isolated pole that approached M^ when N — t- oo. In the worst 
cases, in addition to this we found extraneous poles that either moved away 
as N increased or they tend to be canceled by nearby zeros at A^ — )• oo. 

As an amusement inspired by the Montessus^s theorem and the fact 
that there is only one ressonance a, we have also probed the PA sequence 
P^ which has the same number of inputs as the P^ has for a given N . 
Also has the advantage that to increase one order on the sequence we just 
need one extra parameter while in the P^ case, we need two extra pa- 
rameters. In this new case we have found convergence in both LSM and 
the resonance model presented in the following subsection but slower than 
the Pf^ . For instance, the prediction for the M^ for the first P^ are 
-^ = 0.55,1.47,0.73,1.27,0.81... A criticism that can be done to this se- 
quence is its lack of unitarity, and the impossibility to mimic the brunch cut 
of the original function, in contrast to the other studied sequences. 



5.4 Vector Resonance Model 

In order to broaden our analysis, we consider now a model with just vector 
mesons [158] . It could be derived either from the gauged chiral model 



for the couplings 3gpF^ = Mp, or from resonance chiral theory [32] with 
only vectors and the relation 3Gy = F^. The Trvr-scattering is given in this 



model by 



M2 

A(s,t,u) = — ^ 



+ 



M^-t Mj-u 



(5.19) 



The study of the IJ = 11 partial wave leads to the same conclusions 
found for the LSM. It is at first sight remarkable that, on the contrary to 
the previous case, one recovers Sp = M^ from the first-order approximant 
P^. However, the sequence P^ already worsens at A^ = 2, where the two 
complex-conjugate poles are located at Sp = (0.71ib0.96i) M^ on the physical 
Riemann sheet. On the other hand, Pj^j exactly recovers Sp = M^ for any 
odd A^. For N even, the prediction from P| is a 30% off, but one has again 
a quick convergence to Sp = Mp as N increases: Pf disagrees by less than 
0.1%, Pg disagrees by less than 10~^, etc. 

Furthermore, the P^ and P^ sequences produce, respectively, the same 
analytical structure of poles found for the LSM. This is, P^ generates the 
circular analytical structure of poles of Fig. I5.2[ b. and the sequence P^ 
places one pole at Sp ~ M^ and the remaining ones reproducing the left- 
hand cut in analogy to Fig. 15. 31 b. 
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5.5 The Vacuum Polarization of a Heavy Quark: 
Motivation 

In this second example, the physical motivations is the study of the vac- 
uum polarization function of two electromagnetic currents, in the framework 
of heavy quark physics. This study requires high order perturbative calcu- 
lations which, because of the obvious need to keep a nonzero mass m for 
the quark, become extremely difficult to perform. This is why, while the 
0{a^) and 0{al) contributions have been known for a long time |16Uj . state 
of the art calculations can only produce a result at O(a^) in the form of 
an expansion at low energies (i.e. q^ = 0), at high energies (i.e. q^ — )• oo) 
or at threshold (i.e. q^ = 4m?), but not for the complete function, which 
is still out of calculational reach. In this circumstances, it would of course 
be very interesting to be able to reconstruct this function by some kind of 
interpolation between the former three expansions. 

After the work in refs. |16H I162"t 11631 fTM} I165J . it has become customary 
to attempt this reconstruction of the vacuum polarization function with the 
help of Pade Approximants. Since these approximants are ratios of two 
polynomials in the variable g^, they are very suitable for the matching onto 
the low-energy expansion. This is so because this low-g^ expansion is trul 
an expansion in powers of q'^, as a consequence of the finite energy threshoL 
starting at 4m?. However, it is clear that they cannot fully recover the 
nonanalytic terms which appear, e.g., in the form of logarithms of q'^ in the 
expansion at high energies (or at threshold, where there is also a squared 
root). Therefore, what is really done in practice is to first subtract all these 
logarithmic pieces from the full function (impossible to match exactly with 
a Pade Approximant) with the help of a guess function with the appropriate 
threshold and high-energy behavior, and then apply Pade Approximants to 
the remaining regular expression. In fact, this is done after a conformal 
mapping from z to u whereby all the (cut) complex plane is mapped into 
a circle of unit radius. In this way, the authors of Ref. |166j were able 
to compute, e.g., the value of the physical constant K^"^' appearing in the 
0{a'^) expansion of the vacuum polarization at threshold (Eq. ()5.35p ). which 
has not yet been possible to obtain from a Feynman diagram calculation. 

Although this result is very interesting, the construction is not unique. 



In perturbation theory, this is true so long as purely gluonic intermediate states are 
not considered. Beyond perturbation theory, the threshold occurs at 4m^, where m,^ is 
the pion mass. 
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As recognized in Ref . |166j , some amount of educated guesswork is required 
in order to resolve the inherent ambiguity in the procedure. For instance, 
a certain number of unphysical poles are encountered, and some additional 
criteria have to be imposed in order to decide how to discard these poles. 
Since the resulting ambiguity leads to a systematic error which needs to 
be quantified, this error is then estimated by varying among several of the 
possible arbitrary choices in the construction. Although all these choices are 
made judiciously and in a physically motivated manner, it is very difficult to 
be confident of the error made in the result, which obviously has an impact 
on the value extracted for the constant K^'^> . 

In this section we would like to point out that, regarding the vacuum 
polarization function, one can do away with all the above ambiguities when 
applying Fade Theory. 

Since the vacuum polarization function is a Stieltjes function (as we will 
see in the next subsection), the convergence of all the Fade Approximants 
P]y~^ , {J > —1) is guaranteed. The result of this theorem together with the 
fantastic amount of information obtained on the Taylor expansion around 
q^ = 0, for which 30 terms are known [167j . will allow us to predict a value 
for K^'^' Eq. ()5.35p . As it turns out, our result is very close to that of Ref. 
[166j . although slightly smaller. 

5.5.1 Pade Approximants and the Vacuum Polarization func- 
tion: Definitions 

The vacuum polarization function n(g^) through the correlator of two 
electromagnetic currents j^{x) = qix)j^q{x), is defined as: 

{9^.uq^ - q^.qu) U{q'') = - i ^ dxe^^" (0 |r j^(x)>(0)| 0) , (5.20) 

where q'^ is the external four-momentum. The optical theorem tells us that 
the e~^e~ cross section is proportional to the imaginary part of n(g^). As a 
result, Imn(g^) is a positive definite function, i.e. 

Imn(g2 + ie) > , (5.21) 

a property which will become crucial in what follows. 

In perturbation theory n(g^) may be decomposed to 0{a'l) as 

n(g^) =nW(g2) + (^) n«(,2) + (^)'n(2)(,2) +o{al). (5.22) 



vr / V vr 
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For definiteness, a^ denotes the strong coupling constant in the MS 
scheme at the scale fi = rupoie, but this is not important for the discussion 
which follows. Equation ()5.22p will be understood in the on-shell normaliza- 
tion scheme where a subtraction at zero momentum has been made in such 
a way as to guarantee that 11(0) = 0. In that way, the vacuum polarization 
in Eq. ()5.22p satisfies a once subtracted dispersion relation (see ApplXlfor 
details), i.e. 



n(g2) = g2 / _^ -Imn(i + ie) . (5.23) 



2 / dt 2:^ 

t{t — q^ — ie) vr 

Since all diagrams with intermediate gluon states are absent up to 0{a'^), 
the lower limit for the dispersive integral ()5.23p starts, in fact, at a finite 
value given by the threshold for pair production, i.e. 4m^. This fact only 
carries over to higher orders in Og provided these intermediate gluon states 
are neglected. From now on, we will restrict ourselves to the vacuum polar- 
ization in Eq. ()5.22p to O(a^), neglecting higher orders in ag. 
In terms of the more convenient variable 

z^ ^, (5.24) 

4771^ 

one can rewrite Eq. ()5.23p . after redefining u = Am? /t, aqj 

n(z) = z — - -Imn [4.n?u-^ + ie) . (5.25) 

Jn 1 — tlZ — IE 71 

RecalHng that a Stieltjes function is defined as [HOl [901 IIMl ESI [79] 

f^/^ d(t)(u) 
fiz) = / -r^^^ (5-26) 



1 — uz 
where (p{u) is any nondecreasing function, one sees that the identification 

d(j){u) = -Imn (Am^u'^ + ie) du (5.27) 

TT 

allows one to recognize that the integral in Eq. (|5.25p defines the Stieltjes 
function z^^n(z). It is been common in the literature to work with the 
variable w defined as z = 4a;/(l + uf instead of z [MD [Ml [Ml [Ml [Ml 
I166J . This change of variables allows a mapping to the cut z plane into a 
unit circle in the a; plane. In particular, the three expansions at low energies 
(i.e. z = 0), at high energies (i.e. z — )• cxd) and at threshold (i.e. z = 1) 
translate into expansions at a; = 0, a; = —1 and oj = 1 respectively. This 



^We are simplifying the notation by replacing n(4m z) — > ^{z) 
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conformal mapping from uj to z allows the approximation to be valid also on 
the boundary, providing a smooth approximation of the absorptive part of 
the integral. We note, however, that this change of variables invalidates the 
Stieltjes property and then there is no theorem that we know of ensuring 
the desired convergence for the Pade Approximants in the new conformal 
variable. 

As one can see, the positivity property Eq. (|5.21|) is crucial for the 
identification (j5.27|) to be possible. The connection between the imaginary 
part of the vacuum polarization and the cross section a{e~^e~ — )• had) (a 
clearly positive definite quantity) assures that this property is guaranteed 
to all orders of perturbation theory. 

The representation of the function f{z) in Eq. ()5.26p clearly shows a 
cut in the z complex plane on the positive real axis for R < z < oo. For 
the physical function 11(2:), this of course corresponds to the physical cut in 
momentum for Am? < q^ < 00, i.e. the physical case corresponds to i? = 1 
in Eqs. (J5.26l5.27p . Furthermore, just like the function f{z) in Eq. (j5.26p 
has a power series expansion convergent in the disk \z\ < R, so does the 
function 11(2;) in Eq. (j5.25p have a power series expansion convergent in the 
disk \z\ < 1. 

The position of the poles in the Fade Approximant P^~^ (z) (with J > 
—1) accumulate on the positive real axis starting at threshold, q^ = Am?, 
mimicking the presence of the physical cut in the original function (provided 
that the function is a Stieltjes function). When FAs are applied to the 
vacuum polarization, this means, in particular, that there can be no spurious 
poles outside of the positive real axis in the z plane and, consequently, no 
room for ambiguities. Furthermore, the convergence of the approximation 
(and the error) can be checked as a function of A^, as we will see. 

In Eq. (|5.22p . the full functions 'ny''^'{q^) are known. They are given by 
the following expressions [160| : 



n(o)(2) 
n«(2) 
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(5.28) 
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where 



liz) 



C3+4Li3(-u) + 2Li3(n) -8 2Li2(-u) + Li2(u 



2 ln(l + n) + ln(l - u) 



In u . 
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2u \nu 



u^ 



1 
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(5.29) 



and with (^„, the Riemann zeta function, Li„(ti) the Polylogarithm func- 
tion defined as Li„(z) = Xlfcli F^- 

However, the situation with the function n(^-'(g^) is different. In fact, 
n'^'(g^) is only partially known through its low-energy power series ex- 
pansion around g^ = 0, its high-energy expansion around g'^ — )• oo and its 
threshold expansion around q^ = 4m^ , but the full function has not yet been 
computed. Unlike the latter two expansions, for which only a few terms are 
known, our knowledge of the expansion of Il^'^'{q'^) around g^ = is very 
impressive, after the work of Ref. [167j where 30 terms of this expansion 
were computed. 

Although the full vacuum polarization function n(g^) is Stieltjes, there 
is no reason why all the individual contributions n''''^'^''"-'(g^) should also 
have this property. Amusingly, however, this happens to be true both for 
n(0)(g2) and U^^^q"^) [MD [ISS \M\. The case of n(o)(g2) is trivial as it 
coincides with the full vacuum polarization n(q^) for as=0, and the case 
of n'^^(g^) was verified in [161j . As we will now show, this is no longer 
the case for W'^'^q'^) because its power series expansion around q^ = does 
not satisfy the Hankel determinants' condition Eq. ()3.17p which holds for a 
Stieltjes function. As a reminder of that condition, let f{z) be a Stieltjes 
function with a power expansion around z = as: 



/(^) = E/" 



(5.30) 



n=0 



Let D(m, n) be the determinant constructed with the Taylor coefficients /„ 



D{m,n) 



Jm Jm+1 

fm+1 fm+2 



Jm+n 
Jm+n+1 



Jm+n Jm+n+l • • • Jm+2r, 



(5.31) 



As f{z) is a Stieltjes function by construction, its coefficient must satisfy 
D{m, n) > 0, for ah m, n [801 [901 [Ml [92l [79]. 
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Figure 5.4: Number of decimal places reproduced by the rational approximation in 
Eq. {5.34\l to the function n'^^'(z) as a function of z, in the interval 0.4 < z < 0.9. 
In this plot we show, overlaid, the results for /3 = 0.1, 0.2, 0.3, ..., 1. As one can see, 
the particular value of (i chosen makes no significant impact on the result. 

The particular case of n'^''(g^), however, is different since using the /„ 
coefficients given in Ref. |167j (in the on-shell scheme, with the number of 
hght flavors n^ = 3): 



z"^n(2)(z) Si 0.631107 + 0.616294 z + 0.56596 z^ + 0.520623 z^ + 



one can immediately see that, e.g. 
^(0,1) 



0.631107 0.616294 
0.616294 0.56596 



-0.0226376 < 



(5.32) 



(5.33) 



This proves that the individual function 'ny^'{q^) is, all by itself, not a Stielt- 
jes function, even though the combination n(g^) in Eq. (j5.22p is. Therefore, 
we will now focus on applying the Theory of Fade Approximants to the 
full combination n(q^) in Eq. ()5.22p in order to extract information on the 
individual term n'^)(g^). 



5.5.2 The method 

Since, as it is obvious from Eq. ()5.22p . the function n(g^) depends on the 
value of Os, any Fade Approximant to it will also depend on the value of as, 
i.e. -P/Y^ {z;as)- This means that it is possible to construct a rational ap- 
proximation to the three functions Ii^^'^''^> [q^] from three different sequences 
of Fade Approximants to n(g^) constructed at three arbitrary values of a^, 
let us say a^ = 0, ±/?, with (3 sufficiently small so as to be able to neglect 
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the terms of ©(af) in Eq. ()5.22p . In this way one obtains 

--^nW(z) « ^{p^+^(.;a. = /3)-P^+^(z;a, = -/3)} 

where J > — 1 and N ^ oo. Since the value of /3 chosen is arbitrary, the 
N -^ oo hmit should produce results which are independent of /3, due to the 
convergence of the Pade Approximants to n(g^). Therefore, one should see 
that the three combinations ()5.34p are increasingly independent oi P as N 
grows This is indeed what happens. 

Furthermore, since we know the exact functions 11^'^^ (z) and Il^^'{z), we 
can compare them to the rational approximation on the right hand side 
of the first and second Eqs. ()5.34p in order to test the approximation. 
Actually, the case of Il^^){z) is rather trivial because is independent on a^ 
and also because is a direct application of the theorem of convergence. More 
interesting is the case of II^^' (z) where we can really see how the method 
works. Figure 15.41 shows the number of decimal places reproduced by this 
rational approximation in the interval 0.4 < z < 0.9, when A^ = 14 and 
J = (i.e. the diagonal Fade Approximant -P/f), for values of /3 in the 
interval 0.1 < /3 < 1. As one can see, the dependence on /3 cannot be 
distinguished in the Fig. 15.41 and the accuracy reaches, e.g., ~ 10 decimal 
places at z = 0.9. 

We remark that f3 ~ 0.1 corresponds to /? ~ as{Mz), a value which is 
clearly expected to be perturbative. Of course, the method can be used for 
smaller values of /3 as well but we use what we think are realistic numbers. 
For the values of f3 we choose, we have checked that all the determinantal 
conditions in Eq. (12), which the Taylor coefficients of a Stieltjes function 
must satisfy, are indeed satisfied by the 30 Taylor coefficients made available 
in Ref. |167j . As it turns out, there are 392 determinantal conditions for 
this number of coefficients. 

The third Eq. ()5.34p yields the desired approximation to 11^^^ (z). To be 
precise, it gives us a rational approximation to W'^'{z) in any compact set 
of the z complex plane, away from the cut 1 < z < oo. Since the threshold 



This independence of /3 in the case of fl' ' (g ) is trivially true 



120 



expansion at z ~ 1 can be written as |1661I169] : 

nj^^(z) = • + [0.34375 - 0.0208333 n^] ln2(l - z) 

+ [0.0116822 ne + 1.64058] ln(l - z) + K^^'^ 
+ [-0.721213 - 0.0972614 n^ + 3.05433 ln(l - z)] ^/Y^~z + 0(1 -(5k^5) 

in terms of an unknown constant K^'^'^ our Pade Approximation (|5.34p may 
be used to determine this constant, as we will next discuss. In this threshold 
expansion we take n^ = 3 as the number of light flavors. Even though the 
numerical coefficients have been rounded off for simplicity, they may be 
extracted exactly from the results in Ref. [169j . 

Since PAs are not convergent on the physical cut, it is impossible to 
match the rational approximants (j5.34p to the threshold expansion (j5.35p 
as a function of z. This fact is obvious from the presence of logarithms 
and squared roots in Eq. (|5.35p . However, both approximations (|5.34p and 
()5.35p are valid for values of z at a finite distance from the cut and, in 
particular, in a certain window in the interval < z < 1. Within this 
window, a numerical matching of ()5.34p and (|5.35p is possible and will in 
fact allow us to determine the unknown constant K^'^'. 

In order to determine this window, we make the following observations. 
First, although the rational approximation (j5.34p is convergent as A^ — t- oo 
in the interval < z < 1, it is more accurate the closer one gets to z = 
in this interval, for a given value of N . On the other hand, the threshold 
expansion (|5.35p is more accurate the closer one gets to the branching point 
at z = 1. From these two competing effects it is possible to determine an 
optimal window in z by minimizing a combined error function. We will call 
this error function £-(^z\ 

The function <?(z) has to take into account the error from the Pade 
Approximants as well as the error from the threshold expansion. To estimate 
the error from the PAs, we consider the difference between two consecutive 
elements in the same sequence, i.e. \P^^ — P;sjZ\^ I- As to the threshold 
expansion, we estimate its error as |1 — z|, since the expression ()5.35p is 
accurate up to terms of 0(1 — z). Therefore, in order to avoid possible 
accidental cancelations between the two errors, we define our combined error 
function as the following sum: 



E{z\- 



tt'^z 



^^^{PN^'i^'^s = (3) + Pi:;+'{z;as = -^) -2Pl:l+''{z;as = 0)] - {n^N-i} 
+ \l-z\ . (5.36) 



Minimizing £{z) with respect to z in the interval < z < 1, for every given 
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Figure 5.5: Predicted value for K^"^' from the sequence of diagonal, Pj^ (left panel), 
and first paradiagonal, Pj^ ~ (right panel), Pade Approximants. This figure corre- 
sponds to P = 0.5. 
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Figure 5. 6: Matching of the rational approximant in Eq. \5.34^ (solid blue line) to 
the threshold expansion II^j^ (z) in i5.35\) (solid-dashed red line) for N = 14, J — 
and the value of K^"^' in i5.3l^) . This figure shows the result for different values 
of /3 in the range < /3 < 1 , but the dependence on /3 is so small that cannot be 
discerned. 



values of N and /3, we may determine a value of z at the minimum, namely 
z* . This z* is then the one used to determine the constant K^'^' as 



i^{2) 



TT^Z* 
2/32 



dN+J 

N 



{z*;as 



fi) + P^+\z*;a, 



-/3) 



2P^^-\z* 



a 



s = o)}-n 

(5.37) 



;2(^*) 



for the given A^ and j3. In Eq. (|5.37p . nj:^^ (z*) stands for the expression 
in Eq. ()5.35p without the constant X'^) and, of course, without the term 
0(1 — z), evaluated at z = z* . The knowledge of 30 terms from the low- 
energy expansion gives us enough information to be able to construct up to 
the Fade Approximant Plf from the diagonal sequence, and up to the Fade 
Approximant P/5 from the first paradiagonal sequence. This corresponds 
to J = and J = — 1 in Eq. (|5.37p . The theorem of convergence for 
Stieltjes functions, sec l3.1.3l guarantees convergence for all the sequences 
P^"*" , J > — 1. We focus our effort on the sequences which have more 
elements, i.e., the diagonal and the paradiagonal once, to obtain the largest 
possible sequence. J = — 1 and J = translate into 14 elements for each 
sequence, while for bigger values of J, the number of elements decrease at the 
same rate as J increase, with the limit of J = 29 been the Taylor expansion. 
In all cases considered we have varied /3 in the generous range < /3 < 1, 
but our results are insensitive to this variation within errors, as expected. 
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5.5.3 Result and checks 

Our results for the constant K^"^' are shown in Fig. 15. 5[ This figure shows 
the convergence of the diagonal sequence and the first paradiagonal sequence 
as a function of the order in the PA. As one can see, we find a very nice 
convergence in both cases, with compatible results. Based on this analysis, 
we obtain the following value of K^'^': 

K(2) =3.71 ±0.03 . (5.38) 

This result is very close to, although slightly smaller than, the value obtained 
in Ref. |T66J, i.e. i^^^) ^ 3 g^ ^ o.02. 

The error bars shown on Fig. I5.5l have been calculated as ziz£{z*). Look- 
ing at the figure we see that the change in the value of K^"^' from one element 
of the sequence to the next is of the order of the errors shown, which is a 
good sign that the estimate for the error we have made is correct. In fact, 
we have explored what happens if the |1 — zj term in the error function 
in Eq. ()5.36p is multiplied by an arbitrary constant c and this constant 
is varied within a generous range, say, between c = 10 and c = 1/10. In 
this case the particular form shown in Eq. (|5.36p would correspond to the 
(natural) choice c = 1 (see Eq. ()5.35p ). While the central value for z* and, 
correspondingly, K^"^' are rather insensitive to these changes in c (within the 
errors quoted), the error function of course does change. However, one sees 
that for c = 10 the error would become much larger than the change in the 
central value between two subsequent PAs, whereas for c = 1/10 the error 
would become much smaller. The two cases would correspond either to an 
overestimation of errors (the former) or an underestimation of errors (the 
latter). We conclude, therefore, that our estimate of errors in Eq. ()5.38p is 
reasonable, and essentially the best we could make. Eq. ()5.37p seems rather 
sensitive to the particular value of z* . We have already said that variations 
of the parameter c do not cause a major impact in the determination of 
z* (which moves between the window 0.976-0.987 depending on c). What 
would be, however, the impact on z* of variations of < /3 < 1? For the PA 
P^^ that variation on /3 implies a variation on z* from 0.975 to 0.981, which 
translates to a variation of 5 x 10~^ of the value of K^"^' . 

Although we have taken symmetric errors for simplicity, it is also clear 
from the figure that the approach to the true value is made from below, 
so that a slightly more accurate determination could be achieved with the 
use of an asymmetric error. Apart from that, given the present knowledge 
of the expansions at low energy ()5.32p and at threshold ()5.35p . we find it 
difficult to believe any error estimate which could significantly go below our 
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Figure 5.7: Location of the poles in the Fade P^^ in the complex plane. Notice the 
accumulation of poles at z > 1, simulating the physical cut. 

figure in Eq. ()5.38p . Of course, should more terms in either expansion be 
known, a rerun of our analysis could immediately produce a more precise 
determination of K^"^' . 

Figure [5^ shows the matching of the rational approximant (i.e. the right 
hand side of the third of the Eqs. ()5.34p ) to the threshold expansion given 
by nj^^(z) in Eq. (lOHD . for A^ = 14 and J = 0, i.e. with the PA P^^, 
and for the value of K^'^' we have obtained. As one can see, this PA is 
able to reproduce, with high accuracy, the threshold expansion behavior in 
a window 0.92 < z < 1. Atz = l and above, the two lines in Fig. 15.61 will 
again diverge from each other, just as they do at low z. The value of z* 
minimizing the error function £{z) in (j5.36p was found at z* ~ 0.98 in this 
particular case. 

For illustration, in Fig. 15.81 we show the position of the poles in the 
PA P^^. As one can see, all the poles are sitting on the positive real axis 
above z = 1, as it should be. Notice how they accumulate in the region 
z > 1. This is how PAs approximate the physical cut present in the original 
function. As ensured from Pade Theory, this behavior was found in all the 
PAs considered. 

In order to further illustrate the reliability of our method we will include 
here several other examples. 

Let us consider the function G{z) in Eq. (j5.29p which, as z — )• 1, has the 
same behavior as the function n(^'(2:) we are interested in. One obtains 

G{z)\thr = ^^^=^+K(^^ + \7rVl^z-^il-z) + [{1 - zf^) , (5.39) 

where the constant K^ ' is of course exactly known, i.e. K^ ' = — 1. We 
find that the PA P^^ constructed from the Taylor expansion of G{z) yields 
a very accurate prediction for K^ ' , namely K^ ' = —0.999 it 0.005. The 
accuracy increases as the order of the Pade gets larger. For instance, the 
PA P|o° yields K^^') = -1.0001 ±0.0003. These error estimates are obtained 
on the basis of an error function of the form (j5.36p . but with |1 — z\^''^ 
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instead of |1 — z| in accord with the threshold expansion of G{z) above. 
And, of course, the result can be improved if instead of an error |1 — zp^ 
one uses |1 — zp. In this last case, with the PA P^^ the prediction for the 
constant K^"^) is K^'-^^ = -1.001 ± 0.001. For the PA Pj^?, the prediction 
reads K^'^^ = -1.00003 ± 0.00003. 

One can repeat the same exercise with the function n*-"-* (z) itself and get 
similar results. That function, at threshold, reads: 



3Vr^7 l-z 5(l-z)3/2 ,0) , (1 



.2 



Il('Hz)W = --^^^^ + ^ - ''';/' + K(^^ + '-^^^ , (5.40) 



Stt tt^ 16tt vr^ 



where /s:(°) = A = 0.0675475. 



In that case, the P/f (z) gives a prediction for K(°) = 0.0674 ± 0.0001 
and the PA P50 (^) gives a prediction for i^^^) = 0.06754 ± 0.00002. 

A further example is obtained with the function n*^^-* (z) which, at thresh- 
old, reads 



^^^Hz)\thr = 0.477465^1 -z + K«- 4 log(l-^) 

16 



+ (1-^) 



- log(l -z) + 0.354325 



O ((1 - zf')5,Al) 



where again the constant K^^> is known, i.e. K^^> = —0.314871. A 
rerun of our analysis with the PA P^^ gives in this case the prediction 
ii-Ci) = -0.31493 ± 0.00060, from the second expression in Eq. (fSTMD . In 
Fig. one can see the rate of convergence if the PAs P^ to the constant 
K^^' . Alternatively, one could also construct the PAs directly from the 
Taylor expansion of n'^^(z). Within the error quoted, the two numbers 
agree. These examples confirm that our estimation of errors is reasonable. 

Furthermore, since 30 Taylor coefficients of n'^^(z) are known ( |167j ). 
one can entertain oneself by constructing a PA of low enough order that does 
not require all of the 30 terms in its construction, in order to then predict 
the rest of the Taylor coefficients and compare with their exact result. For 
instance, with the ffist 21 Taylor coefficients we have constructed a PA P^q 
for Il^'^'{z) using the last expression in our Eq. (|5.34|) with A^ = 10 and 
J = which, upon re-expansion about z = 0, predicts the value for the 
coefficient of the term z" of Il^'^'{z) for n = 22,23, ...,30. Using the same 
choice of values for /3 as before, i.e. 0.1 < /? < 1, the order of magnitude 

is given in Table [5Tl As one can 



in the relative error, e = ^exad 

see, the results turn out to be very accurate 
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Figure 5.8: Sequence of approximants to the constant K^^' in Eq \5.41[ together with 
the exact result as a red line. 



n 


22 


23 


24 


25 


26 


27 


28 


29 


30 


e 


10-12 


10-11 


10-10 


10-9 


10-9 


10-8 


10-8 


10-7 


10-7 



Table 5.1: Order of magnitude in the relative error, defined as e = 
, of the prediction for the coefficient c„ of the term c„z" 



in the Taylor expansion of Il^'^'{z) from the PA P^q in Eq. (j5.34p . 
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Finally, as a further test of our method, we have calculated the value of 
the constants Hq and H^ which appear in the large- 2; expansion of the 
function 11^'^' (z) (we take n/ = n£ + 1 = 4 in the following expression): 



r(2) 



li-High-zi^) = (0.034829 - 0.0021109 n/) ln^(-4z) + (-0.050299 + 0.0029205 n/) ln(-4z) 

(2) , ,ln^(-4z) , ,ln(-4z) 
+ i7^ -• + (0.18048 - 0.0063326 n/) — ^ + (-0.59843 + 0.027441 nj)— - 



+ fl^+ o(z-Hn\-z)j . (5.42) 

Using our method, we find Hq = —0.582 ± 0.008. This result is to be 

compared to the true value Hh = —0.5857 |165j . If we now input this exact 

(2) 
value of Hk , by a rerun of the method, we may then determine the value of 

(2) (2) 

H^ . In this way, we find H) = —0.194 it 0.033, which is to be compared to 

the exact value H^ = —0.1872 |165j . Again, we find this agreement rather 

reassuring. 
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Chapter 6 



Pade Theory and 
Experimental Data 



6.1 Motivation 

The fourth chapter of the present work has extensively talked about the 
application of Pade Theory to meromorphic functions (but not Stieltjes func- 
tions) in the context of two-point Green's functions. We focused on mero- 
morphic functions due to the approach we used, the Large- A'^c limit. In that 
limit, Green's functions are meromorphic functions (sec. 12. 4p and then the 
Pommerenke's theorem can be applied. While the limit of large- A'^c is close 
to realistic [30], it is, in itself, an approximation. 

The fifth chapter has focused on Stieltjes functions and meromorphic 
functions of Stieltjes type. In both cases, a theorem of convergence for PA 
sequences is been applied to the low energy amplitudes in the Linear Sigma 
Model and for a Vector Meson Model, and also to the vacuum polarization 
of a heavy quark. And in both cases we obtain nice convergence patterns. 

In the present chapter, we would like to apply Pade Theory to a more 
realistic case, using now experimental data and nothing more, neither Large- 
Nc limit nor chiral limit. 

The particular case we would like to study is the pion vector form factor 
(defined in sec. I2.2.2p in the space-like region, which has been known for 
along time that is very well described by a monopole ansatz of the type 
given by Vector Meson Dominance (VMD) in terms of the p-meson (M^). 
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While since its first phenomenological evidence, some effort has been done 
to extended |170] . it has remained unclear whether there is a good reason 
for this from QCD or it is just a mere coincidence and, consequently, it is 
not known how to go about improving on this ansatz, [171j . [172j . 

6.1.1 The Pion Vector Form Factor 

To begin our discussion, let us define the vector form factor (VFF), F(Q^), 
by the matrix element 

(7r+(p')l I ui^u-^ drd-^ s^s \n+(j})) = {p + p'Y F{Q^) , (6.1) 

where Q^ = —{p'— p)^, such that Q^ > corresponds to space-like data. 
Since the spectral function for the corresponding dispersive integral for 
F{Q'^) starts at twice the pion mass, the form factor can be approximated 
by a Taylor expansion in powers of the momentum for IQ"^] < (2171^^)'^. At 
low momentum, Chiral Perturbation Theory is the best tool for organizing 
the pion interaction in a systematic expansion in powers of momenta and 
quark masses (sec. 12. 2p . This means that the coefficients in the Taylor ex- 
pansion can be expressed in terms of LECs and powers of the quark masses 
(sec. [2X2D . 

In principle, the coefficients in the Taylor expansion may be obtained h 
means of a polynomial fit to the experimental data in the space-like regio: 
below \Q'^\ = (2?7i7r)^- However, such a polynomial fit implies a tradeoff. 
Although, in order to decrease the (systematic) error of the truncated Taylor 
expansion, it is clearly better to go to a low-momentum region, this also 
downsizes the set of data points included in the fit which, in turn, increases 
the (statistical) error. In order to achieve a smaller statistical error one 
would have to include experimental data from higher energies, i.e. from 
IQ'^I > (2?Ti-7r)'^- Since the Taylor expansion can not be applied at that 
high Q^, one can not recover correctly its coefficients through a polynomial 
fit, then the use of alternative mathematical descriptions may be a better 
strategy. 

One such description, which includes time-like data as well, is based on 
the use of the Roy equations and Omnes dispersion relations. This is the 
avenue followed by |173l I174^ 1175) , which has already produced interesting 



by 



^Time-like data is provided by tttv production experiments and, consequently, they 
necessarily correspond to values of the momentum above the tttt cut, i.e. \Q^\ > 4m^ with 
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results on the scalar channel [176j . and which can also be applied to the 
vector channel. Other procedures have relied on conformal transformations 
for the joint analysis of both time-like and space-like data \n7\ I178J . or 
subtracted Omnes relations [1791 I180J . Further analyses may be found in 
Ref. [18n[T82] . 

On the other hand, as already mentioned above, one may also consider 
an ansatz of the type 



Q' Y 



\ VMD/ 

to describe the VFF. Even though the simplicity of the form of Eq. ()6.2p 
is quite astonishing, it reproduces the space-like data rather well, even for 
a range of momentum of the order of a few GeV, i.e. Q^ 2> Am^. If this 
fact is not merely a fluke, it could certainly be interesting to consider the 
form (j6.2p as the first step in a systematic approximation, which would then 
allow improvement on this VMD ansatz. 

The presence of a resonance such as Myj^^ means that the Taylor expan- 
sion of F{Q'^)vMD has a radius of convergence below \Myj^jj\. For values 
of IQ^I > \Myj^jj^\ this representation will diverge. We would like to point 
out that the pion vector form factor is a good laboratory to apply the Fade 
method in front of the usual Taylor expansion. PA allow the inclusion of low, 
medium and high energy information in a rather simple way which will let 
us cover a larger range of experimental data that with a simple polynomial. 

Also, the VMD ansatz can be viewed as the first element in a sequence 
of Fade Approximants which can be constructed in a systematic way. By 
considering higher-order FA in the sequence, one may be able to describe 
the space-like data with an increasing level of accuracy □. 

In the present chapter we will show that, due to the precision allowed 
by the experimental data, there are sequences of FAs which improve on 
the lowest order VMD result in a rather systematic way. This has allowed 
us to extract the values of the lowest-order coefficients of the low-energy 
expansion with an associated error. 

However, precisely these lowest-order coefficients of the low-energy ex- 
pansion that we want to extract is what one needs as an input in doing FAs 
(see sec. 14. ip . Then, instead of the typical procedure, one can consider to 



^Obviously, unlike the space-like data, one should not expect to reproduce the time- 
like data since a Pade Approximant contains only isolated poles and cannot reproduce a 
time-like cut. 
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apply a A^— point Fade using each experimental data. This means to build 
a rational function of degree L and M where L + M + 1 = N , being A^ the 
amount of available data. Nevertheless, this leads to a theoretical uncer- 
tainty. What is the best PA sequence to apply? The parameter L runs from 
L = OtoL = A^ — 1 while the parameter M runs from M = A^ — 1 to M = 0. 
We then have A^ different possible configurations and no way to improve the 
result after choosing a particular one. Also, as one can easily notice, this 
implies a large computation due to the amount of available data. Instead 
of this procedure, our strategy will consist in determining these coefficients 
by a least-squares fit of a Pade Approximant to the vector form factor data 
in the space-like region. This will let as improve the approximation and 
estimate the systematic error. 

There are also several types and sequences of PAs that may be considered 
in a fit procedure. In order to achieve a fast numerical convergence, the 
choice of which one to use is largely determined by the analytic properties of 
the function to be approximated (see, for example, the discussion in section 
14. 1|) . In this regard, a glance at the time-like data of the pion form factor 
makes it obvious that the form factor is clearly dominated by the p-meson 
contribution. The effect of higher resonance states, although present, is 
much more suppressed. In these circumstances the natural choice is a Pi 
Pade sequence [SOl [79], i.e. the ratio of a polynomial of degree L over a 
polynomial of degree one u. Notice that, from this perspective, the VMD 
ansatz in ()6.2p is nothing but the Pi Pade Approximant. 

In that point we may say that strictly speaking there are no mathemat- 
ical theorem that guarantees convergence of the sequence Pi{Q^) to the 
function F{Q'^) with the fit procedure. However, a hint to understand why 
our sequence will converge may be as follows. Since Tp/Mp <C 1 one can 
thing of F(Q'^) to be quasi-meromorphic function with only one simple pole 
at Mp. We have already shown how to deal with meromorphic functions. 
Had this the case, we would have applied the Montessus^s theorem and the 
convergence would have been guaranteed when L — )• oo. Since this is not 
the case, this theorem can not be directly applied and only remains as a 
mathematical inspiration to understand the results. 

To test then the aforementioned single-pole dominance, one should check 
the degree to which the contribution from resonances other than the p-meson 
may be neglected. Consequently, we have also considered the sequence P2, 
and the results confirm those found with the PAs Pi- Furthermore, for 



^Conventionally, without loss of generality, the polynomial in the denominator is nor- 
malized to unity at the origin. 



132 



completeness, we have also considered Fade- Type approximants (PTs) [1681 
I87j . with the feel that one can trust the single-pole dominance of the p- 
resonance. Finally, we have also considered an intermediate case, the Partial- 
Pade approximants (PPs) |1681 [87] , We have fitted all these versions of 
rational approximants to all the available pion VFF space-like data |1831 
[IHllIMllIMllIHIllIHHllIHHllIgOllIIIllinS]. The result of the fit, table[621 is 
rather independent of the kind of rational approximant sequence used and 
all the results show consistency among themselves. 



6.2 Evaluation of the systematic error 

In the previous section we have presented our strategy of fitting to the 
experimental data. What is new in that case is the way to estimate the 
systematic error of the approximation when dealing with real data. In the 
present section we will illustrate the usefulness of the PAs as fitting functions 
in the way we propose here, by using, first, a phenomenological model as 
a theoretical laboratory. Furthermore, the model will also give us an idea 
about the size of possible systematic uncertainties. 

To consider the method as quite realistic, we use a once-subtracted 
Omnes relation to recover the form- factor, 

F(g^)^expr^£.*A}. (,.3, 




The 5{t) function 



5{t) = tan 






(6.4) 



plays the role of the vector form factor phase-shift |1791 11801 1193^ 1194] 
and its t-dependence comes from the width (cr(t) = ^Jl — Arri^/t): 

f„(t)=ro f^l ^^eit-Arhl). (6.5) 

To emphasize that the model should be considered realistic, we chose the 
input parameters close the physical ones: 

Fo = 0.15 GeV , Mp^ = 0.6 GeV^ , Ami = 0-1 GeV^ . (6.6) 

Actually, this model has been used in Ref. [1791 11801 1193^ I194| to extract 
the values for the physical mass and width of the p meson through a direct 
fit to the (timelike) experimental data. 
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po 


PI 


P? 


Pf 


Pf 


Pf 


Pf 


F(Q2) (exact) 


oi(GeV-2) 
a2(GeV-4) 

a3(GeV-6) 


1.549 
2.399 
3.717 


1.615 
2.679 

4AU 


1.639 
2.809 
4.823 


1.651 
2.892 
5.097 


1.660 
2.967 
5.368 


1.665 
3.020 

5.579 


1.670 
3.074 
5.817 


1.685 
3.331 

7.898 


5p(GeV2) 


0.646 


0.603 


0.582 


0.567 


0.552 


0.540 


0.526 


0.6 



Table 6.1: Results of the various fits to the form factor F{Q'^) in the model, 
Eq. (j6.3|) . The exact values for the coefficients ai in Eq. (|6.7p are given on 
the last column. The last row shows the predictions for the corresponding 
pole for each Fade (sp), to be compared to the true mass Af^ = 0.6 GeV^ 
in the model. 



Expanding F{Q'^) in Eq. 
expansion for the VFF 



in powers of Q^ we readily obtain a Taylor 



F{Q' 



I - aiQ^ + a2Q^ - as Q^ + ... , 



(6.7) 



with known values for the coefficients ao. The goal is to obtain these a, 
parameters trough expanding the rational functions after fitting to the ex- 
perimental data. 

The real experimental data are obtained from |183t I184t I185t I186[ I187t 
I188t I189t I190t I19H 1192] and to work close to the real case, we will simulate 
the same amount of data as in the physical casqj. Then, we generate fifty 
"data" points in the region 0.01 < Q^ < 0.25, thirty data points in the 
interval 0.25 < Q^ < 3, and seven points for 3 < Q^ < 10 (all these momenta 
in units of GeV^) 

In Table [6TT] we show the predicted coefficients Oj from the Pi{Q'^) se- 
quence we have fitted. The last PA we have obtained to these data is Pf . 
As expected, the pole position of the Fades differs from the true mass of the 
model Eq. (j6.6p . since the Montessus de Ballore^s theorem is not fulfilled. 

A quick look at Table 16.11 shows that the sequence seems to converge 
to the exact result, although in a hierarchical way, i.e. much faster for ai 



*ao = 1 to resemble the vector current conservation condition. 

^To estimate the systematic error derived purely from our approximate description of 

the form factor, these points will be taken with vanishing error bars. 
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than for 02, and this one much faster than 03, etc, showing exactly the 
same convergence behavior than in the case worked in section I4.4.11 The 
relative error achieved in determining the coefficients a, by the last Pade, 
Pf, is respectively 0.9%, 8% and 26% for 01,02 and 03. Naively, one would 
expect these results to improve as the resonance width decreases since the P^ 
contains only a simple pole, and this is indeed what happens. Repeating this 
exercise with the model, but with a Fq = 0.015 GeV (10 times smaller than 
the previous one) , the relative error achieved by Pf for the same coefficients 
as before is 0.12%, 1.1% and 4.7%. On the other hand, a model with Fq five 
times bigger than the first one produces, respectively, differences of 2.1%, 
14.4% and 37.8%. 

To complete the study, we have also studied the convergence of Pade- 
Type approximants with the model. Thus, in this case, we have placed the 
T/" pole at Sp = Mp and found a similar pattern as in Table [HTTl For Tf, 
the Pade- Type coefficient oi differs a 2.5% from its exact value, 02 by 16% 
and 03 by 40%. 

This procedure gives us a rough estimate of the systematic uncertainties 
when fitting to experimental data. Since, as we will see, the best fit to the 
experimental data comes from the Pade Pf, we will take the error in Table 
16.11 from this Pade as a reasonable estimate and add to the final error an 
extra systematic uncertainty of 1.5% and 10% for ai and 02 (respectively). 



6.3 Fitting the pion Vector Form Factor 

We will use all the available experimental data in the space-like region, 
which may be found in Refs. [Ml [Ml [Ml [Ml [Ml [Ml [Ml [Ml [Ml [M] • 
These data range in momentum from Q^ = 0.015 up to 10 GeV^. 

As discussed previously, the prominent role of the /9-meson contribution 
motivates that we start with the P^ Pade sequence. 



6.3.1 Pade Approximants P/" 

Without any loss of generality, a P^ Pade Approximant is given by 

P[{Q') = 1 + X;«fc(-Q')' + (-Qy .^^^2 ' (6-8) 

fc=l ^+ ai ^ 

where -P^ (0) = 1 has been imposed (in fact, the vector current conservation 
condition) and the coefficients a^ are the low-energy coefficients of the cor- 
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Q^FiQ^) (GeV") 

0.8 




Q\Ge\) 



Figure 6.1: The sequence of P/" PAs is compared to the available space- 
like data nSa dSl ISSl lESl [Ml [Ml [Ml [Ml [Ml [M]: P? (green 
thick-dashed), P^ (brown dashed), P^ (orange dot-dashed), Pf (blue long- 
dashed), Pf (red sohd). 
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responding Taylor expansion of the VFF (compare with ()6.7p for the case of 
a model). 

The fit of P[ to the space-like data points in Refs. ^SS\ [IH [Ml [Ml 
[T871 [THHl [T89l [T90l [TOTl [T92] determines the coefficients a^ that best inter- 
polate them. According to Ref. |1951 11961 I197J . the form factor is supposed 
to fall off as 1/Q'^ (up to logarithms) at large values of Q^. This means that, 
for any value of L, one may expect to obtain a good fit only up to a finite 
value of Q^, but not for asymptotically large momentum. This is clearly 
seen in Fig. 16.11 where the Fade sequence P^" is compared to the data up 
to L = 4. 

Fig. 16.21 shows the evolution of the fit results for the Taylor coefficients 
ai and 02 for the P^ PA from L = up to L = 4. As one can see, after 
a few Fades these coefficients become stable. Since the experimental data 
have non zero error it is only possible to fit a P/' FA up to a certain value 
for L. From this order on, the large error bars in the highest coefficient in 
the numerator polynomial make it compatible with zero and, therefore, it 
no longer makes sense to talk about a new element in the sequence. For 
the data in Refs. [1551 [TM [TB51 fTMl flFTl flBSl \Tm [TMl [TUTl \m2\ . this 
happened at L = 4 and this is why our plots stop at this value. Therefore, 
from the PA Pf we obtain our best fit and, upon expansion around Q^ = 0, 
this yields 

ai = 1.92 ± 0.03 GeV~2 , 02 = 3.49 ± 0.26 GeV"^ ; (6.9) 

with a xVdof = 117/90. 

Eq. (|6.8|) shows that the pole of each P^ PA is determined by the ratio 
Sp = ai/aL+i- This ratio is shown in Fig. 16.31 together with a gray band 
whose width is given by zbMpFp for comparison. From this figure one can 
see that the position of the pole of the PA is close to the physical value 
-^p [S]) although it does not necessarily agree with it, as we already saw in 
the model of the previous section. 

6.3.2 Comment on P2 Fades 

Although the time-like data of the pion form factor is clearly dominated 
by the p(770) contribution, consideration of two-pole P^ PAs will give us a 
way to assess any possible systematic bias in our previous analysis, which 
was limited to only single-pole PAs. 

We have found that the results of the fits of P2 PAs to the data tend to 
reproduce the VMD pattern found for the P^ PAs in the previous section. 
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fli (GeV"^) 



2.00 



1.95 



1.90 
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Figure 6.2: ai and 02 Taylor coefficients for the P^ PA sequence. 

The P2 PAs place the first of the two poles around the /o-mass, while the 
second wanders around the complex momentum plane together with a close- 
by zero in the numerator, what we have already explained to be a defect 
[8OI I198J . sec. 13.1.31 As a local perturbation, at any finite distance from it 
the effect is essentially negligible. This has the net effect that the P2 Pade 
in the Euclidean region looks just like a P(" approximant and, therefore, 
yields essentially the same results. For example, for the P|, one gets 



ai = 1.924 ±0.029 GeV" 
with a xVdof = 120/92. 



02 = 3.50 ± 0.14 GeV" 



(6.10) 



6.3.3 Pade Type and Partial Pade Approximants 

As in the case of the model, one may consider other kinds of rational 
approximants, such as Pade Type and Partial Pade Approximants. Since 
the value of the physical p-mass is known {Mp = 775.5 MeV), it is natural to 
attempt a fit of PTAs to the data with a pole fixed at that mass. Actually, 
the results found with the model and also whit Pf and P| inspire this way 
to proceed. The corresponding sequence will be called T/". This has the 
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s„ (GeV") 



0.7 ^ 
0.6 _ 
0.5 
0.4 



Figure 6.3: Position Sp of the pole for the different P^. The range with the 
physical values Mp it MpTp is shown (gray band) for comparison. 

fli (GeV"^) 

1.95 



1.90 



1.85 



1.80 



Figure 6.4: Low-energy coefficient ai from the Tf" Fade- Type sequence 

obvious advantage that the number of parameters in the fit decreases by 
one and allows one to go a little further in the sequence. Our best value 
is then given by the Fade Type Approximant Tf, whose expansion around 
Q^ = yields the following values for the Taylor coefficients: 



ai = 1.90 ± 0.03 GeV-2 , 



02 = 3.28 ± 0.09 GeV" 



(6.11) 



with a xVdof = 118/90. 

The previous analysis of FTAs may be extended by making further use of 
our knowledge of the vector spectroscopy [5]. For instance, by taking Mp = 
775.5 MeV, Mp> = 1459 MeV and Mp// = 1720 MevH we may construct 
further Fade- Type sequences of the form T2 and Tg . 

In the FTA sequence T2 one needs to provide the value of two poles. For 



As will be seen, results do not depend on the precise value chosen for these masses. 
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fli (GeV~^) 
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1.95 



1.90 



1.85 



ai (GeV"^) 
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1.95 



1.90 
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Figure 6.5: Low energy coefficient oi for the T2 Pade-Type sequence with Alp and 
Alp' (left), and with Alp and AIp^ (right). 



the first pole, the natural choice is M^. For the second pole, we found that 
choosing either M^, or M^„ (the second vector excitation) does not make 
any difference. Both outcomes are compared in Fig. (j6.5p . Using M^,, we 
found that the T2 PTA yields the best values as 



ai = 1.902 ± 0.024 GeV-^ , 



a2 = 3.29 ± 0.07 GeV""^ , (6.12) 



with a xVdof = 118/92. 

Using M^„ as the second pole one also gets the best value from the T^ 
PTA, with the following results: 



ai = 1.899 ±0.023 GeV'^ , 



02 = 3.27 ±0.06 GeV"'', (6.13) 



with a x^/dof = 119/92. We find the stability of the results for the coeffi- 
cients ai^2 quite reassuring. 

We have also performed an analysis of the PTA sequence T3 , with similar 
conclusions. From the T3 we obtain the following values for the coefficients: 



ai = 1.904 ± 0.023 GeV'^ , 



a2 = 3.29 ± 0.09 GeV""^ , (6.14) 
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Fade 


ai (xGey2) 


aa (xGey^) 


Pf 


1.92 ±0.03 


3.49 ±0.26 


Pi 


1.924 ±0.029 


3.50 ±0.14 


Tf 


1.90 ±0.03 


3.28 ±0.09 


Ti 


1.902 ± 0.024 


3.29 ±0.07 


n 


1.899 ± 0.023 


3.27 ±0.06 


n 


1.904 ±0.023 


3.29 ±0.09 


Ki 


1.902 ± 0.029 


3.28 ±0.09 


final 


1.91 ±0.01 


3.30 ± 0.03 



Table 6.2: Summary of results found for the parameters ai and 02 defined 
in Eq. ()6.7p for all the PAs we have worked. See the text for details on the 
notation. 

with a xVdof = 119/92. 

Finally, to complete our analysis, we will also consider Partial Pade 
Approximants, in which only part of the denominator is given in advance. 
In particular, we study the PPA sequence P^^ 13 in which the first pole is 
given by M^ and the other is left free. The best determination of the Taylor 
coefficients is given by -Pf 1, and they yield 

ai = 1.902 ± 0.029 GeV-^ , 02 = 3.28 ± 0.09 GeV""^ , 



(6.15) 



with the free pole of the PPA given by M'j^^^ = (1.6 ± 0.4 GeV)^ and a 
xVdof = 119/92. 



6.4 Combined Results 

The combination of all the previous rational approximants results in an 
average given by 

ai = 1.907±0.010stat±0.03syst GeV'^ , 02 = 3.30±0.03stat±0.33syst GeV"^ 



(6.16) 



'^See sec. 13. 1.41 for notatfon. 
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The first error conies from combining the results of the different fits by means 
of a weighted average. On top of that, we have added what we beheve to 
be a conservative estimate of the theoretical (i.e. systematic) error based 
on the analysis of the VFF model in Sec. 16.21 We expect the latter to give 
an estimate for the systematic uncertainty due to the approximation of the 
physical form factor with rational functions. For comparison with previous 
analysis, we also provide in Table [U^ the value of the quadratic radius, which 
is given by (r^) = 6ai , that would mean Lg = (6.80 ± 0.04) x 10~^, |171| . 
Since the value that we obtain for the quadratic radius is slightly smaller 
than Amendolia'86 result ([183]), our result for Lg is also smaller. 

We also provide in Table 16.31 other determinations of ai and 02 (this sec- 
ond one when available) found in the literature. Refs. [1731 11741 [T75] . named 
CGL, are based on the use if the Roy equations and Omnes dispersion rela- 
tions. Refs. |177lll78j . named TY, rely on conformal transformations for the 
joint analysis of both time-like and space-like data. Ref.[1991, called BCT, 
computed the 0{p^) Vector Form Factor in Chiral Perturbation Theory and 
then they fit to experimental data to extract the LECs. Finally, Ref. [18U] . 
named PP, used a model independent parametrization for the Vector Form 
Factor constrained by unitarity and analyticity. We also include the PDG 
[5] value based on a weighted average of some available experimental data 
such as [183], [EZ], [192] . 

Recently, two lattice collaboration were able to compute on the lattice 
the Pion Vector Form Factor. In Ref.[2D0], the RBC-UKQCD collabora- 
tion study that Form Factor in a gauge configurations with Domain Wall 
Fermions and the Iwasaki gauge action. In Ref. [201] . the JLQCD-TWQCD 
collaboration used two-flavor lattice QCD with exact chiral symmetry and 
employing overlap quark action. We would like to remark that this last 
Ref. [201j found that the best fit to their lattice data could we obtain using 
a VMD-like hypothesis: 



P{Q') = a^{-Q^) + a^i-Q^f + asi-Q'f + ^-^ , (6.17) 

1 + ^ 

where M^ is the value of the p-meson on the lattice. It is really interest- 
ing to notice how close this expression Eq. (|6.17p is to our Eq. (|6.8p for L = 4, 
which is, actually, our best P^ approximant to the experimental data. 
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(r2) (fm2) 


aa (GeV-4) 


PDG \hj 


0.452 ±0.016 




This work falso |171p 


0.445 ± 0.002stat ± 0.007syst 


3.30 ± 0.03stat ± 0.33syst 


cGLimiiiziiins] 


0.435 ± 0.005 




TY 1177111781 


0.432 ± 0.001 


3.84 ±0.02 


BCT tl99] 


0.437 ±0.016 


3.85 ± 0.60 


PP[I80] 


0.430 ± 0.012 


3.79 ± 0.04 


RBC-UKQCD \m)\ 


0.418 ±0.031 




JLQCD-TWQCD [201J 


0.409 ± 0.044 


3.22 ±0.42 



Table 6.3: Our results for the quadratic radius (r^) and second derivative 
02 are compared to other determinations [I73l HZl HTSl EZl EHl HSU fTMl 
I200tl201j . Our first error is statistical. The second one is systematic, based 
on the analysis of the VFF model in section 16. 2i We also include the PDG 
[5] value based on a weighted average of some available experimental data 
such as [183], [W], [T92] . 
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Chapter 7 



Conclusions 



In this work we have exploited several uses of Pade Approximants as 
approximations to QCD. In particular in the chapter H] we saw that the 
resonance saturation in the limit of large- A'^c rnay be reinterpreted within the 
mathematical Theory of Pade Approximants to meromorphic functions. We 
saw that one may expect convergence of a sequence of Pade Approximants 
to QCD Green's functions in the large- Ac limit in any compact region of the 
complex Q^ plane except in a region of zero-measure. As the order of the PA 
grows, the convergence properties guarantees that any given artificial pole 
of the approximation goes to infinity or is almost canceled by a nearby zero. 
With the help of a {VV — AA) model we have reviewed the main results 
of this theory where we have seen that while in the Minkowsky region the 
rational approximant creates the expected artificial poles, in the Euclidean 
region the description is quite accurate. And this happens hierarchically: 
although the first poles and residues in a PA may be used to describe the 
physical mass and decay constants reasonably well, the same is not true 
for the last ones. In general it will be unreliable to extract properties of 
individual mesons from an approximation to large-Ac QCD with only a 
finite number of states. 

Studying form factors one can find particular manifestation of that prob- 
lem. Since a form factor, like a decay constant, is obtained as the residue 
of a Green's function at the corresponding pole(s), this also means that 
one may not extract a meson form factor from a rational approximant to 
a 3-point Green's function, in agreement with Ref . |55j . This observation 
may explain why the analysis of Ref. [49j . which is based on an extraction 
of matrix elements such as (7r|S'|P) and {'k\P\S) from the 3-point function 
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{SPP), finds values for the Ki^ form factor which are different from those 
obtained in other analyses |202j or lattice results |203j . 

In spite of all the above problems related to the Minkowski region, our 
model shows how Pade Approximants may nevertheless be a useful tool in 
other regions of momentum space. We think that this is also true in the 
real case of QCD in the large-A'^c limit. In this case one may use the first 
few terms of the chiral and operator product expansions of a given Green's 
function to construct a Pade Approximant which should yield a reasonable 
description of this function in those regions of momentum space which are 
free of poles. In this construction, Pade Approximants containing complex 
poles, if they appear, should not be dismissed. We exemplify this procedure 
with the estimate of the LECs at 0{p^) of the correlator I[lii{Q^). 

However, if not all the residues and masses in a rational approximant 
are physical, this poses a challenge to any attempt to use a Lagrangian with 
a finite number of resonances such as, for example, the ones in Ref. |4HI42j. 
for describing Green's functions in the large- A'^c limit of QCD. Even if these 
Lagrangians are interpreted in terms of PTAs, with the poles fixed at the 
physical value of the meson masses, we have seen how the residues then get 
very large corrections with respect to their physical counterparts. 

Besides, we also addressed on the different determinations of the mass 
of first axial state mA_ using resonance saturation in the large- Ac limit of 
QCD. We show how these different determinations may be understood as 
the position of the pole of a Pade Approximant and then, nothing to do 
with the physical lowest axial state appearing in the Particle Data Book 
[5]. We saw that PAs and PTs leads basically to the same results when the 
position of their poles may be completely different, even complex for the 
case of the lowest PA (depending of the value of the combination Q. In 
passing we also suggested a way to disentangle the discussion about the sign 
of the condensate of 0{p~^) and its relation with the low-energy constant 
Lio appearing in the Chiral Lagrangian. 

We divided the chapter 5 in two parts to address the issue about Stieltjes 
functions and Pade Approximants. Firstly we tested the reliability of the 
unitarization of low-energy amplitudes in the context of tha Linear Sigma 
Model. We saw that the predictions of the sigma mass and width in the 
LSM with the lAM method are too much out from the original ones. We 
exploited the possibility of lAM being a Pade Approximant sequence which 
lie out of the application of the theorem of convergence and then leading to 
inappropriate results. On the contrary, the sequence P^ displays a quick 
convergence behavior with reliable results for the sigma pole and width. 
For all this, we suggested the use of the P^ sequence rather than P^. 
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Unfortunately, the study on broad resonances requires to go beyond the 
tree-level approximation. Thus, in the real world, the IJ = 00 channel 
needs the inclusion of loops, precluding by now the extension to PAs beyond 
0{p^) in the chiral expansion, i.e. Pi- There is also a clear limitation on 
our experimental knowledge of the low-energy couplings, which barely goes 
beyond 0{p^). However, our proposal should be still suitable for the analysis 
of theories with narrow resonances and a relatively good knowledge of the 
experimental low-energy amplitudes. 

In a second part of the chapter we discussed the extraction of parame- 
ters from the vacuum polarization function of a heavy quark through Pade 
Approximants. We presented a clean way to handle the possible ambigu- 
ities of an approximation to the vacuum polarization function of a heavy 
quark at 0{a1). Since the vacuum polarization function, when considered 
in terms of the external momentum g^, is a function of the Stieltjes type, 
the theorem of convergence guarantees, for example, that the poles of your 
approximation will lie on the real axis mimicking the branch cut present 
on the original function. Also, with the amount of information obtained 
on the Taylor expansion around g^ = for that function we estimate the 
value for the threshold parameter K^"^' which can not be yet computed with 
Feynmann diagrams. 

Finally, in chapter [6] we have explored the reliability of fitting the pion 
vector form factor with rational approximations. Because these approxi- 
mants are capable of describing the region of large momentum, they may 
be better suited than polynomials for a description of the space-like data. 
As our results in Table 16.31 show, the errors achieved with these approxi- 
mants are competitive with previous analysis existing in the literature. As 
a straightforward method, can be easily applied for analyzing other form 
factors in the spacelike region and also to analyze lattice data which usually 
covers large energy region. 
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Appendix A 



Dispersion relations 



A dispersion relation is as has been shown by Kahen and Lehmann [2U4j . 
|205j quite a long time ago a relation that two-point functions have to obey. 
That dispersion relation follows from the analyticity properties of n(q'^) as a 
complex function of g^, the only energy- momentum invariant which appears 
in a two-point function. This tool will be interesting in our calculation of 
the {VV — AA) correlation function in chapter [H In full generality n(q^) is 
an analytic function in the complex g^-plane but for a cut in the real axis 
< q^ < oo. It then follows that ( |206) for a demonstration): 



1 . 2 



n(g^) = / dt 5 lmn{t)+a + bq^ + ■■■ (A.l) 

Jo t-q^ - zevr vr 

where the degree of the arbitrary polynomial in the r.h.s. depends on 
the convergence properties of Imn(t) when t — )• oo. The interest of this 
representation is that -lmn(t) in the integrand, which is usually called 
spectral function p{t), is a physical cross-section. For example, with 

J/^(x) = -u{x)-i^'u{x) - -d{x)j''d{x) - -s(x)7^s(x) (A.2) 

o o o 

the electromagnetic hadronic current of light quarks, the relation to the 
total e~^e~ annihilation cross-section into hadrons is 



<7(9 )e+e-^hadrons = —J2~^ -lmI\.EM{q ) (A.3) 



.2^ _ '*_^^ 2 

and 
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(A.4) 
where the sum is extended to all possible physical states, on-shell states, 

with an integration over their corresponding phase space understood. 

The physical meaning of the coefficients of the arbitrary polynomial in 

the r.h.s. of the Eq. (jA.ip depends on the choice of the local operator J{x) 

in the two-point function. In some cases the coefficients in question are fixed 

by low-energy theorems. 

If n(0) is known, we can trade the constant a in Eq. (jA.ip : 

n(g2) = n(0)+/ ^- ^— Imn(t) + 6g2 + . . . (A.S) 

Jo t t — q"^ — leir vr 

In other cases the constants are absorbed by renormalization constants. 
In general, it is always possible to get rid of the polynomial terms by taking 
an appropriate number of derivatives with respect to q'^, which is a substrac- 
tion process. Eq. ()A.5p . due to the derivative, is called one-substraction 
relation. 
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Appendix B 

Weinberg Sum Rules (WSR) 

In the chiral limit, the first Weinberg sum rule [207] 

/ dtlmULRit) = , (B.l) 

Jo 

implies that 

^/2m2-^/1m1 = /2; (B.2) 

V A 

and the second Weinberg sum rule )207j 

/•oo 

/ dttImIiLR{t) = 0, (B.3) 

Jo 

furthermore implies that 

J2fvK-EfAMl = 0. (B.4) 

V A 

Assuming enough convergence in the ressonance region, in QCD|j, Eq. (jB.2p 
follows from the fact that there is no local order parameter of dimension 
d = 2 and Eq. ()B.4p from the absence of a local order parameter of dimen- 
sion d = A. The first possibly non trivial contribution comes from local order 
parameters of dimension d = 6, as shown in Eq. (j4.5p . 



^For a discussion of the Weinberg sum rules in QCD in the presence of explicit chiral 
symmetry breaking see ref. |208j 
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